Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,269,829 G→T Q168K (CAG→AAG)  murB ← UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,269,8290GT98.2% 199.0 / ‑3.6 56Q168K (CAG→AAG) murBUDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (55/0);  total (56/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGC  >  minE/2269815‑2269861
              |                                
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTcgcg       >  1:11620/2‑42 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGAt     >  1:394363/2‑44 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAg   >  1:693763/2‑46 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:467956/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1016442/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:203082/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:22249/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:28541/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:299101/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:32123/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:352433/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:36992/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:373436/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:442395/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:448235/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:448326/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:166368/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:611028/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:616694/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:66709/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:708480/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:744062/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:75511/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:790319/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:807958/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:809834/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:924419/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:981421/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1557823/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:103582/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1037318/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1052516/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:107247/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1119485/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:112277/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1134654/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1141212/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1164945/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1201144/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1207010/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1252285/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1283738/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1365121/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1399340/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1436724/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1454794/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:14911/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:151669/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1525274/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1526084/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1527870/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1545540/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:1559488/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTTAAAAAATACTGTCGCGATAGc  >  1:481786/2‑47 (MQ=255)
ggCGAAGCGATCCTTGTATTCATGTGTAAAAATACTGTCGCGATAGc  >  1:1068699/2‑47 (MQ=255)
ggCGAAGCGATCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGc  >  1:818533/2‑47 (MQ=255)
              |                                
AGCGAAGCGGTCCTGGTATTCATGTTTAAAAATACTGTCGCGATAGC  >  minE/2269815‑2269861

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