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
                                                                                                                                                                                                           |||                                                                                                                                                                                                                         
aatgggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTGAAAGCCTCCGAACTGGCGTGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACgg                                                                                                                                                                             <  1:242787/251‑1 (MQ=255)
             cgccgcGTTCGTATTTTTGCTGGAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATAGCAATACGTTTTTTCGGCCcg                                                                                                                                                                 <  2:233755/250‑1 (MQ=255)
                                                           gTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCg                                                                                                                  >  1:92364/1‑251 (MQ=255)
                                                                      tttCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATGTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGCTGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTCCACCCACAGTAAAGAAGAGTTCCCCGAAtttccaccc                                                                                                       >  2:175429/1‑245 (MQ=255)
                                                                               gAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGttt                                                                                              >  1:54550/1‑251 (MQ=255)
                                                                                                          gtggATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTGGCCGGGTTTCGCGTGTTTAAAGCCCCCGAACTGTCGCGCCGTGATTAAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGGTGATGATACATATCAATACGGTTTGGCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACGTTACCGCCCTCTGGTGTGATACGGTAAAACCGGGTGCATCTAcc                                                                   <  2:225699/248‑1 (MQ=255)
                                                                                                             gagaCAAAACAATGTGAAGGGCAAAGGGTGCTCGGCTTTGGCGGTTTTCGGGTGTTTAAAGCCCCCGACCTGGCGCGCCGTGATGTAGTATCATCAACGGCGTGGGTTATTTCCGCGCCGTGGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTGACACCGGCAGTAAAGAAGAGTTCCCCGACTTTACCGCTTTCTGGTTTGATACGGTAAAACCGGGGGCAACTTCCttt                                                                <  2:20119/248‑1 (MQ=255)
                                                                                                             gATACAACGCAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTACAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGCGTTATTTCCGCGCCGTTGATGATACATCTCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGCCAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCttt                                                                >  2:34887/1‑251 (MQ=255)
                                                                                                                     aCAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGttt                                                         >  1:217245/1‑250 (MQ=255)
                                                                                                                              aGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGcc                                                                                                       <  1:84057/195‑1 (MQ=255)
                                                                                                                              aGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGcc                                                                                                       >  2:84057/1‑195 (MQ=255)
                                                                                                                               gggaaaaGCCATCTGGGCTTTGCCGGTTTTGGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGATTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTc                                              <  1:37736/247‑1 (MQ=255)
                                                                                                                                     agcgagttcGGCTTGGCCGGTTTTCGCGTTTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATAATAAACGGCGCGAGTTATTGCCGCGCCGTTGATGATACATATCAATACGGTTTGGCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTGTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCGCGATAgc                                        <  2:197359/244‑1 (MQ=255)
                                                                                                                                                       ggTTTTCGCGTGTTTAAAGCCCCCGAACTTGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGCGCGACACTTACACCGACAGTACAGAAGAGTTCCTCGACTTTACCGCCTTCTGGTTTGATACGGCAAAACCGGGCGCAACTACCTTTACCGTTTATGCGTCGCTCGATCACGCCAGCATTACTGGCGcc                      >  1:166438/1‑251 (MQ=255)
                                                                                                                                                                 tgtTTAAAGCCCCCGAACTGGCGAGCCGTGATGTAGTATCCTCAACGGTGCGAGTTATTTCCGCGACGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACAATTACACCGCCAGTACAGAAGACTTCACCGACTTTAACGCCTTCTGGTTTGATAAGGTCAAACCGGGGGCAACTACCTTTACCGTTTTTGCGTTGCTAGATAGAGCACGCACTACTGGTGCCTATAAGTTc             >  1:21015/1‑250 (MQ=255)
                                                                                                                                                                           ccccGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACCCTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATtg  >  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: