Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 334,871 C→T G51D (GGT→GAT)  yahC ← uncharacterized protein YahC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913334,8710CT91.7% 24.2 / ‑5.9 12G51D (GGT→GAT) yahCuncharacterized protein YahC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (7/4);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGCGTCATGCCGAGAAAGCCACCGGGAATAAAGGACAGTAACCAGAAACGCCCTTGCCAGACCACTACCGCCGAGAAGGGTACGGTAGTAATTACGCTCGCCCAGAACGCGGACAGTCCACTTACAGAGGCCAGCCAACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCCGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCG  >  NC_000913/334652‑335112
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
aaCAGCGTCATGCGGAGAAAGCCACCGGGAAAAAAGGACAGTAACCAGAAACGCCCTTGCCAGACCACGACCGCCGAGAAGGGTACGGTAGTAATTACGCTCGCCCAGAACGCGGACAGTCCACTTACAGAGGCCAGCCAACCAGCAGGCAGAGCGGCGAAAAACGCCCACACCAACCGGGAAGAATTATCCGCCAGACGTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGc                                                                                                                                                                                                                    <  1:219214/251‑1 (MQ=255)
                                           aCCAGAAACGCCCTTGCCAGACCACTACCGCCGAGAAGGGTACGGTAGTAATTACGCTCGCCCAGAACGCGGCCAGTCCACTTACAGAGGCCAGCCAACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATACTTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGc                                                                                                                                                                         >  1:44369/1‑251 (MQ=255)
                                                                                          gtaATTACGCTCGCCCAGAACGCGGACAGTCCACTTACAGAGGCCAGCCAACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGAc                                                                                                                          >  1:7953/1‑251 (MQ=255)
                                                                                                                cgGACAGTCCACTTACAGAGGCCAGCCAACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGac                                                                                                                    >  1:59815/1‑235 (MQ=255)
                                                                                                                                    gccagccaACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATGCCCATCCTTGAGAGAGCCCGa                                                                                                                                                                     <  1:9082/166‑1 (MQ=255)
                                                                                                                                      cagccaACCAGCAGTCAGCGCGGCGAAACACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATCAAAACTGCAAAAGGCGACGAAGCCAATGGTTNCTANNNNNTcnn                                                                                                                                                                                       >  2:157796/1‑144 (MQ=255)
                                                                                                                                      cagccaACCACAAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGAAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTgn                                                                                                                                                 <  1:157796/184‑2 (MQ=255)
                                                                                                                                                          cgacgAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACAGGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGGTGCCGTGAGTCCGTTCATCt                                                                                          <  1:90400/216‑1 (MQ=255)
                                                                                                                                                            gcgAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATAGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGACGCCAATGATTCATACTAATTCACATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACACATACCGACGGTGACACCCGTTGcc                                                                                                          >  1:35221/1‑201 (MQ=255)
                                                                                                                                                              naaaaaCGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTATCGCCGCCTGCGGCATaa                                                                                                                                                                                                                                  >  1:156216/2‑79 (MQ=255)
                                                                                                                                                              gAAAAACGCCCACCCCACCCCGGAATAACTATCCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATaa                                                                                                                                                                                                                                  <  2:156216/79‑1 (MQ=255)
                                                                                                                                                                                                                  cccGATTTATCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCCGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCg  >  1:92964/1‑251 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
AACAGCGTCATGCCGAGAAAGCCACCGGGAATAAAGGACAGTAACCAGAAACGCCCTTGCCAGACCACTACCGCCGAGAAGGGTACGGTAGTAATTACGCTCGCCCAGAACGCGGACAGTCCACTTACAGAGGCCAGCCAACCAGCAGTCAGCGCGGCGAAAAACGCCCACACCATCCCGGAATAATTAACCGCCAGACTTCTGATAAACCCCGATTTACCGCCGCCTGCGGCATAAAAACTGCAAAAGGCGACGAAGCCAATGGTTCCTAATAATTCCCATCCTTGAGAGAGCCCGACATGTGAAGACACTAACTGCCAAAGCCCTGCACAAATACCGACGGTGACACCCGTTGCCGTGAGTCCGTTCATCTTGTTCCCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCG  >  NC_000913/334652‑335112

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