Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,665,014 C→T A170V (GCA→GTA)  dmsD → redox enzyme maturation protein DmsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,665,0140CT87.5% 15.7 / NA 8A170V (GCA→GTA) dmsDredox enzyme maturation protein DmsD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (3/4);  minor base G (0/1);  total (3/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCT  >  NC_000913/1664891‑1665231
                                                                                                                           |                                                                                                                                                                                                                         
nttCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGTAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAAcc                                                                                                                              >  1:197963/2‑217 (MQ=255)
                  cTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTTGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGGAGATCACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAAc                                                                                                                               <  1:301034/198‑1 (MQ=255)
                                                  gaATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGTAGAACACCCTTTCTACCGTGCCCTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCAc                                                                                                           >  1:137694/1‑186 (MQ=255)
                                                   aaTGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGTAGAACACCCTTTCTACCg                                                                                                                                                                                                         <  1:241395/91‑1 (MQ=255)
                                                                       nngcnNNNNTTTNCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGTAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAAcc                                                                                                                              <  2:197963/144‑1 (MQ=255)
                                                                                   ccGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGTAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGt                                                                                                                                          <  1:249092/122‑1 (MQ=255)
                                                                                          cAACACGTTTTCTTGATGTTTTTATCGAAAAAGTAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGTTTACAGTGTCACGTTTTGTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct  <  1:49624/251‑1 (MQ=255)
                                                                                                                    naaaaaGTAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGTAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTa                                                                                >  1:60874/2‑147 (MQ=255)
                                                                                                                           |                                                                                                                                                                                                                         
GTTCACTGCTGTTGATGGCTGCGTGGTTGGCGGAGAATGGTCGCCAGACAGAATGCGAAGAATTGCTGGCATGGCACCTTTTTCCGTGGTCAACACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCT  >  NC_000913/1664891‑1665231

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