Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,243,278 G→A 18.2% S431S (AGC→AGT dhaM ← putative dihydroxyacetone‑specific PTS enzymes: HPr component/EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,243,2780GA18.2% 871.2 / 166.4 336S431S (AGC→AGTdhaMputative dihydroxyacetone‑specific PTS enzymes: HPr component/EI component
Reads supporting (aligned to +/- strand):  ref base G (126/148);  new base A (28/33);  total (154/182)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

CGATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTG  >  CP009273/1243131‑1243427
                                                                                                                                                   |                                                                                                                                                     
 gATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCtg                                                                                                                                                    >  2:434567/1‑150 (MQ=255)
  aTGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCtgt                                                                                                                                                   <  1:788932/150‑1 (MQ=255)
    gCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTgg                                                                                                                                                 <  1:4153963/150‑1 (MQ=255)
     ccATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTATGGCACTGTggg                                                                                                                                                >  2:2977869/1‑150 (MQ=255)
     ccATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTggg                                                                                                                                                >  2:1192677/1‑150 (MQ=255)
        tCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATa                                                                                                                                             >  1:384679/1‑150 (MQ=255)
         ccGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGAGTAGGGCGCTGTGGGAtcc                                                                                                                                            >  2:5496197/1‑148 (MQ=255)
           gAACAGTGGCTTAACCCTGTCGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTTGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCg                                                                                                                                          <  1:2703295/150‑1 (MQ=255)
            aaCAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCgg                                                                                                                                         >  1:1420293/1‑150 (MQ=255)
            aaCAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCgg                                                                                                                                         >  2:173052/1‑150 (MQ=255)
              cAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGAc                                                                                                                                       <  1:118021/150‑1 (MQ=255)
               aGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACt                                                                                                                                      >  1:3520106/1‑150 (MQ=255)
                gTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACtt                                                                                                                                     >  2:2676715/1‑150 (MQ=255)
                  ggCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTcc                                                                                                                                   <  2:373461/150‑1 (MQ=255)
                  ggCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTcc                                                                                                                                   >  1:4628556/1‑150 (MQ=255)
                  ggCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGAGTAGGGCGCTGTGGGATACCGGACTTcc                                                                                                                                   >  2:4820001/1‑150 (MQ=255)
                   gCTTAACCCAGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCg                                                                                                                                  <  2:4274876/150‑1 (MQ=255)
                      tAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGgcg                                                                                                                               >  2:1011731/1‑150 (MQ=255)
                        aCCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCt                                                                                                                             >  2:3352218/1‑150 (MQ=255)
                         cccTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTa                                                                                                                            >  1:5458998/1‑150 (MQ=255)
                          ccTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTaa                                                                                                                           <  1:4769351/150‑1 (MQ=255)
                           cTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAg                                                                                                                          <  2:2359093/150‑1 (MQ=255)
                                  ttGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCACATCCAGCCAATCCCCAGTTCACGGGCGCTTAGGGCGCTGGGGGATACCGGACTTCCGGCGCTAAGGCCGATa                                                                                                                   >  2:3195634/1‑150 (MQ=255)
                                   tGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATAc                                                                                                                  >  2:5277076/1‑150 (MQ=255)
                                   tGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCGTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATAc                                                                                                                  <  1:2272418/150‑1 (MQ=255)
                                         gTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTa                                                                                                            >  1:3749316/1‑150 (MQ=255)
                                         gTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTa                                                                                                            >  2:1947973/1‑150 (MQ=255)
                                          ttGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTAc                                                                                                           >  1:2728492/1‑150 (MQ=255)
                                             cGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAAc                                                                                                        <  2:1081658/150‑1 (MQ=255)
                                              gatttAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGGTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAAcc                                                                                                       <  2:4385061/148‑1 (MQ=255)
                                               ttttAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAAccg                                                                                                      <  2:4070560/150‑1 (MQ=255)
                                                  tAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATTCCTTTTACAAccgccg                                                                                                   <  1:4031648/150‑1 (MQ=255)
                                                   aaCGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCgg                                                                                                  <  2:288525/150‑1 (MQ=255)
                                                   aaCGTCCAGCGTGAGCGTTGCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCgg                                                                                                  <  2:3897183/150‑1 (MQ=255)
                                                     cGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGAt                                                                                                >  1:2012566/1‑150 (MQ=255)
                                                      gTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATc                                                                                               >  1:1437745/1‑150 (MQ=255)
                                                         cAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAg                                                                                            >  1:4973660/1‑150 (MQ=255)
                                                           gcgttagcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGtt                                                                                          <  1:4091175/150‑1 (MQ=255)
                                                            cgttagcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTg                                                                                         <  1:1746437/150‑1 (MQ=255)
                                                             gttagcgttTCTTCGGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGc                                                                                        <  2:3918397/150‑1 (MQ=255)
                                                                   gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCCGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  <  2:120255/150‑1 (MQ=255)
                                                                   gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  >  2:4168760/1‑150 (MQ=255)
                                                                   gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  >  2:5018914/1‑150 (MQ=255)
                                                                   gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  <  2:1464805/150‑1 (MQ=255)
                                                                   gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  <  1:4550073/150‑1 (MQ=255)
                                                                   gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  <  1:5417874/150‑1 (MQ=255)
                                                                      tcttctGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACAGCCGGATCCAGTTGCAGTACtgtg                                                                               >  1:345941/1‑150 (MQ=255)
                                                                           tGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAgg                                                                          <  1:1003409/150‑1 (MQ=255)
                                                                             gTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGAt                                                                        <  1:4897978/150‑1 (MQ=255)
                                                                             gTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGAt                                                                        <  1:54722/150‑1 (MQ=255)
                                                                               tGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATaa                                                                      >  1:3959575/1‑150 (MQ=255)
                                                                                 tATCTCATACAGGTTCTCACCCTGCTGGCAAATCCAGCCAAGCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCTGACTTCCGGCGCTAAGGCAGATAACTTTTACAACCGCCGGATCCAGGTGCAGTACTGTGGAAGGATAAAt                                                                    <  2:5388520/150‑1 (MQ=255)
                                                                                   tCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGt                                                                  <  2:1382891/150‑1 (MQ=255)
                                                                                    cGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGtt                                                                 >  1:599118/1‑150 (MQ=255)
                                                                                    cGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGtt                                                                 <  2:3404312/150‑1 (MQ=255)
                                                                                    cGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGtt                                                                 >  1:1778993/1‑150 (MQ=255)
                                                                                      cATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtct                                                               >  1:936348/1‑150 (MQ=255)
                                                                                      cATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtct                                                               <  1:4150723/150‑1 (MQ=255)
                                                                                      cATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtct                                                               >  1:3781945/1‑150 (MQ=255)
                                                                                      cATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtct                                                               <  2:2864476/150‑1 (MQ=255)
                                                                                        tACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTcc                                                             <  1:1908051/150‑1 (MQ=255)
                                                                                        tACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTcc                                                             <  2:5163689/150‑1 (MQ=255)
                                                                                          cAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGc                                                           >  1:86551/1‑150 (MQ=255)
                                                                                           aGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGcc                                                          <  1:2413619/150‑1 (MQ=255)
                                                                                            gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCa                                                         <  2:5452691/150‑1 (MQ=255)
                                                                                            gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCa                                                         >  2:3556736/1‑150 (MQ=255)
                                                                                            gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCa                                                         <  2:4476034/150‑1 (MQ=255)
                                                                                            gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCa                                                         <  2:2642390/150‑1 (MQ=255)
                                                                                            gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCa                                                         <  2:3151224/150‑1 (MQ=255)
                                                                                             tttCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAg                                                        <  1:2335496/150‑1 (MQ=255)
                                                                                               tctcACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTa                                                      <  2:493630/150‑1 (MQ=255)
                                                                                                ctcACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAg                                                     >  1:3165251/1‑150 (MQ=255)
                                                                                                   aCCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAAt                                                  >  1:4447498/1‑150 (MQ=255)
                                                                                                   aCCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAAt                                                  >  2:3903276/1‑150 (MQ=255)
                                                                                                    cccTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATa                                                 <  1:4862978/150‑1 (MQ=255)
                                                                                                          tgGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGAACCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGc                                           >  1:4146316/1‑150 (MQ=255)
                                                                                                          tgGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGc                                           <  1:1947974/150‑1 (MQ=255)
                                                                                                              aaaTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGt                                       <  1:1584590/150‑1 (MQ=255)
                                                                                                              aaaTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGt                                       <  2:2465460/150‑1 (MQ=255)
                                                                                                              aaaTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGt                                       <  2:3233737/150‑1 (MQ=255)
                                                                                                              aaaTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGt                                       >  2:2226660/1‑150 (MQ=255)
                                                                                                                aTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTa                                     >  1:2158969/1‑150 (MQ=255)
                                                                                                                 tCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTaa                                    <  1:3546841/150‑1 (MQ=255)
                                                                                                                   cAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAAc                                  <  1:700318/150‑1 (MQ=255)
                                                                                                                    aGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACt                                 <  1:566354/150‑1 (MQ=255)
                                                                                                                        aaTCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCgg                             <  2:826480/150‑1 (MQ=255)
                                                                                                                        aaTCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCgg                             >  1:296764/1‑150 (MQ=255)
                                                                                                                          tCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGa                           >  1:1599476/1‑150 (MQ=255)
                                                                                                                               aGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttct                      <  2:4790217/150‑1 (MQ=255)
                                                                                                                               aGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttct                      <  2:5511836/150‑1 (MQ=255)
                                                                                                                               aGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttct                      >  2:302090/1‑150 (MQ=255)
                                                                                                                                gTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATCGTTGGCGAGTTAAACTGCGGGAGttctt                     >  2:4556810/1‑150 (MQ=255)
                                                                                                                                gTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttctt                     >  1:2082088/1‑150 (MQ=255)
                                                                                                                                 ttCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttcttc                    <  2:784197/150‑1 (MQ=255)
                                                                                                                                   cACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTcttctt                  <  1:3097273/150‑1 (MQ=255)
                                                                                                                                   cACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTcttctt                  >  2:1617491/1‑150 (MQ=255)
                                                                                                                                    aCGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCttt                 >  1:5266067/1‑150 (MQ=255)
                                                                                                                                        gCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGtt             <  1:3851189/150‑1 (MQ=255)
                                                                                                                                         cGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGttt            >  1:5383788/1‑150 (MQ=255)
                                                                                                                                           aTTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTgg          >  1:5179527/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTG  >  CP009273/1243131‑1243427

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: