Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,326,272	G→A	D35D (GAC→GAT)	phnG ←	carbon‑phosphorus lyase complex subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,326,272	0	G	A	90.5%	115.6 / 3.9	42	D35D (GAC→GAT)	phnG	carbon‑phosphorus lyase complex subunit
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (38/0);  total (39/3)
Fisher's exact test for biased strand distribution p-value = 3.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGC  >  W3110S.gb/4326236‑4326316                                     |                                             caCCAGGCCAGTTTCATCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2729445/2‑59 (MQ=255) caCCAGGCCAGTTTCATCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:1312875/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGtt                                >  1:118655/2‑51 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGt                                 >  1:1398750/2‑50 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:547451/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:290454/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2909585/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2930417/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2950406/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2998977/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:3058396/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:469216/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:473128/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:546563/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:991530/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:636555/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:685951/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:755054/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:853262/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:960348/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:1054139/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2274207/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:1156590/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:1867664/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2115271/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2157655/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2270107/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2677047/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2316772/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:231950/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:235667/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:238568/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:242253/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2491288/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2494442/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2546985/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                        >  1:2557674/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTAAGCGCGt                        >  1:916582/2‑59 (MQ=255)   ccAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGtt                       <  1:740838/58‑1 (MQ=255)   ccAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGtt                       <  1:80745/58‑1 (MQ=255)   ccAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGtt                       <  1:1780953/58‑1 (MQ=255)                    gcgcGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGc  >  1:567627/1‑62 (MQ=255)                                     |                                             TACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGC  >  W3110S.gb/4326236‑4326316

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