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: