Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,886,725 G→T 100% I256I (ATC→ATA yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,886,7250GT97.7% 699.0 / 5.8 220I256I (ATC→ATAyoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (111/104);  total (114/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 220 total aligned reads displayed.

TGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT  >  CP009273/1886579‑1886870
                                                                                                                                                  |                                                                                                                                                 
tGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGAt                                                                                                                                                 >  1:3829582/1‑149 (MQ=255)
 gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATg                                                                                                                                                <  1:2201676/149‑1 (MQ=255)
   aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGt                                                                                                                                              >  1:89891/1‑149 (MQ=255)
   aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGt                                                                                                                                              >  2:3526024/1‑149 (MQ=255)
   aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGt                                                                                                                                              <  1:2120292/149‑1 (MQ=255)
         ggATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTc                                                                                                                                        <  1:223561/149‑1 (MQ=255)
           aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCtt                                                                                                                                      <  2:3744165/149‑1 (MQ=255)
            ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACTATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                     <  1:2308805/149‑1 (MQ=255)
            ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                     <  1:232962/149‑1 (MQ=255)
            ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                     <  2:684734/149‑1 (MQ=255)
              aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTg                                                                                                                                   >  1:4069341/1‑149 (MQ=255)
                 tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCa                                                                                                                                <  2:4017218/149‑1 (MQ=255)
                   aCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAgg                                                                                                                              <  2:85934/149‑1 (MQ=255)
                         cACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGc                                                                                                                        <  2:711869/149‑1 (MQ=255)
                           cGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAg                                                                                                                      >  1:581468/1‑149 (MQ=255)
                            gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGt                                                                                                                     <  2:332088/149‑1 (MQ=255)
                               ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGt                                                                                                                  >  2:3622927/1‑149 (MQ=255)
                               ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGt                                                                                                                  >  2:3248980/1‑149 (MQ=255)
                                 ttACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGTCAGGTCGAGCAGTTGTCg                                                                                                                >  1:1567116/1‑149 (MQ=255)
                                 ttACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCg                                                                                                                >  2:83360/1‑149 (MQ=255)
                                  tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGa                                                                                                               >  2:2415614/1‑149 (MQ=255)
                                         aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                        >  1:2968425/1‑149 (MQ=255)
                                         aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                        >  2:2040207/1‑149 (MQ=255)
                                         aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                        >  1:1227848/1‑149 (MQ=255)
                                          tAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGagag                                                                                                       <  1:3196259/149‑1 (MQ=255)
                                            aCCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTNGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtg                                                                                                     <  1:2085621/149‑1 (MQ=255)
                                             ccTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                    <  2:1958962/149‑1 (MQ=255)
                                               tGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgact                                                                                                  >  1:1587764/1‑149 (MQ=255)
                                               tGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgact                                                                                                  >  1:2126643/1‑149 (MQ=255)
                                                ggCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactg                                                                                                 >  1:4038158/1‑149 (MQ=255)
                                                 gCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactga                                                                                                <  2:932409/149‑1 (MQ=255)
                                                 gCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactga                                                                                                <  1:2724410/149‑1 (MQ=255)
                                                  cTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactgac                                                                                               >  1:2413227/1‑149 (MQ=255)
                                                     aCCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTGTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaa                                                                                            >  1:3498026/1‑149 (MQ=255)
                                                      ccGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaaa                                                                                           >  2:2801358/1‑149 (MQ=255)
                                                       cGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaaaa                                                                                          <  2:1876398/149‑1 (MQ=255)
                                                       cGAGTTGCAGACGAAAATCCTGCGCACTCTGGGAAAGACGATCAGCGCACTTTTGTAACTGCTGTGTGTCTTTTAATTCGGTGCGTTAGGCTATGGTGATGTCTTTTGACAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaaaa                                                                                          >  2:2919552/1‑149 (MQ=255)
                                                           ttGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATAct                                                                                      >  2:1123695/1‑149 (MQ=255)
                                                            tGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActg                                                                                     <  1:1939263/149‑1 (MQ=255)
                                                             gCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActgg                                                                                    <  1:443073/149‑1 (MQ=255)
                                                              cAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggc                                                                                   <  1:782546/149‑1 (MQ=255)
                                                              cAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggc                                                                                   <  1:755991/149‑1 (MQ=255)
                                                               aGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggct                                                                                  >  1:2267373/1‑149 (MQ=255)
                                                                 aCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggctgg                                                                                <  1:2437426/149‑1 (MQ=255)
                                                                 aCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggctgg                                                                                >  1:3858298/1‑149 (MQ=255)
                                                                  cGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACtggctggc                                                                               >  1:532432/1‑149 (MQ=255)
                                                                    aaaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                             >  1:3074108/1‑149 (MQ=255)
                                                                       aTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGt                                                                          <  1:2744964/149‑1 (MQ=255)
                                                                        tCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTc                                                                         >  2:2169669/1‑149 (MQ=255)
                                                                        tCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTc                                                                         <  2:91366/149‑1 (MQ=255)
                                                                                    tGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCtggtg                                                             <  2:3594630/149‑1 (MQ=255)
                                                                                     gggCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCtggtgg                                                            >  1:4076139/1‑149 (MQ=255)
                                                                                     gggCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCtggtgg                                                            >  1:2487380/1‑149 (MQ=255)
                                                                                     gggCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCtggtgg                                                            <  1:2843785/149‑1 (MQ=255)
                                                                                       gCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGc                                                          <  1:1584556/149‑1 (MQ=255)
                                                                                       gCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGc                                                          <  2:311024/149‑1 (MQ=255)
                                                                                         aaGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCtt                                                        <  2:3007681/149‑1 (MQ=255)
                                                                                          aGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtc                                                       <  2:2329071/149‑1 (MQ=255)
                                                                                           gACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtcg                                                      <  1:3849321/149‑1 (MQ=255)
                                                                                             cGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtcgtc                                                    <  2:66165/149‑1 (MQ=255)
                                                                                               aTCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGa                                                  >  2:714679/1‑149 (MQ=255)
                                                                                               aTCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGa                                                  >  2:2382767/1‑149 (MQ=255)
                                                                                                  aGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGa                                               <  1:3332329/149‑1 (MQ=255)
                                                                                                  aGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGa                                               <  2:2562056/149‑1 (MQ=255)
                                                                                                   gcgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAt                                              >  1:808020/1‑149 (MQ=255)
                                                                                                   gcgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAt                                              >  2:3916055/1‑149 (MQ=255)
                                                                                                     gcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCa                                            <  1:3851518/149‑1 (MQ=255)
                                                                                                        cTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGa                                         <  1:1448575/149‑1 (MQ=255)
                                                                                                        cTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGa                                         <  1:1033469/149‑1 (MQ=255)
                                                                                                                  tgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTc                               <  2:338178/149‑1 (MQ=255)
                                                                                                                   gctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTcc                              >  1:3799829/1‑149 (MQ=255)
                                                                                                                    ctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCg                             <  2:1902909/149‑1 (MQ=255)
                                                                                                                     tgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAGGATCATGACGTCCGCTTCCgg                            >  1:1804615/1‑149 (MQ=255)
                                                                                                                     tgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCgg                            >  1:2525070/1‑149 (MQ=255)
                                                                                                                     tgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCgg                            >  1:2013101/1‑149 (MQ=255)
                                                                                                                       ggTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGa                          >  2:1165645/1‑149 (MQ=255)
                                                                                                                        gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGAt                         >  1:3659293/1‑149 (MQ=255)
                                                                                                                        gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGAt                         >  2:512550/1‑149 (MQ=255)
                                                                                                                          gtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCa                       >  2:6759/1‑149 (MQ=255)
                                                                                                                          gtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCa                       >  2:1729682/1‑149 (MQ=255)
                                                                                                                           tCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAg                      <  2:4062352/149‑1 (MQ=255)
                                                                                                                            cTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGt                     <  1:2585933/149‑1 (MQ=255)
                                                                                                                             ttttAATTCGGTGCGGTAGGCTTTGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCGGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGtt                    <  2:1087746/149‑1 (MQ=255)
                                                                                                                             ttttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCAGCTTCCGGGATCAGtt                    >  1:3315238/1‑149 (MQ=255)
                                                                                                                              tttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTc                   >  1:2657318/1‑149 (MQ=255)
                                                                                                                               ttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCg                  >  2:260210/1‑149 (MQ=255)
                                                                                                                                 aaTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGcc                <  1:3622927/149‑1 (MQ=255)
                                                                                                                                 aaTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGcc                <  2:912816/149‑1 (MQ=255)
                                                                                                                                  aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCa               >  1:696746/1‑149 (MQ=255)
                                                                                                                                  aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCa               >  2:1785301/1‑149 (MQ=255)
                                                                                                                                   ttCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCaa              >  1:3184871/1‑149 (MQ=255)
                                                                                                                                      ggTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATc           >  2:2139392/1‑149 (MQ=255)
                                                                                                                                        tGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCa         <  1:3248980/149‑1 (MQ=255)
                                                                                                                                         gCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCAc        <  1:1213018/149‑1 (MQ=255)
                                                                                                                                            gTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct     >  2:1254219/1‑149 (MQ=255)
                                                                                                                                            gTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct     >  1:1510616/1‑149 (MQ=255)
                                                                                                                                              aGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct   <  1:540522/149‑1 (MQ=255)
                                                                                                                                              aGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct   >  1:1164389/1‑149 (MQ=255)
                                                                                                                                              aGGCTATGGTGATGTATTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGATTCCGGGATCAGTTCGCCAAATCCActctct   >  2:1866002/1‑149 (MQ=255)
                                                                                                                                               ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCtt  >  2:1482720/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT  >  CP009273/1886579‑1886870

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