Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,566 A→T 6.7% F183L (TTT→TTA dhaM ← putative dihydroxyacetone‑specific PTS enzymes: HPr, EI components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,5660AT 6.7% 277.0 / 6.5 98F183L (TTT→TTAdhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base T (2/5);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 104 total aligned reads displayed.

TGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGT  >  NC_000913/1248417‑1248716
                                                                                                                                                     |                                                                                                                                                      
tGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTa                                                                                                                                                        >  2:346129/1‑150 (MQ=255)
tGACGGAAAGCGAGCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTa                                                                                                                                                        <  1:1266223/150‑1 (MQ=255)
 gACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaa                                                                                                                                                       >  2:120157/1‑150 (MQ=255)
  aCGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTaaa                                                                                                                                                      <  1:155294/150‑1 (MQ=255)
   cGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATTCTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATAGGCATTAAAt                                                                                                                                                     >  2:286689/1‑150 (MQ=255)
   cGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACACATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAAt                                                                                                                                                     >  1:241846/1‑150 (MQ=255)
     gAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATgc                                                                                                                                                   <  1:1726298/150‑2 (MQ=255)
      aaaGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTc                                                                                                                                                  <  1:437974/150‑1 (MQ=255)
         gCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGAt                                                                                                                                               >  2:732464/1‑150 (MQ=255)
          cGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATa                                                                                                                                              <  2:588145/150‑1 (MQ=255)
           gATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATaa                                                                                                                                             >  1:1694804/1‑150 (MQ=255)
           gATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATat                                                                                                                                             <  1:760024/150‑2 (MQ=255)
            aTCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                                                                           <  1:13301/150‑1 (MQ=255)
             tCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATAATTCGCATCGTTTCGTTATAGCGAACTTGTAGTATAGCAATCTGGTTAATACTCTCTGGTGTGACTCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAgg                                                                                                                                           <  1:195983/150‑1 (MQ=255)
             tCAGTGCCTCATCAGAATCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGAGTAATTGTAGTAACTCAATCTGTTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTACCAGCAACATATCGACATTTAATGCCGATatgt                                                                                                                                           <  2:227050/150‑4 (MQ=255)
               aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtt                                                                                                                                         >  2:868878/1‑149 (MQ=255)
               aGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttg                                                                                                                                         >  2:165193/1‑145 (MQ=255)
                gTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtat                                                                                                                                        >  2:1126718/1‑150 (MQ=255)
                 tGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACTCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGtata                                                                                                                                       <  1:104219/150‑1 (MQ=255)
                 tGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTTAATGCCGATatgttgct                                                                                                                                       <  2:936297/150‑8 (MQ=255)
                    ctctTCAGCTTCTGGCGCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCACTGGTGTGACGCATTTGCCGTTTTTTTCCAGCATCATATCGGCATTAAATGTCGATAAGGTATAAAc                                                                                                                                    <  2:337293/150‑1 (MQ=255)
                     tctTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCTAACTTGAAGTAACGCAATCTGGTTAATACTCTCTGGTGTTACGCATTTGCCGTTTTTTTCCAGCAACATCTCGGCATTAAATGTCGATAAGGTATAAAcc                                                                                                                                   <  1:337848/150‑1 (MQ=255)
                      ctTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTACAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCa                                                                                                                                  >  2:1401461/1‑150 (MQ=255)
                       ttCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCac                                                                                                                                 >  2:1275236/1‑149 (MQ=255)
                         cAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTAAATGTCGATAAGGTATAAACCAGcc                                                                                                                               >  2:821471/1‑150 (MQ=255)
                           gCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCgg                                                                                                                             <  1:623804/150‑1 (MQ=255)
                                   cccTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCgg                                                                                                                     >  1:1120346/1‑150 (MQ=255)
                                    ccTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGa                                                                                                                    <  1:732464/150‑1 (MQ=255)
                                     cTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                                                                                   <  2:643552/150‑1 (MQ=255)
                                     cTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGAc                                                                                                                   <  2:1184119/150‑1 (MQ=255)
                                      tttCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACg                                                                                                                  <  1:55425/150‑1 (MQ=255)
                                       ttCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGt                                                                                                                 >  2:1375175/1‑150 (MQ=255)
                                       ttCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGt                                                                                                                 >  2:1368447/1‑150 (MQ=255)
                                         cGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGTCATTAAATGTCGATAAGGTTTAAAACATCCGGGAGGCCGGACGTaa                                                                                                               >  1:1216326/1‑149 (MQ=255)
                                         cGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGTCCGGACGTAc                                                                                                               >  1:213357/1‑150 (MQ=255)
                                         cGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAAAACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTAAAACCCAGCCGGGAGGCCGGACGTAc                                                                                                               >  1:156536/1‑150 (MQ=255)
                                             aTCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGc                                                                                                           >  1:428138/1‑150 (MQ=255)
                                                 ggCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGc                                                                                                       <  1:1654223/150‑1 (MQ=255)
                                                     cAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGtt                                                                                                   >  2:210861/1‑150 (MQ=255)
                                                     cAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGtt                                                                                                   <  1:1265589/150‑1 (MQ=255)
                                                             cGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGtttt                                                                                           <  1:346129/150‑1 (MQ=255)
                                                              gTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTGACGGttttt                                                                                          >  2:1222824/1‑150 (MQ=255)
                                                               ttATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGtttttt                                                                                         <  2:1200753/150‑1 (MQ=255)
                                                                tataGCGAACTTGTAGTAACGCAATCTGGGTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTACATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTa                                                                                        <  1:1087139/150‑1 (MQ=255)
                                                                    gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATTTCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                                    >  2:1516895/1‑150 (MQ=255)
                                                                    gCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGac                                                                                    >  2:1466378/1‑150 (MQ=255)
                                                                     cGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGacc                                                                                   <  1:1691306/150‑1 (MQ=255)
                                                                     cGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGacc                                                                                   <  2:1120346/150‑1 (MQ=255)
                                                                           tGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCAGTTACGGTTTTTTATGACCACCGcc                                                                             >  1:198882/1‑150 (MQ=255)
                                                                            gtagtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCa                                                                            >  2:1335571/1‑150 (MQ=255)
                                                                             tagtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCag                                                                           <  2:1168427/150‑1 (MQ=255)
                                                                              agtaACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATGTCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCaga                                                                          >  2:593274/1‑150 (MQ=255)
                                                                               gtaACGCAATCTGGTTAATACTCTCTGGTGTGACTCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCTGGAGGCCGGCCGTACATGCAGGCCGTTACTGGTGTTTATGCCCACCGCCtgag                                                                         >  2:90172/1‑150 (MQ=255)
                                                                                       aTCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGct                                                                 <  1:1495569/150‑1 (MQ=255)
                                                                                        tCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGctt                                                                >  1:515939/1‑150 (MQ=255)
                                                                                         cTGGTTAATACTCTCTGGTGTGACGCATTTGCAGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCTTTACGGTTTTTTATGACCACCGCCAGACAACGGGctta                                                               >  1:1591330/1‑149 (MQ=255)
                                                                                          tGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGcttct                                                              <  1:405648/150‑1 (MQ=255)
                                                                                            gTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCttcttc                                                            >  2:1037716/1‑150 (MQ=255)
                                                                                               aaTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATc                                                         >  2:75493/1‑150 (MQ=255)
                                                                                                aTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCCGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCCGCCGGGAGGCCGGACGTACATGCAGGCCGATACGGTTTTTTAAGACCACCGCCAGAGAACGGGCTTCTTCATCg                                                        >  1:875278/1‑150 (MQ=255)
                                                                                                aTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCATGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg                                                        <  1:1470193/150‑1 (MQ=255)
                                                                                                aTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg                                                        <  1:1026255/150‑1 (MQ=255)
                                                                                                aTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg                                                        <  1:538650/150‑1 (MQ=255)
                                                                                                aTACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCg                                                        <  1:545543/150‑1 (MQ=255)
                                                                                                   ctctctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTAc                                                     <  2:6676/150‑1 (MQ=255)
                                                                                                     ctctGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCATCCGGGAGTCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTACATCGTACGc                                                   >  1:658533/1‑150 (MQ=255)
                                                                                                         ggTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGa                                               >  2:243480/1‑150 (MQ=255)
                                                                                                         ggTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGa                                               <  1:1332222/150‑1 (MQ=255)
                                                                                                          gggtgACCCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGAc                                              <  1:286689/148‑1 (MQ=255)
                                                                                                           tgtgACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACa                                             <  2:22846/150‑1 (MQ=255)
                                                                                                             tgACGCATTTGCCGTTTTTTTCCAGCAACATATCGACATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtg                                           <  1:821471/150‑1 (MQ=255)
                                                                                                               aCGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAtgtg                                         <  1:1037716/150‑1 (MQ=255)
                                                                                                                     ttGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCagag                                   >  1:576736/1‑150 (MQ=255)
                                                                                                                       gCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGAt                                 >  1:1290413/1‑150 (MQ=255)
                                                                                                                          gTTTTTTTCCAGCAACATGTCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTc                              <  1:593274/150‑1 (MQ=255)
                                                                                                                            ttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                            <  1:939546/150‑1 (MQ=255)
                                                                                                                            ttttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGCCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                            <  2:253926/150‑1 (MQ=255)
                                                                                                                             tttttCCAGCAACTTATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAgt                           <  1:1368447/150‑1 (MQ=255)
                                                                                                                             tttttCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAgt                           <  1:1375175/150‑1 (MQ=255)
                                                                                                                                  ccAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCATGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCgg                      <  1:1126718/150‑1 (MQ=255)
                                                                                                                                  ccAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCgg                      >  1:785348/1‑150 (MQ=255)
                                                                                                                                   cAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCgga                     >  2:121285/1‑150 (MQ=255)
                                                                                                                                        ccttatCGACATTTAATGCCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACg                >  1:189394/4‑150 (MQ=255)
                                                                                                                                          atatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTACAGTGTCGGAGGACGGt              >  2:777830/1‑150 (MQ=255)
                                                                                                                                           tatCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTa             <  2:1173367/150‑1 (MQ=255)
                                                                                                                                              cGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTCACGGTTTGTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTCCGCAGGACATGTGGCAGAGCTGTCAGTGTCGGAGGACGGTaacc          >  1:12887/1‑148 (MQ=255)
                                                                                                                                               ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCTTACGCCGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAAcc         >  1:152840/1‑150 (MQ=255)
                                                                                                                                               ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAAcc         >  1:526909/1‑150 (MQ=255)
                                                                                                                                               ggCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAAcc         <  1:75493/150‑1 (MQ=255)
                                                                                                                                                gCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGACCGTGCTTCTTCAACGTACGCAGGTCATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAAccc        >  2:1093201/1‑150 (MQ=255)
                                                                                                                                                 cATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                            >  2:1408489/1‑129 (MQ=255)
                                                                                                                                                 cATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAg                            <  1:1408489/129‑1 (MQ=255)
                                                                                                                                                 cATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCa       >  1:1562538/1‑150 (MQ=255)
                                                                                                                                                  aatAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTTAGTGTCGCAGGACGGTAAACCCAg      <  2:630179/148‑1 (MQ=255)
                                                                                                                                                  aTTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCGTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAg      <  1:1275236/150‑1 (MQ=255)
                                                                                                                                                  aTTAAATGGCGATAAGGTATAAACCAGCCGGGAGGCCGTACGTACATGCCGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAg      <  2:931301/150‑1 (MQ=255)
                                                                                                                                                   ttAATTGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGTCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATTTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgt     <  1:353837/150‑1 (MQ=255)
                                                                                                                                                   ttAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgt     <  2:1618565/150‑1 (MQ=255)
                                                                                                                                                     aaaTGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCCCCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACAAGTGTCAGAGCTTTCAGTGTCGGAGGACGTTAACCCAggtga  >  2:495969/1‑149 (MQ=255)
                                                                                                                                                     aaaTGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAgttg   <  1:1027804/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGACGGAAAGCGATCAGTGCCTCTTCAGCTTCTGGCCCTTTCGCAATCAGGCGCAGCGTATCGTTATAGCGAACTTGTAGTAACGCAATCTGGTTAATACTCTCTGGTGTGACGCATTTGCCGTTTTTTTCCAGCAACATATCGGCATTAAATGTCGATAAGGTATAAACCAGCCGGGAGGCCGGACGTACATGCAGGCCGTTACGGTTTTTTATGACCACCGCCAGAGAACGGGCTTCTTCATCGTACGCAGGACATGTGTCAGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGT  >  NC_000913/1248417‑1248716

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