Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 245,295 C→T V171M (GTG→ATG)  dpaA ← peptidoglycan meso‑diaminopimelic acid protein amidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913245,2950CT100.0% 16.6 / NA 8V171M (GTG→ATG) dpaApeptidoglycan meso‑diaminopimelic acid protein amidase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (3/5);  total (3/5)

CGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCACTTGCACGCTCGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGCCGATGGAAACACAATCGCCGTGGATCATCAGGTATTTCCCTTCATAACCATGCGCACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCG  >  NC_000913/245062‑245531
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
cgcTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGTGTAAATACTCATTTGCACGCTCGGCTGAc                                                                                                                                                                                                                             <  1:73391/251‑1 (MQ=255)
              cGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCATTTGCACGCTCGGCTGACCAAACACCAgcgc                                                                                                                                                                                                                <  1:94338/250‑1 (MQ=255)
                                    naGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGCAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTACCCCGGCTTCAGTTGATCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCATTTGCACGCACGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATAtc                                                                                                                                                                                          >  1:26636/2‑250 (MQ=255)
                                                                                                 acaGAAACGGTTGGTGGCTTACGGGTTTGCTAAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCATTTGCACGCTCGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGCCGATGGAAACACAATCGCCGTGGatcat                                                                                                                            >  1:89747/1‑251 (MQ=255)
                                                                                                                    tACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATGTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCATGTGCACGCTCGGCTGACCAAACTCCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGcc                                                                                                                                                       <  1:37145/205‑1 (MQ=255)
                                                                                                                                                       aGTTGCTCCCAGAAGTCCTTAAAGTTGAAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCATTTGCACGCTCGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGCCGATGGAAACACAATCGCCGTGGATCATCAGGTATTTcc                                                                                                                 <  1:211536/208‑1 (MQ=255)
                                                                                                                                                                                                                    tGCGGAACGGGTAAATACTCATTTGCACGCTCGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGCCGATGGAAACACAATCGCCGTGGATCATCATGTATTTCCCTTCATAACCATgcnc                                                                                                 <  1:198533/163‑1 (MQ=255)
                                                                                                                                                                                                                           ctggTAAATACTCATTTGCACGCTCGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGCCGATGGAAACACAATCGCCGTGGATCATCAGGTATTTCCCTTCATAACCATGCGCACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCg  >  1:233611/3‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
CGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACCACGTAACGACCATTGACCACAGAAACGGTTGGTGGCTTACGGGTTTGCTCAAAGTAGTCGTAACCCGGCTTCAGTTGCTCCCAGAAGTCCTTAAAGTTGGAATATTTATGGCGCTTCATATTGGCGTCGGTCATGCGGAACGGGTAAATACTCACTTGCACGCTCGGCTGACCAAACACCAGCGCACCAGTAACGAACTGGAATATCTCATCAATACCCTGATTGGTCATTGCGTAGCAGCCGATGGAAACACAATCGCCGTGGATCATCAGGTATTTCCCTTCATAACCATGCGCACGGTCATAGGCATTGGGGAAACCAATATTAATCGCTTTGTAGTAACGGCTGTCTGGTTTTAACTGATTACGCTGGACGCTATAAAACCCTTCCG  >  NC_000913/245062‑245531

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