Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,497,385 4 bp→CAA coding (530‑533/1656 nt) opgD → glucans biosynthesis protein D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,497,3850TC80.0% 20.1 / 1.7 15F177S (TTT→TCT) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (5/7);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01
*NC_0009131,497,3860T.85.7% 34.8 / 3.9 15coding (531/1656 nt)opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (5/7);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01
*NC_0009131,497,3870CA86.7% 32.4 / ‑2.2 15L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (6/7);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01
*NC_0009131,497,3880TA86.7% 30.3 / ‑2.8 15L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (6/7);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAAGTCAGGTG  >  NC_000913/1497140‑1497620
                                                                                                                                                                                                                                                     |||                                                                                                                                                                                                                                         
cGTTGAGATTCGTCAACTGGACGCTTAGTTCTTCCATATGGGAATGGGATTTCGGCGCCGCGTTCGTATGTTTTCTGTATATCCAGCAACACATCTGGCGCGTGAAAGTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGAGACAAAACAATTAGAATGGCAAAGCGATCGCGGCTTTGCGGGTTTTCGCGGGTTTAAAGCCCCCGAGCTGGCGCGCCGTGATGTAGTTTCTTTTCTCg                                                                                                                                                                                                                                        <  2:54933/251‑1 (MQ=255)
                  ggACTCTCAGTTCTTCCATAGGGGAATGGGATTCCGGCGCCGCTTTCGTATGTTTTCTGTAGATCCAGCAACACAGCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCcgc                                                                                                                                                                                                                       <  2:172798/250‑1 (MQ=255)
                                       gggtatgggaTTACGTCGCCGCGTTCGTATGCTTTCTGTAGCTCCAGTAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCTGGTTTGCGCTTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGGGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATa                                                                                                                                                                                                <  1:138702/247‑1 (MQ=255)
                                              ggattccgactccgtgtgtgTATGGTTTTTGTAGTTCCAGTAACACATCTGGCGTGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACGACGATGCAGGTGTTGATACAAAACAATTAGAAGGTAAAAGCGATCTCGTCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCGGTGATGTAGTTTCATCAACGGCGCGAGTTATTTCCGCGCTGTTGATGATACATATTAATACGGTTTg                                                                                                                                                                                         <  2:384475/232‑1 (MQ=255)
                                                                   aTGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATc                                                                                                                                                                    <  1:297399/251‑1 (MQ=255)
                                                                                 tCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTCGTATCATCAACGGCGCGGGTTATTTCCGCGCCGTTGATGATACATATCAATCAGGTTTGTCAGCCCGCGCCCTGTCGCTGGACACTTCCACCg                                                                                                                                                       >  2:9684/1‑250 (MQ=255)
                                                                                    aGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAg                                                                                                                                                   <  1:367928/251‑1 (MQ=255)
                                                                                                                           gttgttCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGCCAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTg                                                                                                             >  1:76995/1‑250 (MQ=255)
                                                                                                                                        aaCGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTCAAGCCCCCGAACTGGCGCGCAGTGATGTAGTATCATCAACGGCGAGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGATTGATACGGTAAAACCg                                                                                               >  1:19632/1‑251 (MQ=255)
                                                                                                                                              gCAGGTGTTGATACAAAACAATTAGAAGGGAAAAGCCATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCCCGCCGTGATGTAGTATCATCAACGGCGCGATTTATTTCCGCTCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGc                                                                                          <  1:343998/250‑1 (MQ=255)
                                                                                                                                                 ggTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACt                                                                                      >  1:353358/1‑251 (MQ=255)
                                                                                                                                                    gtTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGCGTTTAAAGCCCCCGCACCGGCGAGCAGTGATGTCGTATAATCAACCGCGCGAGATATTTCCACGCCGTTGATGATACATATCACTACAGTTTGCCGGCCCGCGGCCAGGAGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGACTTCTGACTTGATACGGCAAAACCgcaggcccatacc                                                                                   >  1:348315/1‑239 (MQ=255)
                                                                                                                                                             aaaCAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTGCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTGTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGtt                                                                          <  1:135992/251‑1 (MQ=255)
                                                                                                                                                                   ttAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGATGATGATACATATCAATACGGGTTGTCGGCCCGCGGCCTGGAGATCGACACTTACACCGACAGTAAAGTAGAGTACCCCGACTTTAACGCCGTCTGGTTTGATACGGTAAAACCGGGGGCAACTAACTTTACCGTTTATg                                                                      >  1:216699/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cGTGATGTAGTATCATCAACGGCGCTAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCGGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCGCGATAGCGCTAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAAGTCAGGTg  <  1:224597/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |||                                                                                                                                                                                                                                         
CGTTGAGAACCGTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAAGTCAGGTG  >  NC_000913/1497140‑1497620

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