Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 263,353 T→A 5.8% Q106L (CAA→CTA) ‡ ykfI ← CP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913263,3530TA 5.8% 497.8 / 13.9 154Q106L (CAA→CTA) ‡ykfICP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (83/62);  new base A (2/7);  total (85/69)
Fisher's exact test for biased strand distribution p-value = 7.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 154 total aligned reads displayed.

GGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTG  >  NC_000913/263205‑263499
                                                                                                                                                    |                                                                                                                                                  
ggTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTg                                                                                                                                                   >  2:1833626/1‑150 (MQ=255)
ggTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTg                                                                                                                                                   >  2:2128472/1‑150 (MQ=255)
     cATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                                                              >  1:1744263/1‑150 (MQ=255)
     cATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                                                              >  1:1744262/1‑150 (MQ=255)
       tGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaac                                                                                                                                            <  1:2367104/150‑1 (MQ=255)
        gATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaaca                                                                                                                                           >  1:2237860/1‑150 (MQ=255)
         aTTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaacaa                                                                                                                                          >  1:1192621/1‑150 (MQ=255)
          ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaac                                                                                                                                         <  2:197263/150‑2 (MQ=255)
           tgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaacg                                                                                                                                        <  1:525724/150‑3 (MQ=255)
                  ccGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgc                                                                                                                                 <  1:415872/150‑1 (MQ=255)
                   cGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcc                                                                                                                                >  1:555706/1‑150 (MQ=255)
                    gTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcct                                                                                                                               >  2:1013764/1‑150 (MQ=255)
                     tAAAAAGATAACGCGCCTGCCGGGAAGTAGCAGGCTCATTACGCATTAGGTAAACAAAGGAGCAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                                                              >  2:453106/1‑150 (MQ=255)
                      aaaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctgc                                                                                                                             >  1:1985040/1‑150 (MQ=255)
                      aaaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctgc                                                                                                                             >  2:1197253/1‑150 (MQ=255)
                       aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctgcc                                                                                                                            >  1:1248280/1‑150 (MQ=255)
                       aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCCTTGCCGCAACAAGCCAGTtgcctgcc                                                                                                                            >  2:2142786/1‑150 (MQ=255)
                          aGATAACGCGCCTGCCGTGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCcgcg                                                                                                                         <  2:1661359/150‑1 (MQ=255)
                             tAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTc                                                                                                                      <  1:1141877/150‑1 (MQ=255)
                                cgcgCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCa                                                                                                                   >  1:661361/1‑150 (MQ=255)
                                  cgcCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGa                                                                                                                 >  2:2016399/1‑150 (MQ=255)
                                    ccTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAAt                                                                                                               >  1:275033/1‑150 (MQ=255)
                                      tGCCGGGTAGTAGCAGGCACATTACGCAATAGGTAAACAAGGGGGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGt                                                                                                             >  2:710041/1‑150 (MQ=255)
                                       gCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTc                                                                                                            >  2:1781928/1‑150 (MQ=255)
                                         cGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCACCAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTg                                                                                                          >  2:1347494/1‑150 (MQ=255)
                                         cGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTg                                                                                                          <  1:1606758/150‑1 (MQ=255)
                                           ggTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCa                                                                                                        >  1:440380/1‑150 (MQ=255)
                                               gtagCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGccc                                                                                                    >  2:1512044/1‑150 (MQ=255)
                                               gtagCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGccc                                                                                                    >  1:1353823/1‑150 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACTGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  2:1209501/150‑1 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  2:1548391/1‑150 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  2:383208/150‑1 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  2:948997/1‑150 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  1:1461071/150‑1 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  1:1281410/1‑150 (MQ=255)
                                                 agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  1:1366531/1‑150 (MQ=255)
                                                  gCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGga                                                                                                 >  2:1580232/1‑150 (MQ=255)
                                                  gCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGga                                                                                                 >  2:863043/1‑150 (MQ=255)
                                                     ggCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGAt                                                                                              >  2:1167053/1‑150 (MQ=255)
                                                      gcgcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATa                                                                                             >  2:626417/1‑150 (MQ=255)
                                                       cgcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                                                            >  2:578381/1‑150 (MQ=255)
                                                       cgcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                                                            >  2:1874634/1‑150 (MQ=255)
                                                         cATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          <  2:491918/150‑1 (MQ=255)
                                                           ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGgag                                                                                        <  2:1985040/150‑1 (MQ=255)
                                                             aCGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGAt                                                                                      <  1:1197253/150‑1 (MQ=255)
                                                               gCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTg                                                                                    <  1:607898/150‑1 (MQ=255)
                                                                cAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGt                                                                                   >  2:818086/1‑150 (MQ=255)
                                                                cAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGt                                                                                   >  2:155083/1‑150 (MQ=255)
                                                                 aaTAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGtt                                                                                  >  2:1719725/1‑150 (MQ=255)
                                                                    aGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCt                                                                               >  1:2188850/1‑150 (MQ=255)
                                                                     ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  1:1373303/1‑150 (MQ=255)
                                                                     ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  1:1016706/1‑150 (MQ=255)
                                                                     ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  2:2043351/1‑150 (MQ=255)
                                                                     ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  1:519465/1‑150 (MQ=255)
                                                                      gTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGc                                                                             >  1:2362402/1‑150 (MQ=255)
                                                                        aaaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCa                                                                           <  2:1614576/150‑1 (MQ=255)
                                                                         aaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAg                                                                          >  2:1285805/1‑150 (MQ=255)
                                                                            aaGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTa                                                                       <  2:661361/150‑1 (MQ=255)
                                                                             aGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTaa                                                                      <  2:1383054/150‑1 (MQ=255)
                                                                             aGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTaa                                                                      >  1:1584955/1‑150 (MQ=255)
                                                                                  ggaAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATccc                                                                 >  2:724043/1‑150 (MQ=255)
                                                                                     aGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTc                                                              >  1:2086717/1‑150 (MQ=255)
                                                                                      gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             <  2:2237860/150‑1 (MQ=255)
                                                                                      gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             >  2:1564106/1‑150 (MQ=255)
                                                                                      gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             <  1:658343/150‑1 (MQ=255)
                                                                                       ttCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCt                                                            >  1:1173640/1‑150 (MQ=255)
                                                                                        tCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTg                                                           >  1:2157049/1‑150 (MQ=255)
                                                                                         cAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGt                                                          >  1:1022195/1‑150 (MQ=255)
                                                                                          aGACATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGGTCGCAGATTGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTc                                                         <  2:741467/150‑1 (MQ=255)
                                                                                          aGAAATGTAAATCGTGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTc                                                         <  2:749808/150‑1 (MQ=255)
                                                                                           gAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCg                                                        <  1:1518154/150‑1 (MQ=255)
                                                                                           gAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCg                                                        >  1:2327100/1‑150 (MQ=255)
                                                                                               tGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGAGCGGAGATAAGTAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                    <  2:1744262/150‑1 (MQ=255)
                                                                                                  aaaTCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGaa                                                 >  1:90719/1‑150 (MQ=255)
                                                                                                     tCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCa                                              >  2:1078729/1‑150 (MQ=255)
                                                                                                             gTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCATCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAt                                      >  2:2183555/1‑150 (MQ=255)
                                                                                                                gTACGCAATGTTCATCGTAGTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCTCAGAATTTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCATTCCTGTCGATACGAACCAGCTCGTAtttt                                   <  2:1345162/150‑1 (MQ=255)
                                                                                                                  aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          <  2:1975558/93‑1 (MQ=255)
                                                                                                                  aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          >  1:1975558/1‑93 (MQ=255)
                                                                                                                  aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTc                                 <  2:785788/150‑1 (MQ=255)
                                                                                                                     cAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTAc                              >  2:1880702/1‑150 (MQ=255)
                                                                                                                           tCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaaa                        >  2:870720/1‑150 (MQ=255)
                                                                                                                             aTCGTACTACGTTGGTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAAt                      >  2:1409098/1‑150 (MQ=255)
                                                                                                                                 tactacTTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATATGGAGATTGTTGCTGCCAGCTTAATCCCTTGCTGACGATACGAACCAGATCGTATTTTTCTACCAGAAAATTCAc                  <  2:2026306/150‑1 (MQ=255)
                                                                                                                                  actacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACg                 >  1:150941/1‑150 (MQ=255)
                                                                                                                                   ctggcttgttgCGGCAAAGCCGTAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                >  2:417379/6‑150 (MQ=255)
                                                                                                                                   ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                <  1:2060308/150‑1 (MQ=255)
                                                                                                                                   ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                >  2:2334551/1‑150 (MQ=255)
                                                                                                                                   ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                <  1:870720/150‑1 (MQ=255)
                                                                                                                                       gttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATc            >  1:2334819/1‑150 (MQ=255)
                                                                                                                                       cttgttgCGGCAAAGCCGTAACAAGCCAGTTTCCTGCCGCGGTCTCAGAATGTCTGGAGACCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGATCCAGCTCGTATTTTTATACCAGAAAATTCACGGCAta            >  2:1641917/2‑149 (MQ=255)
                                                                                                                                          gttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc         <  2:150941/150‑1 (MQ=255)
                                                                                                                                            tACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTa       >  1:911052/1‑150 (MQ=255)
                                                                                                                                              cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCATGGCATCGGCTAgg     <  1:380083/150‑1 (MQ=255)
                                                                                                                                              cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     <  1:994934/150‑1 (MQ=255)
                                                                                                                                              cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     <  1:761184/150‑1 (MQ=255)
                                                                                                                                              cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     <  1:1314957/150‑1 (MQ=255)
                                                                                                                                               ggCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCGCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAggg    <  2:1310435/150‑1 (MQ=255)
                                                                                                                                                gCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGt   >  1:2437436/1‑150 (MQ=255)
                                                                                                                                                 cTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAGAATTCACGGCATCGGCTAGGGTg  >  2:1100999/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTG  >  NC_000913/263205‑263499

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