Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,728,186 C→T 22.2% C253Y (TGC→TAC)  opgH ← glucans biosynthesis glucosyltransferase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,728,1860CT22.2% 29.0 / 5.1 18C253Y (TGC→TAC) opgHglucans biosynthesis glucosyltransferase H
Reads supporting (aligned to +/- strand):  ref base C (5/9);  new base T (2/2);  total (7/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

TGCTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCTTGTATTTGTCACGCCCGGTGAGCAGCTCGAGGAAGCCCATCA  >  NC_002947/5728051‑5728300
                                                                                                                                       |                                                                                                                  
tGCTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGc                                                                                                                    <  2:173058/136‑1 (MQ=255)
tGCTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGc                                                                                                                    <  2:268779/136‑1 (MQ=255)
  cTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGTAg                                                                                                                  >  1:10192/1‑136 (MQ=255)
  cTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGTAg                                                                                                                  >  1:234604/1‑136 (MQ=255)
  cTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAg                                                                                                                  >  1:141519/1‑136 (MQ=255)
  cTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAg                                                                                                                  >  1:175690/1‑136 (MQ=255)
       ccACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCgg                                                                                                             >  2:177930/1‑136 (MQ=255)
           ggCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATc                                                                                                         <  1:79880/136‑1 (MQ=255)
               aTGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCa                                                                                                     <  2:206928/136‑1 (MQ=255)
                    ggTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGcgcg                                                                                                <  1:129241/136‑1 (MQ=255)
                      tgtCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGTAGATCGGCATCACCAGCGCGGt                                                                                              <  2:234604/136‑1 (MQ=255)
                      tgtCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGTAGATCGGCATCACCAGCGCGGt                                                                                              <  2:10192/136‑1 (MQ=255)
                      tgtCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGt                                                                                              <  2:175690/136‑1 (MQ=255)
                                                      aaCCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTg                                                              >  1:60149/1‑136 (MQ=255)
                                                                         tGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCtt                                           >  2:166245/1‑136 (MQ=255)
                                                                                 tgctaCTCTAAGGTTGGGGGCAGGCCGGCGAACTCCGGCGGCTGGTTTTCGTTGCAGATGGGCATCACCTGGGCGGTACGGGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGATATGTTGTATTTGt                                   <  1:35961/132‑1 (MQ=255)
                                                                                                        gCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCTTGTATTTGTCACGCCCGGTGAGCAGCTCGAgg            <  2:65397/136‑1 (MQ=255)
                                                                                                                  acCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCTTGTATTTGTCACGCCCTGTTATAATCTCGAGGAAGCCCTTCa  <  1:113587/136‑1 (MQ=255)
                                                                                                                  acCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCTTGTATTTGTCACGCCCGGTGAGCAGCTCGAGGAAGCCCATCa  <  2:141519/136‑1 (MQ=255)
                                                                                                                  acCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCTTGTATTTGTCACGCCCGGTGAGCAGCTCGAGGAAGCCCATCa  <  2:13618/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                  
TGCTCGGCCACGGCAATGTCGGTGTCGTTGGTGTCGCTGAGCACAAAGAAGTCGAACCGGTCGAGGTTGCCGCTGGCGGCCACCGACTCGAACGTTGCGCGCAGGCCGGCGAACACCCGCGGCACGTCTTCGTTGCAGATCGGCATCACCAGCGCGGTACGCGCCTCGGGCGCGATCGGCTCGTTGCCTGCGCTGCTACCGGAAATCTTGTATTTGTCACGCCCGGTGAGCAGCTCGAGGAAGCCCATCA  >  NC_002947/5728051‑5728300

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