Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,320,034	Δ2 bp	coding (344‑345/600 nt)	yiiQ →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,320,032	0	A	.	93.2%	173.5 / 6.9	44	coding (342/600 nt)	yiiQ	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (41/0);  total (41/3)
Fisher's exact test for biased strand distribution p-value = 7.55e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
*	minE	2,320,033	0	G	.	93.2%	167.2 / 9.9	44	coding (343/600 nt)	yiiQ	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base . (41/0);  total (41/3)
Fisher's exact test for biased strand distribution p-value = 7.55e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCCAGAGCAAAAAGCCGCGAAAGCGAAA  >  minE/2319981‑2320056                                                    ||                        tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCGTCCAGGGGCC‑‑AGCAAAAAGcc              >  1:190666/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:408510/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2621799/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2656446/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2725551/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2745743/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2919801/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:3086683/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:398910/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1075204/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:419623/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:435879/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:531595/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:60029/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:735870/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:827176/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:916467/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2268099/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1098119/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:216228/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1105159/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:2059634/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1238221/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1907401/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1856224/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:125292/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1842663/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1265603/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1837667/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1757580/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1752362/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:1484834/1‑62 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGca              >  1:2433408/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:3003181/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:485344/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:490456/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:2673526/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:2497550/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:839813/1‑61 (MQ=255) tACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGc               >  1:948300/1‑61 (MQ=255) tACCTTAAAAATCAAAACCATTCAAATGACCTGGCTACCCATCCAGGGGCC‑‑AGCAAAAAGcc              >  1:805766/1‑62 (MQ=255)                aaGCATTCAAATGACCTGGCTACCCATCCAGGGGCCAGAGCAAAAAGCCgcgaaagcgaaa  <  1:286029/61‑1 (MQ=255)                aaGCATTCAAATGACCTGGCTACCCATCCAGGGGCCAGAGCAAAAAGCCgcgaaagcgaaa  <  1:2631394/61‑1 (MQ=255)                aaGCATTCAAATGACCTGGCTACCCATCCAGGGGCCAGAGCAAAAAGCCgcgaaagcgaaa  <  1:1030452/61‑1 (MQ=255)                                                    ||                        TACCTTAAAAATCAAAAGCATTCAAATGACCTGGCTACCCATCCAGGGGCCAGAGCAAAAAGCCGCGAAAGCGAAA  >  minE/2319981‑2320056

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