Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,449,251 A→C 14.6% S18R (AGC→CGC)  PP_5556 → bacteriophage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2510AC14.6% 358.9 / 7.0 182S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (80/64);  new base C (9/23);  total (89/93)
Fisher's exact test for biased strand distribution p-value = 6.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCT  >  NC_002947/3449102‑3449399
                                                                                                                                                     |                                                                                                                                                    
aaGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                                      <  1:1728986/150‑1 (MQ=255)
  gggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                                                    >  1:3453745/1‑150 (MQ=255)
  gggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgc                                                                                                                                                    >  2:1759980/1‑150 (MQ=255)
   ggCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCcgcg                                                                                                                                                   >  1:2797136/1‑150 (MQ=255)
         aTCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcg                                                                                                                                             <  1:2122567/150‑1 (MQ=255)
          tCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                                                            <  2:3811090/150‑1 (MQ=255)
           ccTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggc                                                                                                                                           <  1:3037355/150‑1 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGgcggcg                                                                                                                                          >  2:2196911/1‑150 (MQ=255)
            cTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCGCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcggcc                                                                                                                                          <  1:3830881/150‑2 (MQ=255)
             tGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGgcgc                                                                                                                                         >  2:2229055/1‑150 (MQ=255)
                 gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCgg                                                                                                                                     >  1:3000256/1‑150 (MQ=255)
                 gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCgg                                                                                                                                     >  2:3886040/1‑150 (MQ=255)
                  aTCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg                                                                                                                                    <  2:3692060/150‑1 (MQ=255)
                  aTCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCTGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCggg                                                                                                                                    <  2:1415675/150‑1 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGCTGCTGCAGAGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCGGGCAg                                                                                                                                 >  2:3801493/1‑150 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGGGccg                                                                                                                                 >  1:3197924/1‑148 (MQ=255)
                     aCCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGccg                                                                                                                                 >  1:2785024/1‑148 (MQ=255)
                      ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                >  2:2964481/1‑150 (MQ=255)
                      ccAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGc                                                                                                                                >  2:1056347/1‑150 (MQ=255)
                          aCGTATCTTTGGTGGTAGCATGGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGTGGTGCTGCTGTGCGCGGTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTgg                                                                                                                            <  1:3242088/150‑1 (MQ=255)
                          aCGTATCATTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGCGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCGTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTgg                                                                                                                            <  1:1045279/150‑1 (MQ=255)
                              ctCAAGGGTGGTAGCAACGGTCTGGAGGATCGCACGGCCCTCTACAAGCGCGCGCTCGAGGTGCTGCAGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGGGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGgtgt                                                                                                                        <  2:3678875/149‑1 (MQ=255)
                                cattggtggTAGCAAGGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGTGGTGCTGCCGTGAGCGCTTGGGCCCCCCGCCTGGTCGCGGTGGCGCTGGTGCGGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTgg                                                                                                                      <  1:2406279/147‑1 (MQ=255)
                                cAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTgg                                                                                                                      <  2:734844/150‑1 (MQ=255)
                                 aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGc                                                                                                                     >  2:4421779/1‑150 (MQ=255)
                                     gtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                                 >  2:2274432/1‑150 (MQ=255)
                                       ggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCa                                                                                                               >  2:1993253/1‑150 (MQ=255)
                                           gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGc                                                                                                           >  2:4091429/1‑150 (MQ=255)
                                            cAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGGGGCCGGCCAACGcc                                                                                                          >  1:2786832/1‑150 (MQ=255)
                                             aaCGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCt                                                                                                         >  2:4142775/1‑150 (MQ=255)
                                                 gTCTGGTGGATCGCACGGCCCTCCCCACGCGTGCGCTCGTGGTGCTGCAGGGCGCGCTTGGCCCACCCGCGTGGTCGCCGTGGGGGTGGGGGGGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATaa                                                                                                     <  1:3429825/150‑1 (MQ=255)
                                                  tcgggaggaTCGCAAGGCCCTCTACACGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAg                                                                                                    <  2:4309456/147‑1 (MQ=255)
                                                    gggaggaGCGCCCGGCCCTCTACCCGCGCGCGCGGGGGGTGCTGCTGGGCGGGCTGGGGCCCCCCGCCTGGGCGGCGGGGCGGGGCTGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgc                                                                                                  <  1:1273837/149‑1 (MQ=255)
                                                      gaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                                                >  2:728168/1‑150 (MQ=255)
                                                       aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCa                                                                                               >  2:4343984/1‑150 (MQ=255)
                                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCAAAAAGGCGCAGCTCGCTg                                                                                       >  1:1580294/1‑150 (MQ=255)
                                                                 ggCCCCCTACCCGCGCGCGCTCGAGGGGCTGCTGTGAGCGGTTGGGCCACCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGAt                                                                                     <  1:2606422/150‑1 (MQ=255)
                                                                    ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCAAAAAGGCGCAGCGCGCGGATCAg                                                                                  >  1:1884817/1‑150 (MQ=255)
                                                                     ctctACACGCGAGCGCGCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg                                                                                 <  2:1922951/150‑1 (MQ=255)
                                                                      tctACAAGCGAGCGCTCGAGTTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                <  2:815612/150‑1 (MQ=255)
                                                                      tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc                                                                                >  2:3447888/1‑150 (MQ=255)
                                                                        tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCCCCCCGCCAGCCAGGGTGTGGCAGCCCAACCCCGATAAGGCGCAGCTCGCTGAGCAGGccc                                                                              >  1:3235563/1‑149 (MQ=255)
                                                                             gcgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGAc                                                                         >  2:1568254/1‑150 (MQ=255)
                                                                             gcgagcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACCCCGATAATGCGCAGCACGCGGATCAGGCCGATGAc                                                                         >  1:1223504/1‑150 (MQ=255)
                                                                               gagcgGTCGAGGTGCTGCTGGGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGGGCTGCTGCGGGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACtt                                                                       <  2:2378113/150‑1 (MQ=255)
                                                                               gagcgCTCGTGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACtt                                                                       <  1:4510761/150‑1 (MQ=255)
                                                                                   gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCGGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAAGTCATCGGc                                                                   <  2:3002591/150‑1 (MQ=255)
                                                                                     tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                 >  2:3156463/1‑150 (MQ=255)
                                                                                     tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCaa                                                                 >  2:1320658/1‑150 (MQ=255)
                                                                                       gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGATGACGTCGGCAAGc                                                               >  1:2606073/1‑150 (MQ=255)
                                                                                        tggTGCGGCTGTGTGCGCTTGGGCCACCCGCCTGGTCGCGGTGGGGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCa                                                              <  1:4311873/149‑1 (MQ=255)
                                                                                        cggTGCTGCTGGGCGCGCTGGGGCCACCCGCCTGGTCGCGGGGGCGCTGCTGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCa                                                              <  1:2152397/149‑1 (MQ=255)
                                                                                             ctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGGGTGGCCGGCCAACGCCGATAAGGCGCAGCGCGCAGATCAGGCCGATGACTTCGGCAAGCAGGTgg                                                         >  1:131357/1‑150 (MQ=255)
                                                                                              tgcTGTGCGCGCTTGGGCCCCCCGCCTGGTCGCCGTGGCGCGGCTGCTGGCCGCCTGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                        <  1:441073/150‑1 (MQ=255)
                                                                                              tgcCGTGCGCGGTTGGGCCCCCCGGGGGGGCGCCGTGGCGGGGGGGGGGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                        <  1:1970934/150‑1 (MQ=255)
                                                                                              tgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGc                                                        >  2:4193841/1‑150 (MQ=255)
                                                                                               gcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATACGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCg                                                       >  2:2160898/1‑150 (MQ=255)
                                                                                                 tgTGCGCGGTGGGGCCCCCCGCCTGGTCGCGGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGga                                                     <  1:1000297/149‑1 (MQ=255)
                                                                                                   tGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGgaga                                                   <  2:2474798/150‑1 (MQ=255)
                                                                                                      ggggttGGGCCACCCGCGGGGTCGGGGTGGGGCTGGTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAgg                                                <  1:4229479/146‑1 (MQ=255)
                                                                                                      gcgcTGGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCGGGCCAACGCGGATAAGGCGCAGCGGGCTGATCAGGCCGGTGACTTCGGCCAGCGGTTGGCGGAGAAgg                                                >  1:2314718/1‑150 (MQ=255)
                                                                                                       cgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGa                                               <  2:3289513/150‑1 (MQ=255)
                                                                                                       cgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGa                                               >  2:3054732/1‑150 (MQ=255)
                                                                                                         cTTGGGCCACCCGCGTGGTCGCGGGGGGGCGGCGGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                             <  1:3024689/150‑1 (MQ=255)
                                                                                                         cTTGGGCCACCCGCGTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                             <  1:643613/150‑1 (MQ=255)
                                                                                                         cTTGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGGGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATc                                             <  1:28394/150‑1 (MQ=255)
                                                                                                           tGGGCCACCCGCCTGGTCGCCGTGGCGCGGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCCATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgc                                           <  1:3348420/150‑1 (MQ=255)
                                                                                                                   cggcggggTCGCGGGGGCGCTGCTGCGGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                   <  1:4265775/144‑1 (MQ=255)
                                                                                                                   ccGCGTGGTCGGCGTGGGGGTGCTGCTGGGGGGCAGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCGTTTAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCATGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                   <  1:3355918/150‑1 (MQ=255)
                                                                                                                   ccGCCTGGTCGCCGTGGCGGTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCTGGCCGTTGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgg                                   <  2:364967/150‑1 (MQ=255)
                                                                                                                     gcgtGGTCGCGGTGGCGGTGCTGGTGGCCGCCAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTc                                 <  1:3978575/147‑1 (MQ=255)
                                                                                                                      ccTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCTAAAAGGCGCAGCTCGCGGATCAGGCCGAGGACTTCGGCAAGCACTTGGCGGAGAAGGATCGCGCGCACGGTCt                                >  1:1402668/1‑150 (MQ=255)
                                                                                                                       gtGGTCGGGGTGGCGCTGCTGCTGGCCGCCTGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTg                               <  1:3597160/149‑1 (MQ=255)
                                                                                                                         ggTCGGGGTGGGGCTGCTGCTGGCCGCCCGGGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGa                             <  1:3627103/150‑1 (MQ=255)
                                                                                                                          gTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAg                            <  1:330543/150‑1 (MQ=255)
                                                                                                                             gCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGCAAGCACTTGGCGGAGAAGGATAGCGCGCACGGTCTGGAGCgg                         >  1:3481674/1‑150 (MQ=255)
                                                                                                                             gCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCCGGCCGCCTGGGTGTGGCAGGCCAACGCCGCTAAGGCGCAGCTCGCTGATCAGGCCCATGACTTCGGCAACCCGGTGGCGGAGAAGGATCGCGCGCCCgggcgggggcgg                         >  1:2696888/1‑150 (MQ=255)
                                                                                                                               cGTGGCGCTGCTGCTGGCCGCCAGCGCTTGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAAAGCCTATAATGCGCAGCTCGCTGATCAGGCCGATTAATTAGGCAATCAGTTTGCGGAGAAGTATAGCGCGCACGGTCTGGAGCGGga                       >  1:2934528/1‑150 (MQ=255)
                                                                                                                               cGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                       <  2:536794/150‑1 (MQ=255)
                                                                                                                               cGTGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga                       <  2:1172472/150‑1 (MQ=255)
                                                                                                                                  ggCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa                    <  2:1551232/150‑1 (MQ=255)
                                                                                                                                   gcgcTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAg                   >  2:1226981/1‑150 (MQ=255)
                                                                                                                                   gcgcGGCTGCGGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAg                   <  2:574756/150‑1 (MQ=255)
                                                                                                                                    cgcgggtgcggGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGACGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg                  <  1:3149729/144‑1 (MQ=255)
                                                                                                                                    cgcTGCGGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg                  <  2:1661665/150‑1 (MQ=255)
                                                                                                                                     gctgctgctGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGc                 >  2:3219127/1‑150 (MQ=255)
                                                                                                                                      ctgctgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCt                <  1:376615/150‑1 (MQ=255)
                                                                                                                                         gtgctgGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca             <  1:872024/149‑1 (MQ=255)
                                                                                                                                         ctgctgGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGGGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca             <  2:119119/150‑1 (MQ=255)
                                                                                                                                         ctgctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca             >  2:2331707/1‑150 (MQ=255)
                                                                                                                                         ctgcggGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgca             <  1:63249/150‑1 (MQ=255)
                                                                                                                                          ggcgggCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcag            <  1:457656/146‑1 (MQ=255)
                                                                                                                                           gctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCTGGAGAAGGCTgcagc           >  1:92344/1‑150 (MQ=255)
                                                                                                                                           gctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagc           >  2:2144780/1‑150 (MQ=255)
                                                                                                                                            ctgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagca          <  1:694580/150‑1 (MQ=255)
                                                                                                                                             tgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcag         >  1:1706239/1‑150 (MQ=255)
                                                                                                                                             tgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcag         >  2:3745532/1‑150 (MQ=255)
                                                                                                                                             tgGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTgcagcag         >  1:4124730/1‑150 (MQ=255)
                                                                                                                                                 cgccAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCagcagcagc     <  1:2519930/150‑1 (MQ=255)
                                                                                                                                                   ccAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAgcgc   <  1:305912/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGGCGACATCCTGGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCT  >  NC_002947/3449102‑3449399

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: