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,3850TC87.5% 36.2 / ‑4.3 16F177S (TTT→TCT) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (8/6);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,497,3860T.93.3% 51.6 / ‑0.7 16coding (531/1656 nt)opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (8/6);  minor base A (0/1);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,497,3870CA100.0% 41.4 / NA 16L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (9/7);  total (9/7)
*NC_0009131,497,3880TA100.0% 44.4 / NA 16L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/7);  total (9/7)

AATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTG  >  NC_000913/1497182‑1497604
                                                                                                                                                                                                           |||                                                                                                                                                                                                                         
aatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTGAAAGCCTCCGAACTGGCGTGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACgg                                                                                                                                                                             <  1:242787/251‑1 (MQ=255)
             cgccgcGTTCGTATTTTTGCTGGAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATAGCAATACGTTTTTTCGGCCcg                                                                                                                                                                 <  2:233755/250‑1 (MQ=255)
                                                           gTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCg                                                                                                                  >  1:92364/1‑251 (MQ=255)
                                                                      tttCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATGTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGCTGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTCCACCCACAGTAAAGAAGAGTTCCCCGAAtttccaccc                                                                                                       >  2:175429/1‑245 (MQ=255)
                                                                               gAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGttt                                                                                              >  1:54550/1‑251 (MQ=255)
                                                                                                          gtggATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTGGCCGGGTTTCGCGTGTTTAAAGCCCCCGAACTGTCGCGCCGTGATTAAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGGTGATGATACATATCAATACGGTTTGGCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACGTTACCGCCCTCTGGTGTGATACGGTAAAACCGGGTGCATCTAcc                                                                   <  2:225699/248‑1 (MQ=255)
                                                                                                             gagaCAAAACAATGTGAAGGGCAAAGGGTGCTCGGCTTTGGCGGTTTTCGGGTGTTTAAAGCCCCCGACCTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGTGGGTTATTTCCGCGCCGTGGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTGACACCGGCAGTAAAGAAGAGTTCCCCGACTTTACCGCTTTCTGGTTTGATACGGTAAAACCGGGGGCAACTTCCttt                                                                <  2:20119/248‑1 (MQ=255)
                                                                                                             gATACAACGCAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTACAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGCGTTATTTCCGCGCCGTTGATGATACATCTCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGCCAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCttt                                                                >  2:34887/1‑251 (MQ=255)
                                                                                                                     aCAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGttt                                                         >  1:217245/1‑250 (MQ=255)
                                                                                                                              aGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGcc                                                                                                       <  1:84057/195‑1 (MQ=255)
                                                                                                                              aGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGcc                                                                                                       >  2:84057/1‑195 (MQ=255)
                                                                                                                               gggaaaaGCCATCTGGGCTTTGCCGGTTTTGGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGATTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTc                                              <  1:37736/247‑1 (MQ=255)
                                                                                                                                     agcgagttcGGCTTGGCCGGTTTTCGCGTTTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTAT‑AATAAACGGCGCGAGTTATTGCCGCGCCGTTGATGATACATATCAATACGGTTTGGCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTGTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCGCGATAgc                                        <  2:197359/244‑1 (MQ=255)
                                                                                                                                                       ggTTTTCGCGTGTTTAAAGCCCCCGAACTTGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGCGCGACACTTACACCGACAGTACAGAAGAGTTCCTCGACTTTACCGCCTTCTGGTTTGATACGGCAAAACCGGGCGCAACTACCTTTACCGTTTATGCGTCGCTCGATCACGCCAGCATTACTGGCGcc                      >  1:166438/1‑251 (MQ=255)
                                                                                                                                                                 tgtTTAAAGCCCCCGAACTGGCGAGCCGTGATGTAGTATC‑CTCAACGGTGCGAGTTATTTCCGCGACGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACAATTACACCGCCAGTACAGAAGACTTCACCGACTTTAACGCCTTCTGGTTTGATAAGGTCAAACCGGGGGCAACTACCTTTACCGTTTTTGCGTTGCTAGATAGAGCACGCACTACTGGTGCCTATAAGTTc             >  1:21015/1‑250 (MQ=255)
                                                                                                                                                                           ccccGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACCCTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATtg  >  2:103088/1‑251 (MQ=255)
                                                                                                                                                                                                           |||                                                                                                                                                                                                                         
AATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTG  >  NC_000913/1497182‑1497604

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