Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 486,927 G→A R131H (CGC→CAC)  mscK → potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913486,9270GA100.0% 13.3 / NA 6R131H (CGC→CAC) mscKpotassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (2/4);  total (2/4)

AAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCA  >  NC_000913/486697‑487134
                                                                                                                                                                                                                                      |                                                                                                                                                                                                               
aaataaaCAAAAAGATCTTTCTGCTCATGACAAACTGGTGCAGCAGGATCGGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTTAAGCGCCGGAAAAAAGGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCACCAGCTGGAAACTCGCGt                                                                                                                                                                                               <  1:159519/249‑1 (MQ=255)
                                                                                      ncATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCACCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGa                                                                                                        >  1:41888/2‑250 (MQ=255)
                                                                                                                                                               ggcgTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCACCAGCTTGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACATCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATATCGTTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATg                              <  1:88976/251‑1 (MQ=255)
                                                                                                                                                                       cAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAACATTCTGAGCACGCTGTCGTTGCACCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCg                        >  1:27545/1‑249 (MQ=255)
                                                                                                                                                                                           aaCGACAAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCACCAGCTGGAAACTCGCGGTGCCCAGGCGCTGGACGTTTTGCAAAACGCACAAAACGATCTGTCGTCTTATAACAGCCAGCTGTTTTCGTTACAGACTCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGTATCAGCTGCAACAAATTCGCAGTCGTCTGGATGGTACTGATGTCGGCGAGACAGCCTTCCGTCccagcta  <  1:168035/251‑3 (MQ=255)
                                                                                                                                                                                                    naaaCGCGCAAAATTCTGAGNACGCTGTCGNTGCACNNNCNTGAANCTCGCGtt                                                                                                                                                                                              <  2:20287/53‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                               
AAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGGATCTGACAGATACATTAGCCACCCTCGATAAAATCGATCGCATAAAAGAAGAGACAGTTCAGCTACGGCAAAAAGTCGCTGAAGCGCCGGAAAAAATGCGCCAGGCGACCGCGGCGTTAACAGCACTTAGCGATGTCGATAACGACGAAGAAACGCGCAAAATTCTGAGCACGCTGTCGTTGCGCCAGCTGGAAACTCGCGTTGCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCA  >  NC_000913/486697‑487134

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