Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,219,640 A→C 10.3% intergenic (‑159/+52) osmF ← / ← bglX putative ABC transporter periplasmic binding protein/beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,219,6400AC10.3% 494.4 / 42.1 156intergenic (‑159/+52)osmF/bglXputative ABC transporter periplasmic binding protein/beta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (77/63);  new base C (8/8);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 160 total aligned reads displayed.

CTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGTTGAT  >  NC_000913/2219491‑2219789
                                                                                                                                                     |                                                                                                                                                     
cTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGa                                                                                                                                                       >  1:551230/1‑150 (MQ=255)
cTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGa                                                                                                                                                       <  2:1525527/150‑1 (MQ=255)
   aTTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAgcg                                                                                                                                                    >  2:2475042/1‑150 (MQ=255)
       ttttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCtt                                                                                                                                                >  1:2194259/1‑150 (MQ=255)
       ttttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATAAGGCGCTCGCGccgc                                                                                                                                                <  2:1574092/150‑3 (MQ=255)
       ttttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATAAGGCGCTCGCGccgc                                                                                                                                                <  2:599374/150‑3 (MQ=255)
       ttttACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATAAGGCGCTCGCGccgc                                                                                                                                                <  2:599361/150‑3 (MQ=255)
           aCTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                            <  1:2197773/150‑1 (MQ=255)
            cTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCg                                                                                                                                           <  2:150650/150‑1 (MQ=255)
               ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCTGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                        >  1:423333/1‑150 (MQ=255)
               ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGTATAAGGCGCTCGCGCCGCATCCGGca                                                                                                                                        <  1:2324555/150‑2 (MQ=255)
               ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                        >  1:81069/1‑150 (MQ=255)
               ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                        >  1:511615/1‑150 (MQ=255)
               ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                        >  1:321596/1‑150 (MQ=255)
               ttgttTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGcc                                                                                                                                        >  1:321587/1‑150 (MQ=255)
                  ttttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                                                                                                                                     >  1:2191107/1‑150 (MQ=255)
                  ttttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTAc                                                                                                                                     >  2:562716/1‑150 (MQ=255)
                   tttCAGCGAATTTAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCGTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACa                                                                                                                                    >  2:1104581/1‑150 (MQ=255)
                   tttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACa                                                                                                                                    >  2:38758/1‑150 (MQ=255)
                   tttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACa                                                                                                                                    >  2:1242741/1‑150 (MQ=255)
                   tttCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACa                                                                                                                                    >  1:1426435/1‑150 (MQ=255)
                     tCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACaaa                                                                                                                                  <  2:1342501/150‑1 (MQ=255)
                      cAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAg                                                                                                                                 <  2:795159/150‑1 (MQ=255)
                      cAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAg                                                                                                                                 <  2:1886622/150‑1 (MQ=255)
                        gCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAggg                                                                                                                               >  2:2221031/1‑150 (MQ=255)
                           aaTTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGgcgc                                                                                                                            >  2:638506/1‑150 (MQ=255)
                           aaTTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGgcgc                                                                                                                            <  1:554073/150‑1 (MQ=255)
                           aaTTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGgcgc                                                                                                                            >  2:1908114/1‑150 (MQ=255)
                                  agaAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTc                                                                                                                     >  1:736924/1‑150 (MQ=255)
                                  agaAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGAAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTc                                                                                                                     <  1:240055/150‑1 (MQ=255)
                                    aaaaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTcgc                                                                                                                   >  1:600056/1‑150 (MQ=255)
                                      aaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCgcgc                                                                                                                 >  2:319818/1‑150 (MQ=255)
                                      aaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCgcgc                                                                                                                 <  2:1752862/150‑1 (MQ=255)
                                      aaGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCgcgc                                                                                                                 <  2:1952426/150‑1 (MQ=255)
                                          aTAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCtt                                                                                                             >  1:168199/1‑150 (MQ=255)
                                              ttGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTa                                                                                                         >  2:1669124/1‑150 (MQ=255)
                                              ttGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTa                                                                                                         >  1:806106/1‑150 (MQ=255)
                                                gATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTaaa                                                                                                       <  2:2093834/150‑1 (MQ=255)
                                                gATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTaaa                                                                                                       >  2:1390895/1‑150 (MQ=255)
                                                   aaTGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCa                                                                                                    <  1:150311/150‑1 (MQ=255)
                                                       gCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTAc                                                                                                >  2:916913/1‑150 (MQ=255)
                                                        cGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACa                                                                                               <  1:1026105/150‑1 (MQ=255)
                                                        cGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACa                                                                                               <  1:805546/150‑1 (MQ=255)
                                                        cGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACa                                                                                               >  1:1682215/1‑150 (MQ=255)
                                                           ggTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCa                                                                                            <  2:2191111/150‑1 (MQ=255)
                                                            gTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCaa                                                                                           <  2:1286839/150‑1 (MQ=255)
                                                                  ttGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCaaa                                                                                     <  2:256856/150‑1 (MQ=255)
                                                                    gAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACt                                                                                   <  1:210716/150‑1 (MQ=255)
                                                                      gCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCg                                                                                 >  1:477010/1‑150 (MQ=255)
                                                                       cGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACACAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                                <  2:1496146/150‑1 (MQ=255)
                                                                       cGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                                <  2:2155009/150‑1 (MQ=255)
                                                                       cGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                                <  2:1496254/150‑1 (MQ=255)
                                                                       cGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGACGCATATTCGCTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                                <  2:1099899/150‑1 (MQ=255)
                                                                       cGGAAGGAGTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGc                                                                                <  2:2155286/150‑1 (MQ=255)
                                                                        ggAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAGGGGCGCAAACGTCGCGCCCTTACTAAAGCATTCCAGCACCTCAAACTCGcc                                                                               >  2:1896159/1‑150 (MQ=255)
                                                                        ggAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGACCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGcc                                                                               >  1:1164797/1‑150 (MQ=255)
                                                                          aaGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCtt                                                                             >  1:1076074/1‑150 (MQ=255)
                                                                          aaGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCtt                                                                             >  1:1521155/1‑150 (MQ=255)
                                                                          aaGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCGTACTAAAGCATTACAGCCACTCAAACTCGCCtt                                                                             >  1:260781/1‑150 (MQ=255)
                                                                                    ttGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAAcgcg                                                                   >  2:341095/1‑150 (MQ=255)
                                                                                    ttGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAAcgcg                                                                   >  2:2504889/1‑150 (MQ=255)
                                                                                    ttGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAAcgcg                                                                   >  2:650868/1‑150 (MQ=255)
                                                                                       cAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGc                                                                >  1:158406/1‑150 (MQ=255)
                                                                                            tCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAAt                                                           >  1:1928110/1‑150 (MQ=255)
                                                                                            tCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAAt                                                           >  2:1326288/1‑150 (MQ=255)
                                                                                               gATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAg                                                        <  1:916913/150‑1 (MQ=255)
                                                                                                aTAAATACTTAAGAGGAATATTCGGTCCGGCGTTTCATGCCGGATGCGGCGCGAGCGCATTATCCGGCCTACAAAGGGCGCAAACGTCGCGCACTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGt                                                       <  2:1995560/150‑1 (MQ=255)
                                                                                                  aaaTACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGc                                                     <  1:1908114/150‑1 (MQ=255)
                                                                                                           agagGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGa                                            >  1:1758410/1‑150 (MQ=255)
                                                                                                             agGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTAAAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACa                                          >  1:2117226/1‑150 (MQ=255)
                                                                                                                  tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCTTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGaa                                     >  2:2263326/1‑150 (MQ=255)
                                                                                                                  tatTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGaa                                     >  2:757791/1‑150 (MQ=255)
                                                                                                                   atTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAAc                                    <  1:1318937/150‑1 (MQ=255)
                                                                                                                   atTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAAc                                    <  2:171068/150‑1 (MQ=255)
                                                                                                                    ttCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACt                                   >  2:1261209/1‑150 (MQ=255)
                                                                                                                         tCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCa                              <  2:531835/150‑1 (MQ=255)
                                                                                                                         tCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCa                              <  2:581578/150‑1 (MQ=255)
                                                                                                                               gTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCg                        <  2:792760/150‑1 (MQ=255)
                                                                                                                               gTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCg                        <  2:879937/150‑1 (MQ=255)
                                                                                                                                  tCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCg                     <  1:1808478/150‑1 (MQ=255)
                                                                                                                                   ccTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCATTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGt                    >  1:858990/1‑150 (MQ=255)
                                                                                                                                    cTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTc                   >  2:2353190/1‑150 (MQ=255)
                                                                                                                                    cTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTc                   >  2:1823310/1‑150 (MQ=255)
                                                                                                                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAACGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCa                  <  2:103426/150‑1 (MQ=255)
                                                                                                                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAACCTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCa                  <  2:212403/150‑1 (MQ=255)
                                                                                                                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCa                  <  2:871869/150‑1 (MQ=255)
                                                                                                                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCa                  <  2:2370349/150‑1 (MQ=255)
                                                                                                                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCa                  <  2:103372/150‑1 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGCTCGCGCCGCATCCGGCCTACAAAGGGCGCAAACGACGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACAAAGAACTTGCCAGGCTCGGCGTcct                 >  1:845611/1‑148 (MQ=255)
                                                                                                                                           ataaggcgCTCGCGCCGCATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTc            >  1:2320300/5‑150 (MQ=255)
                                                                                                                                           aTGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTc            >  1:771601/1‑150 (MQ=255)
                                                                                                                                            taaggcgCTCGCGCCGCATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCa           >  1:1191585/4‑150 (MQ=255)
                                                                                                                                            tGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCa           >  1:2446046/1‑150 (MQ=255)
                                                                                                                                             gcggcgCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCAt          >  1:1740815/1‑150 (MQ=255)
                                                                                                                                             gcggcgCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCAt          <  2:1682215/150‑1 (MQ=255)
                                                                                                                                               ggcgCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCt        >  2:625986/1‑150 (MQ=255)
                                                                                                                                               ggcgCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCt        >  1:1267118/1‑150 (MQ=255)
                                                                                                                                                gcgcGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTg       >  1:450661/1‑150 (MQ=255)
                                                                                                                                                 cgcgAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGt      <  1:543377/150‑1 (MQ=255)
                                                                                                                                                    gAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGTTGa   <  2:551223/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCATTCTTTTACTGTTGTTTTCAGCGAATTAAGAGAAAAGCATAGTTGATAATGGCGGGGTTAGCTTGAGCGGAAGGATTCAGTTGCAGAATCAGATAAATACTTAAGAGGCATATTCGGTCCGGCGTTTCCTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAAAGGGCGCAAACGTCGCGCCCTTACTAAAGCATTACAGCAACTCAAACTCGCCTTTCTTAACGCGTGCGGAATCAGTGCCGATAAAGACATTGAACTTGCCAGGCTCGGCGTCATATTTCATCTGTTGAT  >  NC_000913/2219491‑2219789

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: