Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 204,796 C→T P102S (CCG→TCG)  rnhB → RNase HII

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913204,7960CT91.7% 22.6 / ‑3.8 12P102S (CCG→TCG) rnhBRNase HII
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/3);  total (8/4)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAGCGAAAAACGCCGTCTGGCGCTCTATGAAGAGATCAAAGAGAAAGCGTTGAGCTGGAGTCTGGGCCGCGCGGAACCCCACGAAATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTGGCTGAACACGGCGCGAC  >  NC_000913/204566‑205031
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
gcgcGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAGCGAACAACGCCGTCTGGCGATCTATGAAGAGATCAAAGAGAAAGCGTTGAGATGGAGTCTGGGCCGCGCGGAGCCCCACGAAATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGATGCATATTGCGTCGGAATATGTGTtgcttgat                                                                                                                                                                                                                         >  1:120524/1‑251 (MQ=255)
                 ncGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCACAAAGCTGAGCGAAAAACGCCGTCTGGCGCTCTATGACGAGATCAAAGAGAAAGCGTTGAGCTGGAGTCTGGGCCGCGCGGACCCCCACGAAATCGACGAGCTGAACATTCTTCCTGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGa                                                                                                                                                                                                         >  1:28963/2‑250 (MQ=255)
                                                                            cGAAAAACGCCGTCTGGCGCTCTATGAAGAGATCAAAGAGAAAGCGTTGAGCTGGAGTCTGGGCCGCGCGGAACCCCACGAAATCGAAGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTg                                                                                                                                             >  1:49571/1‑251 (MQ=255)
                                                                                 aaCGCCGTCTGGCGCTCTATGAAGAGATCAAAGAGAAAGCGTTGAGCTGGAGTCTGGGCCGCGCGGAAACCCACGAAATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATCTTGCGTCGGAATATGTGtt                                                                                                                                                                                                                                >  1:173651/1‑163 (MQ=255)
                                                                                                                                         ggCCGCGCGGAACCCCACGAAATCGACGAGCTGAACATTCTTCAGGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATcn                                                                                                                                 <  1:33817/202‑2 (MQ=255)
                                                                                                                                                     ccccACGAAATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGa                                                                                                                                                     >  1:81993/1‑170 (MQ=255)
                                                                                                                                                       ccACGAAATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCcg                                                                                                                                                              >  1:71746/1‑159 (MQ=255)
                                                                                                                                                         aCGAANNCGNNNAGNNGAACATTCTTCATGCGACCATGCNGGCGATGCAGCGTGCCGTCGCTGGGCTGCNTNTNGCGTNNGAANATGTGTTGATNGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGNANGGCGATAGCCGCGTACCGGAAATCAGTGCCGNGTCTATCCTGGCGAAAGAGACGCGNGACGNNg                                                                                                         >  2:123043/1‑210 (MQ=255)
                                                                                                                                                                                        gACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATTGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTGTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCttttn                                 <  1:20210/251‑2 (MQ=255)
                                                                                                                                                                                         aCCATGCGGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATGACCGATGCCTGCGAGGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGGGCCGCGTCTATCCTGGCGAAAGTTACGCGTGACGCCGAAATGGCGGCGCTGGAAATTGTTTTCCCGCAATAGGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCa                                <  1:156025/251‑1 (MQ=255)
                                                                                                                                                                                                        cAGCGTGCCGTCGCTGGGCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGAGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTGGc                 <  1:123043/251‑1 (MQ=255)
                                                                                                                                                                                                                       gggCTGCATATTGCGTCGGAATATGTGTTGATTGATGGTAACCGCTGCACGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGACGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAACGTGACGCGTGACGCCGACATGGCGGAGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCACCACAAAGGGTACCCAACCGCTTTTCATCTGGAAACAATGGCTGAACACGACGCGAc  >  1:113700/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
GCGCGGTCGTCACCGCTGCGGTGATCCTTGACCCGGCGCGCCCGATTGCCGGGCTGAATGATTCCAAAAAGCTGAGCGAAAAACGCCGTCTGGCGCTCTATGAAGAGATCAAAGAGAAAGCGTTGAGCTGGAGTCTGGGCCGCGCGGAACCCCACGAAATCGACGAGCTGAACATTCTTCATGCGACCATGCTGGCGATGCAGCGTGCCGTCGCTGGGCTGCATATTGCGCCGGAATATGTGTTGATTGATGGTAACCGCTGCCCGAAATTACCGATGCCTGCGATGGCTGTGGTGAAAGGCGATAGCCGCGTACCGGAAATCAGTGCCGCGTCTATCCTGGCGAAAGTGACGCGTGACGCCGAAATGGCGGCGCTGGATATTGTTTTCCCGCAATATGGTTTTGCCCAACACAAAGGGTACCCAACCGCTTTTCATCTGGAAAAACTGGCTGAACACGGCGCGAC  >  NC_000913/204566‑205031

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