Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 521,011 G→T 21.1% 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,0110GT21.1% 745.1 / 134.1 264E765D (GAG→GATrhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (132/76);  new base T (41/15);  total (173/91)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

GGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCC  >  CP009273/520862‑521160
                                                                                                                                                     |                                                                                                                                                     
ggCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGgag                                                                                                                                                       >  1:3635481/1‑150 (MQ=255)
ggCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGgag                                                                                                                                                       >  1:3780300/1‑150 (MQ=255)
 gCGAACCGGGTGAGCTGTTGCAGTTTGATGGCCTCGGCTGGCTGACAGTCATAATCCACCTGAGCGAATGCCACCGTGTTGACGTCCACTATGGATATGACGATAATGGAAGCCTGACCGGCGAATGCCTTACGGTGGAGAACCCGgaga                                                                                                                                                      <  1:4010548/150‑1 (MQ=21)
  cGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGAc                                                                                                                                                     >  2:504579/1‑150 (MQ=255)
  cGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTAGGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGAc                                                                                                                                                     <  2:700313/150‑1 (MQ=255)
   gAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACg                                                                                                                                                    >  1:1826859/1‑150 (MQ=255)
      ccGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCAAGACGGTGGAGAACCCGGAGACgggg                                                                                                                                                 <  1:2982608/150‑1 (MQ=255)
       cGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACggggg                                                                                                                                                >  1:16273/1‑150 (MQ=255)
           agagCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAAct                                                                                                                                            >  2:1199368/1‑150 (MQ=255)
             agCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAActgc                                                                                                                                          >  2:2763011/1‑150 (MQ=255)
              gCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAActgct                                                                                                                                         <  1:4159032/150‑1 (MQ=255)
               cAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAActgctg                                                                                                                                        >  2:1563497/1‑150 (MQ=255)
               cAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAActgctg                                                                                                                                        >  2:893085/1‑150 (MQ=255)
                  tGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTgg                                                                                                                                     >  1:4463452/1‑150 (MQ=255)
                   ggCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGgc                                                                                                                                    >  2:4358087/1‑150 (MQ=255)
                     cAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGgcag                                                                                                                                  >  2:4471710/1‑150 (MQ=255)
                           gatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATgag                                                                                                                            <  1:1683684/150‑1 (MQ=255)
                           gatgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATgag                                                                                                                            <  1:1353395/150‑1 (MQ=255)
                             tgGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGAc                                                                                                                          >  1:2440597/1‑150 (MQ=255)
                                 cACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGaaa                                                                                                                      >  1:3061485/1‑150 (MQ=255)
                                 cACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGaaa                                                                                                                      >  1:1846062/1‑150 (MQ=255)
                                    ggctggctGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAacac                                                                                                                   >  1:2775331/1‑150 (MQ=255)
                                    ggctggctGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAacac                                                                                                                   >  1:1316182/1‑150 (MQ=255)
                                       tggctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCa                                                                                                                >  2:4064995/1‑150 (MQ=255)
                                       tggctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCa                                                                                                                >  2:2478713/1‑150 (MQ=255)
                                        ggctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTCTGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCAt                                                                                                               >  2:2897010/1‑150 (MQ=255)
                                          ctgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATAc                                                                                                             >  2:4415836/1‑150 (MQ=255)
                                           tgACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACa                                                                                                            >  1:2620259/1‑150 (MQ=255)
                                                   aTCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAg                                                                                                    >  2:2245790/1‑150 (MQ=255)
                                                     cAGCCTCCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAggg                                                                                                  <  1:209854/150‑1 (MQ=255)
                                                      aGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAgggg                                                                                                 >  2:1314214/1‑150 (MQ=255)
                                                      aGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAgggg                                                                                                 >  2:4094339/1‑150 (MQ=255)
                                                      aGCCACCTGAGCGAAGGACACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAgggg                                                                                                 >  2:1906982/1‑150 (MQ=255)
                                                       gCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGc                                                                                                <  1:3212422/150‑1 (MQ=255)
                                                        ccACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCt                                                                                               >  1:4066986/1‑150 (MQ=255)
                                                               aGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAAc                                                                                        <  2:3326222/150‑1 (MQ=255)
                                                               aGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAAc                                                                                        >  1:2222929/1‑150 (MQ=255)
                                                               aGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAAc                                                                                        <  1:226888/150‑1 (MQ=255)
                                                                gCGAAGGCCACCGTGTTGCCGTCCACTAGGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAAcc                                                                                       <  2:4456685/150‑1 (MQ=255)
                                                                 cGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcg                                                                                      >  1:1990522/1‑150 (MQ=255)
                                                                 cGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcg                                                                                      <  2:4601189/150‑1 (MQ=255)
                                                                 cGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcg                                                                                      >  2:1666636/1‑150 (MQ=255)
                                                                 cGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcg                                                                                      >  2:1934650/1‑150 (MQ=255)
                                                                  gAAGGCCATCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgc                                                                                     >  2:363933/1‑150 (MQ=255)
                                                                  gAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgc                                                                                     >  1:1951342/1‑150 (MQ=255)
                                                                  gAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgc                                                                                     >  2:4233146/1‑150 (MQ=255)
                                                                   aaGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgcg                                                                                    >  1:423857/1‑150 (MQ=255)
                                                                   aaGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACcgcg                                                                                    >  1:507120/1‑150 (MQ=255)
                                                                      gCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCa                                                                                 >  2:1117583/1‑150 (MQ=255)
                                                                      gCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCa                                                                                 <  1:3432508/150‑1 (MQ=255)
                                                                          ccGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGcc                                                                             <  1:1481375/150‑1 (MQ=255)
                                                                           cGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACTCGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCg                                                                            >  1:1163742/1‑150 (MQ=255)
                                                                               ttGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACa                                                                        <  2:1719681/150‑1 (MQ=255)
                                                                               ttGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACa                                                                        >  2:736974/1‑150 (MQ=255)
                                                                                tGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                       >  1:2701543/1‑150 (MQ=255)
                                                                                tGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAg                                                                       >  1:46156/1‑150 (MQ=255)
                                                                                 gCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGc                                                                      >  1:614434/1‑150 (MQ=255)
                                                                                    gTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGACCTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTg                                                                   >  1:2385151/1‑150 (MQ=255)
                                                                                       cACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTgccg                                                                <  1:653174/150‑1 (MQ=255)
                                                                                          tATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGccgccg                                                             >  1:2877769/1‑150 (MQ=255)
                                                                                             ggCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTg                                                          >  1:1946407/1‑150 (MQ=255)
                                                                                                tATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAg                                                       >  2:1237993/1‑150 (MQ=255)
                                                                                                  tGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTg                                                     >  1:4316107/1‑150 (MQ=255)
                                                                                                      gATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTg                                                 >  1:2310354/1‑150 (MQ=255)
                                                                                                      gATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTg                                                 >  1:2673173/1‑150 (MQ=255)
                                                                                                       aTAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGACCAGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGa                                                >  2:1486538/1‑150 (MQ=255)
                                                                                                       aTAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGa                                                <  1:1103633/150‑1 (MQ=255)
                                                                                                         aaaGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACg                                              >  1:2940887/1‑150 (MQ=255)
                                                                                                         aaaGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACg                                              >  1:1040976/1‑150 (MQ=255)
                                                                                                         aaaGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACg                                              >  1:2607518/1‑150 (MQ=255)
                                                                                                         aaaGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACg                                              >  2:1630873/1‑150 (MQ=255)
                                                                                                          aaGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGt                                             >  1:1769800/1‑150 (MQ=255)
                                                                                                           aGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGACCTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTa                                            >  2:655324/1‑150 (MQ=255)
                                                                                                           aGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTa                                            >  2:2321017/1‑150 (MQ=255)
                                                                                                           aGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTa                                            >  2:2600589/1‑150 (MQ=255)
                                                                                                           aGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTa                                            >  2:3359944/1‑150 (MQ=255)
                                                                                                            ggCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCNCGTCACGCCNGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAt                                           >  1:3376428/1‑150 (MQ=255)
                                                                                                            ggCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAt                                           >  2:2579452/1‑150 (MQ=255)
                                                                                                                   tGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCg                                    <  1:3976388/150‑1 (MQ=255)
                                                                                                                    gACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCgg                                   >  1:809915/1‑150 (MQ=255)
                                                                                                                      ccGGCGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGtt                                 <  1:4395696/150‑1 (MQ=255)
                                                                                                                      ccGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGtt                                 <  1:1451731/150‑1 (MQ=255)
                                                                                                                      ccGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGtt                                 <  1:3198465/150‑1 (MQ=255)
                                                                                                                       cGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTa                                >  1:983816/1‑150 (MQ=255)
                                                                                                                        ggCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTAc                               >  2:4208671/1‑150 (MQ=255)
                                                                                                                          cGAATGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCt                             >  2:2106708/1‑150 (MQ=255)
                                                                                                                             aTGCCAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGc                          >  2:1684098/1‑150 (MQ=255)
                                                                                                                             aTGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTTGc                          <  2:2156331/150‑1 (MQ=255)
                                                                                                                              tGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCg                         <  2:3904389/150‑1 (MQ=255)
                                                                                                                               gCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCgg                        <  2:808602/150‑1 (MQ=255)
                                                                                                                                ccAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCggg                       >  1:1900863/1‑150 (MQ=255)
                                                                                                                                 cAGACGGTGGAGAACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGa                      >  1:696296/1‑150 (MQ=255)
                                                                                                                                 cAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGa                      >  1:2403609/1‑150 (MQ=255)
                                                                                                                                         ggAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCt              >  2:3585173/1‑150 (MQ=255)
                                                                                                                                            gaACCCGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTggg           >  2:2481953/1‑150 (MQ=255)
                                                                                                                                              aCCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCg         <  1:912492/150‑1 (MQ=255)
                                                                                                                                               cccGGATACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCgg        <  2:2999837/150‑1 (MQ=255)
                                                                                                                                                   gagaCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAcg    <  1:3288464/150‑1 (MQ=255)
                                                                                                                                                   gagaCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAcg    <  1:2733557/150‑1 (MQ=255)
                                                                                                                                                    agaCGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGAcgc   >  2:1427380/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCAGCATGAGACGAAACACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCC  >  CP009273/520862‑521160

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: