Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	598,434	A→G	100%	A259A (GCA→GCG)	dmsA →	dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	598,434	0	A	G	81.3%	135.6 / 24.5	64	A259A (GCA→GCG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base A (0/12);  new base G (52/0);  total (52/12)
Fisher's exact test for biased strand distribution p-value = 3.04e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCC  >  minE/598387‑598452                                                |                   gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTATTGcc  <  1:299190/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:2232421/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:579303/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:61781/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:771574/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:1446258/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:1379762/66‑1 (MQ=255) gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:126503/66‑1 (MQ=255)  aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:1641929/65‑1 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2276882/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2265376/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:337608/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:354176/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2256765/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:37404/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:374678/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:554810/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:560007/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2254764/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:630754/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:609314/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:627184/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1072074/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:740877/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:741777/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:775376/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:789759/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:837095/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:851991/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:853751/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:875297/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:879644/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:893930/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2201940/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1016535/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1140764/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1152458/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1212659/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1329780/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1363758/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1475368/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1538960/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:15585/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1690522/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1707289/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:173210/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1738341/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1768389/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1770652/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1783646/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1830440/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1870636/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1960472/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1967294/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1741224/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:1981250/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2153407/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2192152/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2199547/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGaa            >  1:2242929/1‑52 (MQ=255)     cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCGCGCCAGa             >  1:981020/1‑51 (MQ=255)             aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:1454040/54‑1 (MQ=255)             aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:2273579/54‑1 (MQ=255)                          tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGcc  <  1:990099/41‑1 (MQ=255)                                                |                   GAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCC  >  minE/598387‑598452

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