Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,452,283 T→G 13.6% A19A (GCT→GCG PP_5556 → bacteriophage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,452,2830TG13.6% 422.6 / 10.8 215A19A (GCT→GCGPP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (94/82);  new base G (14/22);  total (111/104)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

CGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  NC_002947/3452134‑3452432
                                                                                                                                                     |                                                                                                                                                     
 gACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTg                                                                                                                                                      <  2:5764301/150‑1 (MQ=255)
       cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGACCCCTCCAGGTGCTGCAGTTAGCGCTTGGGCCACCCGCCTGGTCGACGAGGCGCTGCAGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                                >  2:3852814/1‑150 (MQ=255)
       cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                                >  2:430905/1‑150 (MQ=255)
       cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                                >  1:2608060/1‑150 (MQ=255)
        tGGCGATCACGAAACGTATGAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACTAGCGAGCGCTCGAGTTGCTGCTGTGATCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcga                                                                                                                                               <  2:37281/150‑2 (MQ=255)
                aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAg                                                                                                                                       >  1:5085535/1‑150 (MQ=255)
                aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAg                                                                                                                                       >  2:5530336/1‑150 (MQ=255)
                 ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                      >  1:5064156/1‑150 (MQ=255)
                  cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGCGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGcc                                                                                                                                     <  1:2235804/150‑1 (MQ=255)
                   aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCt                                                                                                                                    >  2:5591337/1‑150 (MQ=255)
                     aCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                                  <  2:6832204/150‑1 (MQ=255)
                        tATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                                               >  1:3361181/1‑150 (MQ=255)
                           cAATGGGGGTCGCACCGGTGGGGAGGATCGCAAGGCCCTCTCCACGCGCGCGCTCTAGGGGCTGCAGTGAGCGCTGGGGCCACCCGCCGGGCCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTgg                                                                                                                            <  2:5355480/150‑1 (MQ=255)
                             agggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACCCGCGCGCGCTCGCGGTGCTGCAGTGCGCGCTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                          <  1:2211236/148‑1 (MQ=255)
                             aTGGTGGTAGCAACGGTCTGGAGGATCGCACGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCa                                                                                                                          <  1:6820571/150‑1 (MQ=255)
                              tggtggTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAg                                                                                                                         >  1:3665112/1‑150 (MQ=255)
                                 tggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                                                      >  1:2898611/1‑150 (MQ=255)
                                 tggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGcc                                                                                                                      >  2:3931658/1‑150 (MQ=255)
                                    gaGCAACGGTCTGGAGGATCGCACGGCCCTCTACAAGCGAGCGCTCTAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAc                                                                                                                   <  2:3345966/149‑1 (MQ=255)
                                        aaCGGTCTGGTGGAACGCACGGCCCTCTACACGCGCGCGCTCGCGGTGCTGCTGTGTGCGCTTGGGCCCCCCGCCTGGTCGCCGTGGCGCTGCTGGTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCt                                                                                                               <  1:618003/150‑1 (MQ=255)
                                         aCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCaaacccaa                                                                                                              >  1:2642157/1‑144 (MQ=255)
                                           ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCCGGCCAACGCCgata                                                                                                            >  1:4326636/1‑150 (MQ=255)
                                           ggTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggc                                                                                                                                                    >  1:1159915/1‑110 (MQ=255)
                                            gTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCCGGCCCACGccgagac                                                                                                           >  1:782752/1‑145 (MQ=255)
                                              cTGGAGTATCGCATGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAgg                                                                                                         <  2:4457246/150‑1 (MQ=255)
                                               tGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGAAAAGgc                                                                                                        >  1:3027791/1‑150 (MQ=255)
                                                 gaggaTCGCCAGGCCCTCTACACGCGCGCGCTGGTGGGGCTGCAGTGAGCGCTTGGGCCCCCCGCGGGGTCGCGGGGGCGGTGCTGCGGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                                      <  1:4673728/150‑1 (MQ=255)
                                                 gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                                      >  2:5748149/1‑150 (MQ=255)
                                                  aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCa                                                                                                     >  2:6671251/1‑150 (MQ=255)
                                                   ggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAg                                                                                                    >  1:6536315/1‑150 (MQ=255)
                                                   ggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAg                                                                                                    >  2:6275847/1‑150 (MQ=255)
                                                    gaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                                                   >  2:4468320/1‑150 (MQ=255)
                                                         cacggCCCTCTACCCGCGTGCGGTCGTGGTGCTGCGGTGAGGGCTTGGGCCCCCCGCGTGGTCGCGGTGGGGCTGCTGCGGGCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                              <  1:2426176/147‑1 (MQ=255)
                                                         cAAGGCCCTCTACACGCGAGCGCTGGGGGTGCTGCAGTGAGCGCGTGGGCCACCCGCGTGGGCGCGGGGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                              <  1:4318537/150‑1 (MQ=255)
                                                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                                              >  2:2154632/1‑150 (MQ=255)
                                                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCt                                                                                              >  1:849667/1‑150 (MQ=255)
                                                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCTCGCTg                                                                                             >  2:2329850/1‑150 (MQ=255)
                                                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTCTAAGGCGCAGCTCGCTg                                                                                             >  2:3487993/1‑150 (MQ=255)
                                                             gCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCGCGCTGATc                                                                                          >  2:736148/1‑150 (MQ=255)
                                                               ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCGCGCTGATCAg                                                                                        >  1:4990423/1‑150 (MQ=255)
                                                               ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGGGGCACGCCCACCCCGCTACGGCGCAGCGCGCTGATCAg                                                                                        >  1:5317783/1‑150 (MQ=255)
                                                               ccTCAACACGGGAGCGCTCGCGGTGCTGCCGGGAGCGCTTGGGCCACCCGCGTGGTCGCGGTGGGGCGGCTGCTGGCCGCCTGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAg                                                                                        <  2:878479/150‑1 (MQ=255)
                                                                ctctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg                                                                                       >  1:3871640/1‑150 (MQ=255)
                                                                   tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                                    >  2:2800806/1‑150 (MQ=255)
                                                                   tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCTATAAGGCGCAGCTCGCTGATCAGGCCg                                                                                    >  1:5405244/1‑150 (MQ=255)
                                                                    aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGa                                                                                   >  1:819105/1‑150 (MQ=255)
                                                                    aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGa                                                                                   >  1:5804458/1‑150 (MQ=255)
                                                                     cacgcgcgcgcTCGTGGTGCTGCCGTGAGCGCTTGGGCCACCCGCCTGGTCGCGGGGGCGCTGCTGCTGCCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAt                                                                                  <  1:5557379/147‑1 (MQ=255)
                                                                     cAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAt                                                                                  >  2:5914095/1‑150 (MQ=255)
                                                                      aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATg                                                                                 >  2:1160059/1‑150 (MQ=255)
                                                                       cgcgcgcgcTCGAGGTGCGGCTGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGGGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGa                                                                                <  2:624447/149‑1 (MQ=255)
                                                                         cgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACt                                                                              <  2:425493/150‑1 (MQ=255)
                                                                            gcgcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCg                                                                           >  2:6495325/1‑150 (MQ=255)
                                                                              gctggtggTGCTGCTGTGTGGGGTTGGGCCACCGGCGTGGTCGCGGTGGCGGTGCTGGTGGCCGCCTGGGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGc                                                                         <  1:6118246/144‑1 (MQ=255)
                                                                               cTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCa                                                                        >  2:2595294/1‑150 (MQ=255)
                                                                                tCGCGGTGCTGCCGTGCGCGCTTGGGCCACCCGGCTGGTCGCGGTGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                       <  1:4851279/150‑1 (MQ=255)
                                                                                tCGAGGTGCTGCTGTGAGCGCTTGGGCCACCCGCGTGGTCGCGGTGGGGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                       <  1:1405275/150‑1 (MQ=255)
                                                                                tCGAGGTGCTGCAGTGCGCGCGTGGGCCACCCGCCTGGGCGCCGGGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                       <  2:5984917/150‑1 (MQ=255)
                                                                                 ggcggTGCTGCTGTGCGCGGGGGGGCCACCGGCCTGGTCGCGGGGGCGGTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAg                                                                      <  1:2934558/147‑1 (MQ=255)
                                                                                     gTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                                                  <  1:2758722/150‑1 (MQ=255)
                                                                                       gctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGGCCACGCCTAAAAGGCGCAGCTCGCTGATCAGGCCGATGACTTAGGCAAGCAGTTg                                                                >  1:1955059/1‑150 (MQ=255)
                                                                                            agggAGCGCTTGGGCCACCCGCATGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGTCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                           <  2:4609094/147‑1 (MQ=255)
                                                                                             gTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgag                                                          >  2:1533077/1‑150 (MQ=255)
                                                                                               gcgcgcTGGGGCCACCCGCCTGGTCGCCGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGAGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGaa                                                        <  2:4044288/148‑1 (MQ=255)
                                                                                                 gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAgg                                                      <  2:1047604/150‑1 (MQ=255)
                                                                                                    cTTGGGCCACCCGCCTGGTCGCGGGGGCGGTGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGACGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                                   <  1:989258/150‑1 (MQ=255)
                                                                                                    cTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                                   <  2:3652657/150‑1 (MQ=255)
                                                                                                     ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcg                                                  >  2:2651857/1‑150 (MQ=255)
                                                                                                     ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcg                                                  >  2:4706772/1‑150 (MQ=255)
                                                                                                         gCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGACGGATcgcgcg                                              >  2:4802531/1‑150 (MQ=255)
                                                                                                           cACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCa                                            >  2:6409463/1‑150 (MQ=255)
                                                                                                             cccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACg                                          <  2:2441411/150‑1 (MQ=255)
                                                                                                              ccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACCCCTCTAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGCGAAGGATCGCGCCCACgg                                         >  1:5197241/1‑150 (MQ=255)
                                                                                                              ccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                         >  1:4540038/1‑150 (MQ=255)
                                                                                                              ccGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                         >  2:4429166/1‑150 (MQ=255)
                                                                                                               cGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGGAGAGAAGGATCGCGCGCACggg                                        >  1:3597020/1‑149 (MQ=255)
                                                                                                                 ccTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCt                                      <  2:1670612/150‑1 (MQ=255)
                                                                                                                  cTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg                                     <  1:5672396/150‑1 (MQ=255)
                                                                                                                    ggTCGCCGTGGCGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                                   <  2:5064156/150‑1 (MQ=255)
                                                                                                                    ggTCGCCGGGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                                   <  2:3131484/150‑1 (MQ=255)
                                                                                                                      tCGCCGTGGCGCGGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGc                                 <  2:5302379/150‑1 (MQ=255)
                                                                                                                        gCCGTGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg                               <  2:3615397/150‑1 (MQ=255)
                                                                                                                         cgggggCGCGGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                              <  1:3232150/146‑1 (MQ=255)
                                                                                                                         ccGTGGCGCTGCTGGTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                              <  1:1667900/150‑1 (MQ=255)
                                                                                                                         ccGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCACCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCCATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                              >  1:2495558/1‑150 (MQ=255)
                                                                                                                         ccGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCggg                              <  2:5109440/150‑1 (MQ=255)
                                                                                                                          cGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                             <  2:6221799/150‑1 (MQ=255)
                                                                                                                             ggCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                          <  1:3052650/150‑1 (MQ=255)
                                                                                                                               cgcTGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg                        <  1:100568/150‑1 (MQ=255)
                                                                                                                                gctgctgctGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                       >  2:5220227/1‑150 (MQ=255)
                                                                                                                                  tgctgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTg                     >  2:5188144/1‑150 (MQ=255)
                                                                                                                                    ctgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca                   <  2:2015066/150‑1 (MQ=255)
                                                                                                                                     tgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag                  >  2:4300206/1‑150 (MQ=255)
                                                                                                                                     ggcgggCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag                  <  2:5735457/146‑1 (MQ=255)
                                                                                                                                        gcgccgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCCAGCAGTTGGCGGAGaa                                                        <  2:6157278/107‑1 (MQ=255)
                                                                                                                                              ccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc         >  1:5903793/1‑150 (MQ=255)
                                                                                                                                               aaGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCACGCAGTTGGCGGAGAAGGATCGCGCGCCCGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCt        >  2:5428457/2‑150 (MQ=255)
                                                                                                                                                   gcTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAt    >  2:2847814/1‑150 (MQ=255)
                                                                                                                                                     tGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCa  <  2:3052423/150‑1 (MQ=255)
                                                                                                                                                     tGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCa  >  1:3111898/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCA  >  NC_002947/3452134‑3452432

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