Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,073,671	C→T	100%	G50S (GGC→AGC)	ydhQ ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,073,671	0	C	T	80.0%	69.2 / 12.1	35	G50S (GGC→AGC)	ydhQ	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/0);  new base T (28/0);  total (35/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.86e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTC  >  minE/1073620‑1073695                                                    |                         tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCCCTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:761138/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGGCCCCGTGACCAc           >  1:1540774/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCa            >  1:807026/1‑66 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCa            >  1:583096/1‑66 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCa            >  1:439012/1‑66 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:390528/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:992287/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:975913/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:924267/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:881776/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:751489/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:704590/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:626293/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:591252/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:561107/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1039945/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1949865/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1041353/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1077746/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1079267/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1175399/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1269783/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1321265/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:1945415/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:306307/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:2213946/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:2308406/1‑67 (MQ=255) tACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCTAGCCCCCGTGACCAc           >  1:26975/1‑67 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtttt   >  1:236840/1‑66 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTAtttt   >  1:407997/1‑66 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTc  >  1:2072956/1‑67 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTc  >  1:525616/1‑67 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTc  >  1:117793/1‑67 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTc  >  1:458036/1‑67 (MQ=255)          gtgtgCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTc  >  1:920121/1‑67 (MQ=255)                                                    |                         TACCGTACTGTGTGCATTTGGCGATCCCACCGGTTGCCACTTTCAACCAGCCAGCCCCCGTGACCACGGTATTTTC  >  minE/1073620‑1073695

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