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

ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGC  >  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)
                                                                tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCTCCGAAACTGCTGGGTGGGGCGGGGGAGAAAGCGGGAGCCGCGCCGGAAATAGGGCTTAG‑cc                                                                                    >  1:880667/1‑150 (MQ=255)
                                                                tgGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGTGGGGTCCCTCCGAAACTGCTGGGTGGGCTGGCGAAGAAAGCGGGAGCCGCGCCGTAAATTGGCCCTTG‑ccc                                                                                   >  2:145377/1‑150 (MQ=255)
                                                                 ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCATCAGAAAATGCTGGAAGGGCTCGCGAAAAAAACGGGAGCCGCGCAGTAATTCCGCCTTAG‑ACa                                                                                   >  2:1221829/1‑150 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGGGGAGAAAGGGGGAGGCGCGCAGTAAATCGGCCTTAG‑cc                                                                                    >  1:76012/1‑147 (MQ=255)
                                                                   aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGTCCCTCAAAAACAGCTGGGATGGCTGGGGGAAAAAGGGGGAGCCCCGCAGGAATTTGGccctttcc                                                                                     >  1:325181/1‑141 (MQ=255)
                                                                    tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGGGGGAGGCGCGCCGGAATTTCGCCTTTG‑CCAg                                                                                  >  1:922164/1‑148 (MQ=255)
                                                                     aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGatat                                                                              <  1:1029729/151‑1 (MQ=255)
                                                                      cGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGAtata                                                                             <  2:1199002/151‑1 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAgga                                                                          <  1:902877/150‑3 (MQ=255)
                                                                          ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGt                                                                           <  2:648150/149‑1 (MQ=255)
                                                                              gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGt                                                                     <  1:680942/151‑1 (MQ=255)
                                                                              gATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGt                                                                     <  1:532148/151‑1 (MQ=255)
                                                                               atatGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTa                                                                    <  2:593724/151‑1 (MQ=255)
                                                                                 atGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCAACAAAAAACGCTGGGTGGGGGGGCGGAAAAAGCGGGAGGCGCGCCCTAATTCGGCCttttccc                                                                                    >  1:1185783/1‑129 (MQ=255)
                                                                                   ggtggtGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGcc                                                                <  2:480291/151‑1 (MQ=255)
                                                                                     tggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGGCCATCAGAAACTGCTGGGTGGGCTGGCGCAGAAAGGGGGAGCCGCGCCGTAAATTGGCCTTAG‑CCAGATATAGTATGCGTAGGcc                                                                >  1:60856/1‑149 (MQ=255)
                                                                                       gtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATTCGTAGGCCgg                                                              >  2:1180441/1‑149 (MQ=255)
                                                                                          aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAg                                                        <  2:824719/151‑1 (MQ=255)
                                                                                           aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGCGGGGTCCCTCAGAAAATGCTGGGTGGGCTGGCGCAGAAGGGGGAAGCCGCGCAGTAAATCGGCCTTAG‑CCAGATATATTATGCGTCGGCCGGAT‑AAgg                                                       >  2:1142335/1‑151 (MQ=255)
                                                                                               cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGtt                                                   <  2:880706/151‑1 (MQ=255)
                                                                                               cGCTCAGCCTGGCGACCATATTCTGGCGATGAGCAACGGCGGTTCTGGTGGGGGCCCTCACAAAATGCTGGGTGGGCTGGCGAAGAGGGGGGAAGGCGCGCCGCAATTCGGCCTTAG‑CCAGATATAATATGGGTAGGGCGGAG‑AAggaggt                                                   >  1:821723/1‑147 (MQ=255)
                                                                                                gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTc                                                  <  2:1044429/151‑1 (MQ=255)
                                                                                                   cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAAAAAATGCTGGATGGGGTGGCGGAGAAGGGGGGAGGCGCGCAGGAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGGT‑AAGGCGTTCAc                                                >  2:1091108/1‑150 (MQ=255)
                                                                                                   cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCCTCAGAAACTGCTGGATGGGGTGGGGAAGAAGGCGGAAGCCGCGCCGGAAATCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCAc                                                >  1:874690/1‑150 (MQ=255)
                                                                                                        tGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCACAGAAACTGCTGGGAGGGGCGGGGGAGAAGGCGGGAGGCGCGCAGTAATTTGGCCCTAG‑CCCGATATAGTATGCGCAGGCCGGGTAAAGGGGTTC‑cgccgc                                           >  1:432652/1‑150 (MQ=255)
                                                                                                         ggCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAAATGCTGGATGGGGTGGCGGAGAAGGGGGGAGACGCGCAGTAATTTGGCCTTAG‑CCCGATATAGTATGCGTAGGCCGGGT‑AAGGCGGTCAcgccg                                            >  1:158265/1‑148 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGGTGGGGTGGGGGAGAAAGCGGGAACCCCGCCGGAATTCCGCCTTAG‑CCCGATATAGTATGCGCAGGCCGGGG‑AAGGGGTTCTCGCCCCCTcc                                       >  1:1179955/1‑151 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATcc                                       <  2:989146/151‑1 (MQ=255)
                                                                                                           cGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATcc                                       <  1:1125749/151‑1 (MQ=255)
                                                                                                            gACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTTGGCCCTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAAGGGTTCACGCCGCATcc                                       >  2:392205/1‑150 (MQ=255)
                                                                                                              ccATATTCTGGTGATGAGCAACGGCGGATTTGGGGGGATCCCTCAGAAAATGCTGGAGGGGGTGGCGGAGAAAGGGGGAGACGCGGAGTAAATTGGTCTTAGACCAG‑CATAATAAAGGGAGAGCGGCA‑AAAGCGTTCACGCCGCATCCg                                      >  1:631694/1‑149 (MQ=255)
                                                                                                               cATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGATCCATCAGAAAAAGCTGGGTGGGGTGGCGAAAAAAGGGGGAGCCGCGCAGTAAATCGGCCTTAG‑CCCAATATAGTATGGGTAGCCCGGGT‑AAGGCGGTTACCCCGCCTcccg                                     >  2:700185/1‑147 (MQ=255)
                                                                                                                atatTCTGGTGATGAGCAGCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCa                                   <  2:1095989/150‑1 (MQ=255)
                                                                                                                   ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGGTGGGGCGGGGAAAAAAGCGGAAGCCGCGCAGTAATTCGGCCTTTG‑CCAGATAAAGAATGCGTAGGCCGGGT‑AAGGGGCTCACGCCGCATCCCGCAttt                                >  2:1134038/1‑150 (MQ=255)
                                                                                                                   ttCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGGT‑AAGGCGTTCACGCCGCCTCCGGCAttt                                >  2:669051/1‑150 (MQ=255)
                                                                                                                       ggTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCa                           <  2:328074/151‑1 (MQ=255)
                                                                                                                         tgatgaGCAACGGCGGTTTTTGTGGGGTCCCTCAGAAACTGCTGGGGGGGCTGGGGGAGAAGGCGGGAGCCGCGCAGGAAATCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGGT‑AAGGCGTTCACGCCCCATCCGGCAttt                                >  1:294851/1‑144 (MQ=255)
                                                                                                                                cAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTg                  <  2:845001/151‑1 (MQ=255)
                                                                                                                                  aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt                <  2:12689/151‑1 (MQ=255)
                                                                                                                                  aCGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCAGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt                <  1:128706/151‑1 (MQ=255)
                                                                                                                                     gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg             >  2:248447/1‑151 (MQ=255)
                                                                                                                                       gtttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGa            <  2:602683/148‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc           <  2:644744/151‑1 (MQ=255)
                                                                                                                                       ggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc           <  2:316974/151‑1 (MQ=255)
                                                                                                                                          tttGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc         <  2:1206730/150‑1 (MQ=255)
                                                                                                                                             ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg      >  1:1137292/1‑150 (MQ=255)
                                                                                                                                              gtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGa     >  1:704690/1‑150 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc  >  1:268519/1‑151 (MQ=255)
                                                                                                                                                gggATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGc  >  2:455586/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAG‑CCAGATATAGTATGCGTAGGCCGGAT‑AAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGC  >  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: ‑