Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	802,452	+T	coding (601/1071 nt)	dhaK ←	dihydroxyacetone kinase, N‑terminal domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	802,449	1	.	T	96.7%	220.0 / 0.7	61	coding (604/1071 nt)	dhaK	dihydroxyacetone kinase, N‑terminal domain
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (59/0);  total (61/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.31e-01

AATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGG‑T‑TTGCCCGCGGCAGGAA  >  minE/802402‑802466                                                 |                   aaTGCCGACGCCAAACTCCATCTCATTATCCGCCATGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:902515/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCggaaggaa  >  1:700260/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:6152/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:220249/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1000382/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:291459/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:334642/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:354514/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:378350/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:461311/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:474994/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:503564/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:539767/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:562469/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:593472/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:603812/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:606889/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:26164/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:641225/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:656285/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:670094/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:671851/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:694746/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:796438/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:81139/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:817389/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:88186/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:922871/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:940471/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:979046/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1549302/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1007481/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1010520/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1021524/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1026811/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1057571/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1067664/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1096177/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1148061/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:115263/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1153364/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1157701/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1168104/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1216758/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1224794/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1236886/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1291829/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1297200/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:130246/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1343707/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1356124/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1392780/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1428251/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1477536/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1537159/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1539884/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:1546433/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGaa  >  1:156985/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGGTT‑TTGCCCGCGGCAGGa   >  1:608164/1‑65 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGG‑t                   >  1:476266/1‑49 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGG‑TGTTGCCCGCGGCAGGaa  >  1:1431519/1‑66 (MQ=255) aaTGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGG‑TCTTGCCCGCGGCAGGaa  >  1:430652/1‑66 (MQ=255)                                                 |                   AATGCCGACGCCAAACTCCATCTCATTATCCGCCAGGGTAAAAGAAGG‑T‑TTGCCCGCGGCAGGAA  >  minE/802402‑802466

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: ‑