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,3850TC100.0% 49.9 / NA 15F177S (TTT→TCT) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/9);  total (6/9)
*NC_0009131,497,3860T.93.3% 54.5 / ‑0.7 15coding (531/1656 nt)opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (5/9);  minor base A (1/0);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01
*NC_0009131,497,3870CA93.3% 49.2 / NA 15L178I (CTC→ATC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (5/9);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,497,3880TA93.3% 42.1 / ‑2.4 15coding (533/1656 nt)opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (5/9);  minor base . (1/0);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAA  >  NC_000913/1497165‑1497589
                                                                                                                                                                                                                            |||                                                                                                                                                                                                          
cAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTtgat                                                                                                                                                                               <  2:12386/251‑1 (MQ=255)
cAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTtg                                                                                                                                                                                 <  2:180701/249‑1 (MQ=255)
 agtgcttcCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTtgatg                                                                                                                                                                              <  1:28154/247‑1 (MQ=255)
                                      tCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGc                                                                                                                                         >  1:58554/1‑251 (MQ=255)
                                               ttCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGaca                                                                                                                                 <  1:285121/250‑1 (MQ=255)
                                                                                         tCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTACAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGtt                                                                                                                                                                                  >  2:178834/1‑159 (MQ=255)
                                                                                         tCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGtt                                                                                                                                                                                  <  1:178834/159‑1 (MQ=255)
                                                                                                  gttgttCAACTACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCAGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTCCACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTTGTTTgc                                                                             >  1:217850/1‑250 (MQ=255)
                                                                                                    tgttCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGAt                                                                            <  1:74827/250‑1 (MQ=255)
                                                                                                              cAACGATTCAGGTTTGGATACAAAACAATTAGAAGGGCAAATCGATCTCGGCTTTGCCGGTTTTCGCTTGTTTAAAGCCCCCGAACTGGCTCGCCGGGATGTAGTATCATCAACGGCGCAAGTTATTTCCGCGCCGTTAATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACTGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGGTTGATACGGTAAAAcc                                                                 <  2:252462/251‑1 (MQ=255)
                                                                                                                           gtTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGACGGTATTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGCGTTATTTCCGCGCCGTTGATGCTACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGAAAAGAAGAGCTCCCCGACTTTACCGCCTTCTGGTTTTAAACGGAACAACCGGGGGCAACtccc                                                    >  2:327428/1‑248 (MQ=255)
                                                                                                                             tGATACAACCCAATTAGCAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCACCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCtt                                                  >  2:147553/1‑251 (MQ=255)
                                                                                                                             tGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCt                                                   <  1:271709/250‑1 (MQ=255)
                                                                                                                                                         gATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGcc                      >  1:48032/1‑251 (MQ=255)
                                                                                                                                                                              cgcgTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATCAACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATaa  <  2:11458/250‑1 (MQ=255)
                                                                                                                                                                                                                            |||                                                                                                                                                                                                          
CAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAA  >  NC_000913/1497165‑1497589

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