Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,1050AG39.2% 105.3 / 27.4 74I441V (ATC→GTC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (8/37);  new base G (29/0);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 2.47e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGC  >  NC_000913/4456958‑4457252
                                                                                                                                                   |                                                                                                                                                     
atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATc                                                                                                                                                     <  2:317984/150‑1 (MQ=255)
atTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATc                                                                                                                                                     <  1:676458/150‑1 (MQ=255)
     gtggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAg                                                                                                                                               <  1:470031/151‑1 (MQ=255)
        gcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGaaa                                                                                                                                            <  2:1244116/151‑1 (MQ=255)
              ggcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActgctg                                                                                                                                      <  1:536197/151‑1 (MQ=255)
               gcagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTgg                                                                                                                                     <  1:231985/151‑1 (MQ=255)
                  gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATg                                                                                                                                  <  2:1143580/151‑1 (MQ=255)
                    agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGAGCCATCCGAAACAGCtgggtggg                                                                                                                                >  2:391587/1‑151 (MQ=255)
                      tggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCACCAGAAACTGCTGGCTggg                                                                                                                                >  1:1068834/1‑149 (MQ=255)
                            aagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAAAAGCTGggtggggtgg                                                                                                                            >  2:714621/1‑147 (MQ=255)
                               ccTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaagaa                                                                                                                     <  1:1135022/151‑1 (MQ=255)
                                cTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTCTTTGGGGGGGCCCTCACAAACCGCTCGGTGGGGGGGCggagga                                                                                                                     >  1:906189/1‑148 (MQ=255)
                                  gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGaagaag                                                                                                                    <  1:174863/149‑1 (MQ=255)
                                     ttCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCAGAAACTGCTGGATGGGCTGGCGGAGAAAGCgg                                                                                                                >  2:593242/1‑150 (MQ=255)
                                      tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTGCCTCAGAAAAAGCTTGGTGGGGGGGCGGAAAAGGCgg                                                                                                                >  1:692590/1‑149 (MQ=255)
                                       cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAg                                                                                                             <  2:654010/151‑1 (MQ=255)
                                             gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCcgc                                                                                                         <  2:895725/149‑1 (MQ=255)
                                              cacaCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCATCAGAAAATGCTGGGTGGGGTGGCgaagaaagggg                                                                                                                >  1:1120908/1‑136 (MQ=255)
                                              cacaCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCa                                                                                                      <  1:767641/151‑1 (MQ=255)
                                                  cTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCCGAAAATGCTGGAGGGGGTGGCGGAGAAGGGGGGAGGCGCGCAgga                                                                                                   >  1:1247265/1‑148 (MQ=255)
                                                   tGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGGTGGGGGAGAAAGGGGGAGGCGCGCAGGaa                                                                                                  >  2:177034/1‑150 (MQ=255)
                                                    ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAAtt                                                                                                <  2:520445/151‑1 (MQ=255)
                                                     gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTaa                                                                                                  >  2:1019491/1‑148 (MQ=255)
                                                     gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                               <  2:74857/151‑1 (MQ=255)
                                                     gAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTc                                                                                               <  2:933362/151‑1 (MQ=255)
                                                                tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCCGAAACTGCTGGGTGGGGCGGGGGAGAAAGCGGGAGCCGCGCCGGAAATAGGGCTTAGcc                                                                                    >  1:880667/1‑150 (MQ=255)
                                                                tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGTCCCTCCGAAACTGCTGGGTGGGCTGGCGAAGAAAGCGGGAGCCGCGCCGTAAATTGGCCCTTGccc                                                                                   >  2:145377/1‑150 (MQ=255)
                                                                 ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCATCAGAAAATGCTGGAAGGGCTCGCGAAAAAAACGGGAGCCGCGCAGTAATTCCGCCTTAGACa                                                                                   >  2:1221829/1‑150 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGGGGAGAAAGGGGGAGGCGCGCAGTAAATCGGCCTTAGcc                                                                                    >  1:76012/1‑147 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCTCAAAAACAGCTGGGATGGCTGGGGGAAAAAGGGGGAGCCCCGCAGGAATTTGGccctttcc                                                                                     >  1:325181/1‑141 (MQ=255)
                                                                    tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGGGGGAGGCGCGCCGGAATTTCGCCTTTGCCAg                                                                                  >  1:922164/1‑148 (MQ=255)
                                                                     aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGatat                                                                              <  1:1029729/151‑1 (MQ=255)
                                                                      cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGAtata                                                                             <  2:1199002/151‑1 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAgga                                                                          <  1:902877/150‑3 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGt                                                                           <  2:648150/149‑1 (MQ=255)
                                                                              gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGt                                                                     <  1:680942/151‑1 (MQ=255)
                                                                              gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGt                                                                     <  1:532148/151‑1 (MQ=255)
                                                                               atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTa                                                                    <  2:593724/151‑1 (MQ=255)
                                                                                 atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCAACAAAAAACGCTGGGTGGGGGGGCGGAAAAAGCGGGAGGCGCGCCCTAATTCGGCCttttccc                                                                                    >  1:1185783/1‑129 (MQ=255)
                                                                                   ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                                <  2:480291/151‑1 (MQ=255)
                                                                                     tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTGCTGGGTGGGCTGGCGCAGAAAGGGGGAGCCGCGCCGTAAATTGGCCTTAGCCAGATATAGTATGCGTAGGcc                                                                >  1:60856/1‑149 (MQ=255)
                                                                                       gtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATTCGTAGGCCgg                                                              >  2:1180441/1‑149 (MQ=255)
                                                                                          aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAg                                                        <  2:824719/151‑1 (MQ=255)
                                                                                           aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGCGGGGTCCCTCAGAAAATGCTGGGTGGGCTGGCGCAGAAGGGGGAAGCCGCGCAGTAAATCGGCCTTAGCCAGATATATTATGCGTCGGCCGGATAAgg                                                       >  2:1142335/1‑151 (MQ=255)
                                                                                               cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGtt                                                   <  2:880706/151‑1 (MQ=255)
                                                                                               cGCTCAGCCTGGCGACCATATTCTGGCGATGAGCAACGGCGGTTCTGGTGGGGGCCCTCACAAAATGCTGGGTGGGCTGGCGAAGAGGGGGGAAGGCGCGCCGCAATTCGGCCTTAGCCAGATATAATATGGGTAGGGCGGAGAAggaggt                                                   >  1:821723/1‑147 (MQ=255)
                                                                                                gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTc                                                  <  2:1044429/151‑1 (MQ=255)
                                                                                                   cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAAAAAATGCTGGATGGGGTGGCGGAGAAGGGGGGAGGCGCGCAGGAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCAc                                                >  2:1091108/1‑150 (MQ=255)
                                                                                                   cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCCTCAGAAACTGCTGGATGGGGTGGGGAAGAAGGCGGAAGCCGCGCCGGAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAc                                                >  1:874690/1‑150 (MQ=255)
                                                                                                        tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCACAGAAACTGCTGGGAGGGGCGGGGGAGAAGGCGGGAGGCGCGCAGTAATTTGGCCCTAGCCCGATATAGTATGCGCAGGCCGGGTAAAGGGGTTCcgccgc                                           >  1:432652/1‑150 (MQ=255)
                                                                                                         ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAAATGCTGGATGGGGTGGCGGAGAAGGGGGGAGACGCGCAGTAATTTGGCCTTAGCCCGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgccg                                            >  1:158265/1‑148 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGGTGGGGTGGGGGAGAAAGCGGGAACCCCGCCGGAATTCCGCCTTAGCCCGATATAGTATGCGCAGGCCGGGGAAGGGGTTCTCGCCCCCTcc                                       >  1:1179955/1‑151 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                       <  2:989146/151‑1 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                       <  1:1125749/151‑1 (MQ=255)
                                                                                                            gACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTTGGCCCTAGCCAGATATAGTATGCGTAGGCCGGATAAAGGGTTCACGCCGCATcc                                       >  2:392205/1‑150 (MQ=255)
                                                                                                              ccATATTCTGGTGATGAGCAACGGCGGATTTGGGGGGATCCCTCAGAAAATGCTGGAGGGGGTGGCGGAGAAAGGGGGAGACGCGGAGTAAATTGGTCTTAGACCAGCATAATAAAGGGAGAGCGGCAAAAGCGTTCACGCCGCATCCg                                      >  1:631694/1‑149 (MQ=255)
                                                                                                               cATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATCCATCAGAAAAAGCTGGGTGGGGTGGCGAAAAAAGGGGGAGCCGCGCAGTAAATCGGCCTTAGCCCAATATAGTATGGGTAGCCCGGGTAAGGCGGTTACCCCGCCTcccg                                     >  2:700185/1‑147 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAGCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                   <  2:1095989/150‑1 (MQ=255)
                                                                                                                   ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGCGGGGAAAAAAGCGGAAGCCGCGCAGTAATTCGGCCTTTGCCAGATAAAGAATGCGTAGGCCGGGTAAGGGGCTCACGCCGCATCCCGCAttt                                >  2:1134038/1‑150 (MQ=255)
                                                                                                                   ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCCTCCGGCAttt                                >  2:669051/1‑150 (MQ=255)
                                                                                                                       ggTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCa                           <  2:328074/151‑1 (MQ=255)
                                                                                                                         tgatgaGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGGGGGCTGGGGGAGAAGGCGGGAGCCGCGCAGGAAATCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCCCATCCGGCAttt                                >  1:294851/1‑144 (MQ=255)
                                                                                                                                cAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTg                  <  2:845001/151‑1 (MQ=255)
                                                                                                                                  aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt                <  2:12689/151‑1 (MQ=255)
                                                                                                                                  aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCAGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt                <  1:128706/151‑1 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg             >  2:248447/1‑151 (MQ=255)
                                                                                                                                       gtttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa            <  2:602683/148‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc           <  2:644744/151‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc           <  2:316974/151‑1 (MQ=255)
                                                                                                                                          tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc         <  2:1206730/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg      >  1:1137292/1‑150 (MQ=255)
                                                                                                                                              gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGa     >  1:704690/1‑150 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc  >  1:268519/1‑151 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc  >  2:455586/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGC  >  NC_000913/4456958‑4457252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: