BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_012967 649,391 T→A 100% I471F (ATC→TTC)  mrdA ← peptidoglycan DD‑transpeptidase MrdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_012967649,3910TA100.0% 37.4 / NA 14I471F (ATC→TTC) mrdApeptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/7);  total (7/7)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGG > NC_012967/649362‑649420
|
aaCCTTCACGATACCGTCATTAATCAGGAACATCAg < 1:235356/36‑1 (MQ=37)
aaCCTTCACGATACCGTCACTAATCAGGAacatcaa > 1:451453/1‑35 (MQ=25)
cTTCACGATACCGTCATTAATCAGGAACATCAGTGc < 2:543819/36‑1 (MQ=37)
ttCACGATACCGTCATTAATCAGGAACATCAGTGcc < 1:570901/36‑1 (MQ=37)
ttCACGATACCGTCATTAATCAGGAACATCAGTGcc > 1:900870/1‑36 (MQ=37)
cACGTTACCGTCATTAATCAGGAACATCAGTGCCtt < 1:919479/36‑1 (MQ=25)
cACGATACCGTCATTAATCAGGAACATCAGTGCCtt < 1:138512/36‑1 (MQ=37)
cACGATACCGTCATTAATCAGGAACATCAGCGCCtt > 2:824486/1‑36 (MQ=21)
cGATATCGTCATTAATCAGGAACATCAGTGCCTTAc < 2:727904/36‑1 (MQ=25)
cGATACCGTCATTAGTCAGGAACATCAGTGCCTTAc > 1:465473/1‑36 (MQ=37)
cGATACCGTCATTAATCCGGAACATCAGTGCCttcc > 1:647247/1‑34 (MQ=25)
gATACCGTCATTAATCAGGAACATCAGTGCCTTACt < 1:269849/36‑1 (MQ=37)
ccGTCATTAATCAGGAACATCAGTGCCTTACTCATc > 2:806664/1‑36 (MQ=38)
aaTCAGGAACATCAGTGCCTTACTCATCTCGATTgg > 2:737523/1‑36 (MQ=37)
tCAGGAACATCAGTGCCTTACTCATCTGGATTgggg > 1:279810/1‑36 (MQ=38)
|
AACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGG > NC_012967/649362‑649420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: