Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 521,011 G→T 57.0% E765D (GAG→GAT rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273521,0110GT57.0% ‑0.6 / inf 312E765D (GAG→GATrhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (66/68);  new base T (94/84);  total (160/152)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 320 total aligned reads displayed.

GGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGC  >  CP009273/520862‑521159
                                                                                                                                                     |                                                                                                                                                    
 gCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGgaga                                                                                                                                                     <  1:3650760/150‑1 (MQ=255)
   gAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACg                                                                                                                                                   >  1:1907643/1‑150 (MQ=255)
   gAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACg                                                                                                                                                   >  2:1903098/1‑150 (MQ=255)
    aaCCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACgg                                                                                                                                                  <  2:663799/150‑1 (MQ=255)
    aaCCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACgg                                                                                                                                                  >  2:3818622/1‑150 (MQ=255)
    aaCCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACgg                                                                                                                                                  <  1:2796265/150‑1 (MQ=255)
       cGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACggggg                                                                                                                                               >  1:429428/1‑150 (MQ=255)
        ggCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGa                                                                                                                                              <  2:493160/150‑1 (MQ=255)
        ggCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTCTGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGa                                                                                                                                              >  1:3753450/1‑150 (MQ=255)
        ggCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGa                                                                                                                                              <  1:1228413/150‑1 (MQ=255)
        ggCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGa                                                                                                                                              >  2:1738269/1‑150 (MQ=255)
          cAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAAc                                                                                                                                            <  1:291277/150‑1 (MQ=255)
            gagCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAActg                                                                                                                                          >  2:152140/1‑150 (MQ=255)
             agCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGACctgc                                                                                                                                         >  2:1019968/1‑150 (MQ=255)
                   ggCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGgc                                                                                                                                   <  1:767141/150‑1 (MQ=255)
                      aGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGgcagc                                                                                                                                <  1:3334260/150‑1 (MQ=255)
                         atgatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATg                                                                                                                             >  2:2652323/1‑150 (MQ=255)
                         atgatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATg                                                                                                                             >  1:4211252/1‑150 (MQ=255)
                         atgatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATg                                                                                                                             <  2:2056935/150‑1 (MQ=255)
                          tgatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATga                                                                                                                            <  2:1248278/150‑1 (MQ=255)
                           gatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATgag                                                                                                                           <  1:4479010/150‑1 (MQ=255)
                           gatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATgag                                                                                                                           <  1:1305644/150‑1 (MQ=255)
                             tgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGAc                                                                                                                         <  1:3248527/150‑1 (MQ=255)
                              ggCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACg                                                                                                                        >  2:1994534/1‑150 (MQ=255)
                              ggCCACGGCTGGCTGACAGACATCAGCCAACTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACg                                                                                                                        >  1:2025871/1‑150 (MQ=255)
                                ccACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGaa                                                                                                                      <  2:1421967/150‑1 (MQ=255)
                                 cACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGaaa                                                                                                                     <  1:4552747/150‑1 (MQ=255)
                                   cGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAaca                                                                                                                   <  1:312498/150‑1 (MQ=255)
                                    ggctggctGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAacac                                                                                                                  <  2:3112237/150‑1 (MQ=255)
                                         gctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATa                                                                                                             <  1:3819501/150‑1 (MQ=255)
                                                     cAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAggg                                                                                                 >  1:1698954/1‑150 (MQ=255)
                                                      aGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAgggg                                                                                                <  2:2271819/150‑1 (MQ=255)
                                                      aGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAgggg                                                                                                >  1:2422232/1‑150 (MQ=255)
                                                         cACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTg                                                                                             >  1:4509896/1‑150 (MQ=255)
                                                         cACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTg                                                                                             >  1:2042277/1‑150 (MQ=255)
                                                         cACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTg                                                                                             >  1:1054952/1‑150 (MQ=255)
                                                         cACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTg                                                                                             >  2:3373210/1‑150 (MQ=255)
                                                            cTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCa                                                                                          <  2:960629/150‑1 (MQ=255)
                                                               aGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAAc                                                                                       <  1:3695444/150‑1 (MQ=255)
                                                                 cGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcg                                                                                     >  2:1689930/1‑150 (MQ=255)
                                                                  gAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgc                                                                                    >  2:3741555/1‑150 (MQ=255)
                                                                   aaGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgcg                                                                                   >  1:2627293/1‑150 (MQ=255)
                                                                        cACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACg                                                                              <  1:2201661/150‑1 (MQ=255)
                                                                         aCCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGc                                                                             <  2:984055/150‑1 (MQ=255)
                                                                            gtgtTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCgg                                                                          <  1:1857714/150‑1 (MQ=255)
                                                                               ttGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACa                                                                       <  1:1184159/150‑1 (MQ=255)
                                                                               ttGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACa                                                                       <  1:734625/150‑1 (MQ=255)
                                                                                tGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                      <  1:1014147/150‑1 (MQ=255)
                                                                                  ccGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGcc                                                                    <  2:3014116/150‑1 (MQ=255)
                                                                                  ccGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGcc                                                                    >  1:651518/1‑150 (MQ=255)
                                                                                   cGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCt                                                                   >  1:1607296/1‑150 (MQ=255)
                                                                                    gTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCGGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTg                                                                  >  1:2348281/1‑150 (MQ=255)
                                                                                      ccACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTgcc                                                                >  2:627591/1‑150 (MQ=255)
                                                                                            tGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGt                                                          >  2:2999063/1‑150 (MQ=255)
                                                                                            tGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGt                                                          >  2:1847108/1‑150 (MQ=255)
                                                                                             ggCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTg                                                         <  1:3555443/150‑1 (MQ=255)
                                                                                             ggCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTg                                                         <  1:976752/150‑1 (MQ=255)
                                                                                               cTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGa                                                       <  1:4329989/150‑1 (MQ=255)
                                                                                                 aTGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAGCCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGt                                                     >  1:3444905/1‑150 (MQ=255)
                                                                                                   gACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTgg                                                   <  1:4127158/150‑1 (MQ=255)
                                                                                                     cGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGTGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCt                                                 <  1:2321823/150‑1 (MQ=255)
                                                                                                     cGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCt                                                 >  1:2626462/1‑150 (MQ=255)
                                                                                                      gATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTg                                                >  2:435798/1‑150 (MQ=255)
                                                                                                      caTAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTg                                                >  2:1061439/2‑150 (MQ=255)
                                                                                                            ggCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAt                                          <  1:968694/150‑1 (MQ=255)
                                                                                                            ggCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAt                                          <  2:59838/150‑1 (MQ=255)
                                                                                                            ggCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAt                                          >  2:3377218/1‑150 (MQ=255)
                                                                                                             gccgccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATg                                         <  1:530192/150‑1 (MQ=255)
                                                                                                              ccgccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATgg                                        >  1:2285741/1‑150 (MQ=255)
                                                                                                              ccgccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATgg                                        >  1:164856/1‑150 (MQ=255)
                                                                                                              ccgccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATgg                                        >  2:692906/1‑150 (MQ=255)
                                                                                                                gccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCa                                      <  1:2752604/150‑1 (MQ=255)
                                                                                                                 ccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAg                                     <  1:1550692/150‑1 (MQ=255)
                                                                                                                 ccTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAg                                     <  1:1738269/150‑1 (MQ=255)
                                                                                                                  cTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGc                                    >  1:4593469/1‑150 (MQ=255)
                                                                                                                   tGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCg                                   >  2:2889676/1‑150 (MQ=255)
                                                                                                                   tGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCg                                   >  1:3607881/1‑150 (MQ=255)
                                                                                                                   tGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCg                                   >  1:2068785/1‑150 (MQ=255)
                                                                                                                   tGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCg                                   <  2:94463/150‑1 (MQ=255)
                                                                                                                     aCCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGt                                 >  2:3679792/1‑150 (MQ=255)
                                                                                                                      ccGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGtt                                <  1:1417882/150‑1 (MQ=255)
                                                                                                                      ccGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGtt                                <  1:1617492/150‑1 (MQ=255)
                                                                                                                        ggCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTAc                              >  1:3644146/1‑150 (MQ=255)
                                                                                                                          cGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCt                            <  2:4413605/150‑1 (MQ=255)
                                                                                                                          cGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCt                            <  2:113989/150‑1 (MQ=255)
                                                                                                                           gAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTg                           <  2:807720/150‑1 (MQ=255)
                                                                                                                            aaTGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTgg                          >  1:2228108/1‑150 (MQ=255)
                                                                                                                             aTGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGc                         >  1:854207/1‑150 (MQ=255)
                                                                                                                             aTGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGc                         <  1:1106743/150‑1 (MQ=255)
                                                                                                                                 cAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGa                     >  2:1626915/1‑150 (MQ=255)
                                                                                                                                  aGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGaa                    <  2:2419619/150‑1 (MQ=255)
                                                                                                                                     cGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGa                 <  2:3870981/150‑1 (MQ=255)
                                                                                                                                         ggAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCt             <  2:1338685/150‑1 (MQ=255)
                                                                                                                                            gaACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTggg          >  1:3009492/1‑150 (MQ=255)
                                                                                                                                             aaCCCGGATACGGGGGAACAGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACTCCCGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGc         >  2:3536008/1‑150 (MQ=255)
                                                                                                                                                 cGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGa     <  2:2854549/150‑1 (MQ=255)
                                                                                                                                                  ggATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAc    >  2:4152864/1‑150 (MQ=255)
                                                                                                                                                  ggATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAc    >  2:2569018/1‑150 (MQ=255)
                                                                                                                                                   gagaCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAcg   >  1:61552/1‑150 (MQ=255)
                                                                                                                                                    agaCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAcgc  >  1:1812896/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGC  >  CP009273/520862‑521159

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: