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% 90.1 / NA 30F177S (TTT→TCT) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/17);  total (13/17)
*NC_0009131,497,3860T.100.0% 126.5 / NA 30coding (531/1656 nt)opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (13/17);  total (13/17)
*NC_0009131,497,3870CA96.7% 87.3 / NA 30L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (12/17);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
*NC_0009131,497,3880TA100.0% 89.6 / NA 30L178N (CTC→AAC) opgDglucans biosynthesis protein D
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/17);  total (13/17)

GTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTT‑CGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAA  >  NC_000913/1497151‑1497612
                                                                                                                                                                                                                                           |||                                                                                                                                                                                                                                 
gTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTT‑CGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTAtt                                                                                                                                                                                                                    >  1:337433/1‑251 (MQ=255)
          cGCTCAGTTCTTCAATATGGGAATGGGATTCCGTCGCCGCGTT‑CGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGAGCGCCGTGCTGTAGTATCATC‑CACGGCGCGAGTTATTTCCGCGCCGt                                                                                                                                                                                                          >  2:490087/1‑251 (MQ=255)
           gCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTT‑CGTATGTTTTCTGTAGAACCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGtt                                                                                                                                                                                                         <  2:265888/251‑1 (MQ=255)
                     tcCATATGGGAATGGGATTCCGTCGCCGCGGT‑CGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATAc                                                                                                                                                                                                 <  1:313098/249‑1 (MQ=255)
                                  gggaTTCCGGCGCC‑CGTTGAGTATGTTTTCTGGAGAGCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGttt                                                                                                                                                                                  <  2:2542/251‑1 (MQ=255)
                                                    t‑CGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTgg                                                                                                                                                                 <  2:133502/250‑1 (MQ=255)
                                                      cGTATGTTTTCTGTAGAGCCAACAACACATCGGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTGGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTGCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTATTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGc                                                                                                                                                                <  2:55401/250‑1 (MQ=255)
                                                              ttCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGaca                                                                                                                                                        <  1:395436/250‑1 (MQ=255)
                                                                     gATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCAGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGAGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACAcc                                                                                                                                                >  1:469312/1‑251 (MQ=255)
                                                                      aTCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGCTGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCg                                                                                                                                               >  2:243647/1‑251 (MQ=255)
                                                                        ccAGCAACACATCGGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGTCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGAc                                                                                                                                             <  1:382206/251‑1 (MQ=255)
                                                                          aGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTACCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGAc                                                                                                                                             >  2:391014/1‑249 (MQ=255)
                                                                             aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACACCGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGCATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACCTATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGCCACTTACAACGACAGTa                                                                                                                                         >  1:414745/1‑250 (MQ=255)
                                                                               cacaTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTaa                                                                                                                                        >  1:337448/1‑249 (MQ=255)
                                                                                      ggCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGt                                                                                                                                <  1:372535/250‑1 (MQ=255)
                                                                                                             cGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGt                                                                                                        >  1:41249/1‑251 (MQ=255)
                                                                                                                              aaCGATGCGGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGACCTGGCGCGCCGTGAGGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCg                                                                                       <  2:262705/251‑1 (MQ=255)
                                                                                                                                    gCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGc                                                                                  >  1:259527/1‑250 (MQ=255)
                                                                                                                                          gtTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTAcc                                                                           <  1:492471/251‑1 (MQ=255)
                                                                                                                                                aCAAAACAATTAGAAGGAAAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTAc                                                                      <  1:303659/250‑1 (MQ=255)
                                                                                                                                                   aaaCAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGtt                                                                  >  1:16372/1‑251 (MQ=255)
                                                                                                                                                       aaTTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATg                                                              >  1:353254/1‑251 (MQ=255)
                                                                                                                                                                   aaaGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATa                                                  <  1:391014/251‑1 (MQ=255)
                                                                                                                                                                   aaaGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGAGGGCAACTACCTTTACCGTTTATGCGTTGCTCGAt                                                   >  2:396395/1‑250 (MQ=255)
                                                                                                                                                                                  ttGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACGGGCGCGCCGTGATGTAGTACCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTAc                                     <  2:258332/249‑1 (MQ=255)
                                                                                                                                                                                    gCCGGTTGTCGCGGGTTAAAAGCCCCCGAGCTGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTCCCGCGCCGTTGATGATACATACCAATACGGTTTGGCGGCCCGCGGCCTGGCGATCGACTCTTACATCGACAGTAAAGAAGAGGTCCCCGACTTTACCGCCGTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTgg                                  <  2:112156/250‑1 (MQ=255)
                                                                                                                                                                                           ttCGTGTGTTTAAAGCCCCCGAACGGGCGCGCCGTGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCtat                           <  1:80525/250‑1 (MQ=255)
                                                                                                                                                                                           ttCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGCTGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAACGACGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCCGGGGCCACTACCTTTACCGTTTCTGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCtat                           >  2:28998/1‑250 (MQ=255)
                                                                                                                                                                                                gtgtTTAAAGCCCCCAAACGGTCGCGCCGGGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTGGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTc                     <  2:41249/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCGCGCCGGGATGTAGTATCATC‑AACGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGaaaa  <  1:445939/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |||                                                                                                                                                                                                                                 
GTCAACTGGACGCTCAGTTCTTCCATATGGGAATGGGATTCCGTCGCCGCGTT‑CGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCCCGCGGCCTGGCGATCGACACTTACACCGACAGTAAAGAAGAGTTCCCCGACTTTACCGCCTTCTGGTTTGATACGGTAAAACCGGGGGCAACTACCTTTACCGTTTATGCGTTGCTCGATAGCGCCAGCATTACTGGTGCCTATAAGTTCACTATCCATTGTGAGAAAA  >  NC_000913/1497151‑1497612

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