Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 280,145 G→A 21.5% E412K (GAG→AAG)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273280,1450GA21.5% 1293.7 / inf 564E412K (GAG→AAG) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (189/253);  new base A (58/63);  total (247/317)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 573 total aligned reads displayed.

CGGGCGGCGTGCCGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTG  >  CP009273/279996‑280294
                                                                                                                                                     |                                                                                                                                                     
cgggcggCGTGCCGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGcgcg                                                                                                                                                       >  2:2045169/1‑150 (MQ=255)
  ggcggcGTGCCGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAg                                                                                                                                                     <  1:279706/150‑1 (MQ=255)
        gTGCCGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGAc                                                                                                                                               <  2:1610613/150‑1 (MQ=255)
           ccGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGc                                                                                                                                            <  1:4367823/150‑1 (MQ=255)
            cGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCg                                                                                                                                           <  1:1317053/150‑1 (MQ=255)
             ggAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGt                                                                                                                                          >  2:2379699/1‑150 (MQ=255)
                 gTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTagag                                                                                                                                      <  1:1225791/150‑1 (MQ=255)
                 gTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTagag                                                                                                                                      >  1:2519626/1‑150 (MQ=255)
                    aTGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCg                                                                                                                                   >  2:2304496/1‑150 (MQ=255)
                    aTGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCg                                                                                                                                   >  2:2702338/1‑150 (MQ=255)
                     tGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCgg                                                                                                                                  <  2:4840517/150‑1 (MQ=255)
                     tGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCgg                                                                                                                                  <  2:52487/150‑1 (MQ=255)
                       cTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGAt                                                                                                                                >  2:256002/1‑150 (MQ=255)
                              tGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGa                                                                                                                         <  2:2072472/150‑1 (MQ=255)
                                 gcgACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGAtac                                                                                                                      <  2:2235547/150‑3 (MQ=255)
                                  cgACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCt                                                                                                                     <  1:4228836/150‑1 (MQ=255)
                                   gACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTg                                                                                                                    <  1:1639930/150‑1 (MQ=255)
                                   gACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTg                                                                                                                    >  2:3209032/1‑150 (MQ=255)
                                    aCCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTgc                                                                                                                   >  2:4398582/1‑150 (MQ=255)
                                     ccTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGATTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTgcc                                                                                                                  <  2:415230/150‑1 (MQ=255)
                                     ccTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTgcc                                                                                                                  <  1:2243808/150‑1 (MQ=255)
                                     ccTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTgcc                                                                                                                  <  1:4708467/150‑1 (MQ=255)
                                      cTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTgccg                                                                                                                 <  2:4104862/150‑1 (MQ=255)
                                         ggCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGccgccg                                                                                                              >  1:1905784/1‑150 (MQ=255)
                                          gCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCgg                                                                                                             <  1:3645232/150‑1 (MQ=255)
                                            ctgctgCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGgag                                                                                                           >  1:2426758/1‑150 (MQ=255)
                                             tgctgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGgaga                                                                                                          <  1:823729/150‑1 (MQ=255)
                                                tgcATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAgg                                                                                                       >  2:1744701/1‑150 (MQ=255)
                                                    tCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCaaa                                                                                                   >  1:2960555/1‑150 (MQ=255)
                                                      tGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAg                                                                                                 >  1:3743571/1‑150 (MQ=255)
                                                      tGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAg                                                                                                 <  2:906503/150‑1 (MQ=255)
                                                       ggACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGc                                                                                                <  2:1127860/150‑1 (MQ=255)
                                                         aCGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGa                                                                                              <  1:4801583/150‑1 (MQ=255)
                                                          cGCCTTGACCGTGACCGGCCAGTCGGTGGGCGAGAACATTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGaa                                                                                             <  2:3830856/150‑1 (MQ=255)
                                                          cGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGaa                                                                                             <  1:2551422/150‑1 (MQ=255)
                                                             cATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAgg                                                                                          <  1:3366079/150‑1 (MQ=255)
                                                                 accgtgacggGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGGGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTg                                                                                      <  1:4431935/150‑1 (MQ=255)
                                                                 accgtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTg                                                                                      <  2:4578867/150‑1 (MQ=255)
                                                                  ccgtgaccgGCCAGACGGGGGGCGAGAACCTTGATTGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGa                                                                                     <  1:2929920/150‑1 (MQ=255)
                                                                   cgtgaccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGAc                                                                                    <  2:4855556/150‑1 (MQ=255)
                                                                        ccgGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGa                                                                               <  2:760322/150‑1 (MQ=255)
                                                                          ggCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACg                                                                             <  1:309194/150‑1 (MQ=255)
                                                                          ggCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACg                                                                             <  1:1780768/150‑1 (MQ=255)
                                                                          ggCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACg                                                                             >  2:2887414/1‑150 (MQ=255)
                                                                           gCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACgg                                                                            <  1:2702338/150‑1 (MQ=255)
                                                                                aCGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGc                                                                       <  1:1044091/150‑1 (MQ=255)
                                                                                aCGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGc                                                                       <  1:4127528/150‑1 (MQ=255)
                                                                                 cGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCt                                                                      <  2:1905784/150‑1 (MQ=255)
                                                                                      ggCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCg                                                                 <  2:310505/150‑1 (MQ=255)
                                                                                      ggCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCg                                                                 <  1:2122830/150‑1 (MQ=255)
                                                                                        cGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGAc                                                               <  1:2489039/150‑1 (MQ=255)
                                                                                         gagaACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACg                                                              >  2:3594475/1‑150 (MQ=255)
                                                                                         gagaACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACg                                                              >  1:3539586/1‑150 (MQ=255)
                                                                                         gagaACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACg                                                              <  2:3169102/150‑1 (MQ=255)
                                                                                            aaCCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGc                                                           >  2:1288847/1‑150 (MQ=255)
                                                                                               cTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAAc                                                        >  2:3324778/1‑150 (MQ=255)
                                                                                               cTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAAc                                                        >  1:1287131/1‑150 (MQ=255)
                                                                                                 tGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACAt                                                      >  2:1477294/1‑150 (MQ=255)
                                                                                                     tggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCt                                                  <  1:3954468/150‑1 (MQ=255)
                                                                                                     tggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCt                                                  >  2:3730750/1‑150 (MQ=255)
                                                                                                      ggtggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTc                                                 >  1:1441059/1‑150 (MQ=255)
                                                                                                        tggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCg                                               >  2:3501205/1‑150 (MQ=255)
                                                                                                         ggCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCgg                                              >  1:1387165/1‑150 (MQ=255)
                                                                                                           cAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACAGCGCTCCGGac                                            >  2:3249231/1‑149 (MQ=255)
                                                                                                            aGGCGTCCGAGCGCCGGGCGCGCTTCCGCCATTGCCTGCGCGAGCAGGACGGCTTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAg                                           <  2:1201665/150‑1 (MQ=255)
                                                                                                            aGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAg                                           <  2:3994991/150‑1 (MQ=255)
                                                                                                            aGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAg                                           >  2:587966/1‑150 (MQ=255)
                                                                                                             ggCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAgg                                          <  2:435406/150‑1 (MQ=255)
                                                                                                              gCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGt                                         <  2:272307/150‑1 (MQ=255)
                                                                                                              gCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGt                                         <  2:4554921/150‑1 (MQ=255)
                                                                                                                gTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTc                                       >  2:2348275/1‑150 (MQ=255)
                                                                                                                gTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTc                                       >  1:526953/1‑150 (MQ=255)
                                                                                                                 tCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCg                                      <  1:2009879/150‑1 (MQ=255)
                                                                                                                 tCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCg                                      <  1:2074569/150‑1 (MQ=255)
                                                                                                                 tCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCg                                      <  1:3209032/150‑1 (MQ=255)
                                                                                                                   cGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGt                                    >  2:2122925/1‑150 (MQ=255)
                                                                                                                   cGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGt                                    <  2:4046697/150‑1 (MQ=255)
                                                                                                                    gAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                   <  2:3390769/150‑1 (MQ=255)
                                                                                                                    gAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                   <  1:1099871/150‑1 (MQ=255)
                                                                                                                    gAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTg                                   <  2:2757122/150‑1 (MQ=255)
                                                                                                                     aGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGa                                  <  2:2632098/150‑1 (MQ=255)
                                                                                                                           gggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAgg                            <  1:1744644/150‑1 (MQ=255)
                                                                                                                            ggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGc                           <  1:2628844/150‑1 (MQ=255)
                                                                                                                            ggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGc                           >  2:4432197/1‑150 (MQ=255)
                                                                                                                            ggCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGc                           <  1:604653/150‑1 (MQ=255)
                                                                                                                                cgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACg                       <  1:2332746/150‑1 (MQ=255)
                                                                                                                                cgcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACg                       >  2:4276864/1‑150 (MQ=255)
                                                                                                                                 gcTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACgg                      >  2:2996905/1‑150 (MQ=255)
                                                                                                                                  cTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGc                     <  1:4345148/150‑1 (MQ=255)
                                                                                                                                  cTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGACTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGc                     >  1:3494971/1‑150 (MQ=255)
                                                                                                                                   ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                    >  1:1672692/1‑150 (MQ=255)
                                                                                                                                   ttCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCg                    >  1:1502137/1‑150 (MQ=255)
                                                                                                                                     ccGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGAt                  <  2:1950271/150‑1 (MQ=255)
                                                                                                                                      cGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATc                 <  1:2578362/150‑1 (MQ=255)
                                                                                                                                      cGCCAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATc                 <  1:3854951/150‑1 (MQ=255)
                                                                                                                                        ccAGTGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGa               >  1:821787/1‑150 (MQ=255)
                                                                                                                                         cAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAc              >  2:1591094/1‑150 (MQ=255)
                                                                                                                                          aGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAcc             >  1:4297137/1‑150 (MQ=255)
                                                                                                                                            tGCCTGCGCAAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGAAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACTTCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCg           <  1:4876025/150‑1 (MQ=255)
                                                                                                                                                     gAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAACGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGgtggtg  <  1:954954/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGCGGCGTGCCGGAGGTGATGCTCCACCTGCGCGACCTCGGCCTGCTGCATCTGGACGCCATGACCGTGACCGGCCAGACGGTGGGCGAGAACCTTGAATGGTGGCAGGCGTCCGAGCGCCGGGCGCGCTTCCGCCAGTGCCTGCGCGAGCAGGACGGCGTAGAGCCGGATGACGTGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTG  >  CP009273/279996‑280294

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: