Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,551,604 G→T 100% A71S (GCT→TCT)  yibD → predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,551,6040GT100.0% 147.4 / NA 43A71S (GCT→TCT) yibDpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (43/0);  total (43/0)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTGCATCAGGCGAATGCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGAA  >  minE/2551589‑2551639
               |                                   
ttgCATCAGGCGAATTCTGGCGCATCGGTGGGGCGTAATCGTGGGATTGaa  >  1:1613549/1‑51 (MQ=39)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:424279/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1991818/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1992191/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:2058383/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:216815/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:2184091/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:2218197/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:230190/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:302396/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:376340/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:42245/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1838835/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:501789/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:636831/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:694220/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:753839/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:787941/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:939186/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:95051/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:962543/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:977522/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1832174/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1130583/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1173830/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1273879/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1308254/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1371363/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1378528/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1570083/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1601033/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1657163/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1664757/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1665894/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1685299/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1694162/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1715753/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1739513/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1742242/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1764702/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:1023107/1‑51 (MQ=255)
ttgCATCAGGCGAATTCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGa   >  1:92410/1‑50 (MQ=255)
ttgCATCAGGCGAATTCTGGAGCATCGGTGGCGCGTAATCGTGGGATTGaa  >  1:166411/1‑51 (MQ=255)
               |                                   
TTGCATCAGGCGAATGCTGGCGCATCGGTGGCGCGTAATCGTGGGATTGAA  >  minE/2551589‑2551639

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