Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,325,522 T→G 100% K71T (AAA→ACA)  wcaL ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,325,5220TG100.0% 57.2 / NA 16K71T (AAA→ACA) wcaLpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (0/16);  total (0/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGT  >  minE/1325468‑1325538
                                                      |                
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:1009614/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:1054722/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:1455894/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:1806970/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:1854123/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2017513/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2148852/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2166052/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2205175/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2268067/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2547001/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:2561940/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:284004/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:389887/71‑1 (MQ=255)
cAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:805205/71‑1 (MQ=255)
 aGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGt  <  1:1915132/70‑1 (MQ=255)
                                                      |                
CAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGT  >  minE/1325468‑1325538

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