Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA39.5% 92.0 / 50.0 77intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (41/5);  new base A (0/30);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 2.51e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5009‑5307
                                                                                                                                                          |                                                                                                                                                      
aaTCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag                                                                                                                                                        <  2:416890/151‑1 (MQ=255)
    atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                     >  1:950295/1‑150 (MQ=255)
     tcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                                    >  2:600316/1‑150 (MQ=255)
       aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTACCTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                                   <  2:709201/149‑1 (MQ=255)
        gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaaca                                                                                                                                                >  1:1103269/1‑150 (MQ=255)
           aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                           <  2:435664/151‑1 (MQ=255)
           aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                            >  2:1219026/1‑150 (MQ=255)
            aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                          >  1:143988/1‑151 (MQ=255)
            aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                          >  1:936209/1‑151 (MQ=255)
                 aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                         <  1:530734/147‑1 (MQ=255)
                    cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                   >  2:750092/1‑150 (MQ=255)
                    cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                   >  2:562460/1‑150 (MQ=255)
                       tATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc                                                                                                                                >  1:135554/1‑150 (MQ=255)
                           tcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                            >  2:1095974/1‑150 (MQ=255)
                              aTCAGGCCGGGTTTGTTTTTATGCACCCGGTTTTTTTTAAAAAAAAATTATGGAAAAAAATAACAGGAAAAAAGGAAAAATTCTCAAAAAATGCGGTAACTTAGATATTAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                         <  2:781090/150‑1 (MQ=255)
                               tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGAGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                       >  1:923392/1‑151 (MQ=255)
                                cAGGCCGTTTTTTTTTTTACCCGGCCCGTTTTTTTAAGAAAAAAATTGGGGAAAAAAATGACAGGAAAAAAGAAAAATTTCAAAATAAAGGCTTTAACTTAGAAATTAGGATGGGAGAAAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                      <  1:473426/151‑1 (MQ=255)
                                cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                      >  2:519731/1‑151 (MQ=255)
                                 aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTa                                                                                                                           >  1:592917/1‑145 (MQ=255)
                                  ggCGGGTTTTGTTTTTAGCCACCCCGGTTTTTTTAGAAAAAAATTAGGGAAAAAAATGAAGGAAAAAAGGGAAAAATCCAAAATAAATGGGTTAATTTAGAGATAGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTa                                                                                                                           <  2:592934/144‑1 (MQ=255)
                                  ggCCGGTTTTGCTTTTATGCAGCCCGTTTTTTTTAAAAAAAAATTAGGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTCAAAAAATGCGGTAACTAAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                       <  2:407603/148‑1 (MQ=255)
                                     gggttttGTTTTTATGCAGCCCGGTTTTTTTATGAAAAAATTAGGAAAAAAAATGACGGGAAAAAAGAAAAAATTCTAAAAAAGGGGGTTATCTTAAGAATGGGGATGGGGAAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  1:426974/145‑1 (MQ=255)
                                     gggttttGCTTTTAGGCACCCCGTTTTTTTTATAAAAAATTTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAATAAATGCGGTAATTTAGAAATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 <  1:1087142/147‑1 (MQ=255)
                                     cGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 >  2:82040/1‑151 (MQ=255)
                                       ggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               >  2:195721/1‑151 (MQ=255)
                                         tttttttttATGCACCCCGTTTTTTTTAGAAAAAAATTAGGAAAAAAAATGACAGAAAAAAAGGAAAAATTCAAAAAAAATGGGATAATTTAGAAATTGGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               <  2:526032/144‑1 (MQ=255)
                                         tttGCTTTTATGCAGCCCGGCTTTTTTATGAAGCAATTATGGAGAAAACTGACAGGGAAAACGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACTGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             >  1:159802/1‑151 (MQ=255)
                                          ttGTTTTTTTGCAGCCGGATTTTTTTAGGAAAAAATATTGGAAAAAAAGGATGGGAAAAAAGGAAAATTTTCAAATAAAGCCGTAAACTAAGAGGTTAGGAGGGGGGAAAATAACCGCCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             <  2:344631/150‑1 (MQ=255)
                                          ttGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             >  2:84216/1‑150 (MQ=255)
                                            gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                          >  2:236152/1‑151 (MQ=255)
                                               tttATGCAGCCCGTTTTTTTTAAAAAAAAATTATGGAAAAAAATGAAAGGAAAAAAGAAAAAATCCAAAAAAAATGCGTTACTTAAGATATAAGGTTTGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                           <  1:300369/147‑1 (MQ=255)
                                               tttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                        >  1:1170052/1‑150 (MQ=255)
                                                 tATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      >  2:939796/1‑150 (MQ=255)
                                                  aggccgCCGGGTTTTTTAAAAAAAAAATGATGGAAAAAAAGAAAAGGAAAAAAGAAAAATTTCCAAATAAATGCGGTAATTTAGAAATTAGGTTTGCGGAAAAAACCAACCGCCGTTCTCATCTAGTAATCTCCGGATATCGACCCATAAc                                                                                                    <  1:74861/146‑1 (MQ=255)
                                                   gccgccccGGTTTTTTTATGAAAAAATTATGGAAAAAAATAACGGGAAAAAAGGAGAATTCCTAAAAAACTGGGGAAATTAAGGGTGGAGATGGGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                   <  2:592356/146‑1 (MQ=255)
                                                        cccgttttttttATAAAAAATTTAGGGAAAAAAATGACAGGAAAAAAGGAGAAATTCTAAATAAATGGGTAAATTTAAAAATAAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                <  2:1120901/143‑1 (MQ=255)
                                                        cccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                               >  2:74219/1‑150 (MQ=255)
                                                         ccctgttttttAAAAAAAAAATAAGGAAAAAAAAGGAAAGAAAAAAAGAAAAAATTCTAAAAAAATGCGATAATTTAGAAATTAGGATTGGGAAAAAAAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                               <  1:658213/145‑1 (MQ=255)
                                                           ggtttttttAAAAAAAAATTAGGAAAAAAAAGGACAGGAAAAAGGAAAAAATTCAAAAAAAATGCGGTAATTTAGAGATTGGGTTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                           <  1:1199939/148‑1 (MQ=255)
                                                                           aattaatgaaaaaaaaGGACAGAAAAAAAGGAAAAATTCTAAAAAAAGGGGTTAATTAAGAGATGAGGATGGGGGAAAAAACAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                            <  1:398745/140‑1 (MQ=255)
                                                                           aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  1:176675/1‑151 (MQ=255)
                                                                           aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  2:68572/1‑151 (MQ=255)
                                                                           aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  2:494773/1‑151 (MQ=255)
                                                                             ttATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         >  1:1023394/1‑151 (MQ=255)
                                                                                  gagaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                     >  1:850862/1‑150 (MQ=255)
                                                                                  gaaaaaaaaTGAAAGGAAAAAGGAAAAAATTCAAAAAAAAGGCGGAATTTAAGAAATTAGATTGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         <  1:452393/143‑1 (MQ=255)
                                                                                  gaaaaaaaaGGAAAGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGTAACTTAGATTTTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        <  1:303507/144‑1 (MQ=255)
                                                                                        aaaagagagaaaaaaaGGAAAAATTCTAAAAAAAGGGGGTAATTTAGAGATTAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               <  2:1061715/140‑1 (MQ=255)
                                                                                        aaaTGAGGGGAAAAAGGGAAAATTTCAAAAAAAATGCGTAATCTAAGAAATTAGTATGGGGAAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               <  2:234417/151‑1 (MQ=255)
                                                                                         aaggaaaggaaaaaaggaaaaaTTCAAAAAAAATGCGTTAACTAAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                              <  2:686359/145‑1 (MQ=255)
                                                                                         aaaaaCAGGAAAAAAGGAAAAATTCAAAAAAATGGCGGAAACTGAGAGTTTAGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATgaa                                                              <  2:1103335/147‑3 (MQ=255)
                                                                                          ataaaaGGAAAAAAGAAAAAATTCTAAAAAAAGGCGGTAATTTAGAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                             <  2:1041196/146‑1 (MQ=255)
                                                                                                ggaaaaaaGGAAAAATTCAAAAAAAAGGCGATAACTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                         <  2:809433/146‑1 (MQ=255)
                                                                                                  gAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     >  1:664692/1‑151 (MQ=255)
                                                                                                  aaaaaaGGAAAAATTCTAAAAAAAGGCGGTAACTTAGAGATAAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     <  2:143994/150‑1 (MQ=255)
                                                                                                      aaGGAAAAATTCTAAAAAAATGCGGTAATTAAGAGATTGGGATTGGGGAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                 <  2:472912/151‑1 (MQ=255)
                                                                                                           aaaaaTTCTAAAAAAAGGGGTTAACTTAGAGATTAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                               <  1:609469/147‑1 (MQ=255)
                                                                                                            gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                            >  1:818467/1‑151 (MQ=255)
                                                                                                                    taaaaaaaGGCGGAAATTAAGAAATTGGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                        <  1:689757/143‑1 (MQ=255)
                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                  >  2:405747/1‑151 (MQ=255)
                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                  >  2:608897/1‑151 (MQ=255)
                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   >  1:1171886/1‑150 (MQ=255)
                                                                                                                         aaaTGCGGAAAGTAAGAAATAGGAATGGGGAAAAAAACAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                <  1:936058/151‑1 (MQ=255)
                                                                                                                          aaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               >  1:640778/1‑151 (MQ=255)
                                                                                                                             gggaaaattaaaggaTGAGTATGGCGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               <  1:737209/135‑1 (MQ=255)
                                                                                                                             gcttaaagtaagagATTGGAATGGGGAAAAAAAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            <  2:1208732/141‑1 (MQ=255)
                                                                                                                             gCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            >  2:1218336/1‑151 (MQ=255)
                                                                                                                               ggTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGc                          >  2:20646/1‑151 (MQ=255)
                                                                                                                                  aaCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                       >  2:585722/1‑151 (MQ=255)
                                                                                                                                   aCTTAGTGATTAGGATTGCGGAGAATAACAACTGCCGTTCTCATCGAGTAATTTCCGGATATCGACCCATAACGGGCAATGATAAAACGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCGTGGTTCTGGTCGCTc                        >  1:14245/1‑149 (MQ=255)
                                                                                                                                     ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    >  1:343785/1‑151 (MQ=255)
                                                                                                                                     taagatATTGGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                           <  2:1063203/142‑1 (MQ=255)
                                                                                                                                             tAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac            >  2:461839/1‑151 (MQ=255)
                                                                                                                                                   tGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg      >  2:633790/1‑151 (MQ=255)
                                                                                                                                                    ggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTGTCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc     <  2:762209/149‑1 (MQ=255)
                                                                                                                                                      ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg   >  1:31546/1‑151 (MQ=255)
                                                                                                                                                      ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTATGGTCGCTCCCATGGCAGCACAGGCTGCgg   >  2:690120/1‑151 (MQ=255)
                                                                                                                                                      ggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg   <  2:224328/151‑1 (MQ=255)
                                                                                                                                                         agaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  1:1127758/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5009‑5307

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