Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,641,728	A→G	N64N (AAT→AAC)	dapF ←	diaminopimelate epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,641,728	0	A	G	91.7%	31.1 / ‑3.6	12	N64N (AAT→AAC)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (0/11);  total (1/11)
Fisher's exact test for biased strand distribution p-value = 8.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCAGCGCGCACCGTTGCCGCACTGCGCCACTTCACTGCCATCAGCATTGAAAATGCGATAGTGAAAATCCAGTTCAGGATCATACGGCGGCTCAACCACCAGCAG  >  W3110S.gb/3641683‑3641787                                              |                                                            gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:1004937/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:1606750/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:1631987/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:205114/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:2442614/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:279859/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:2979038/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:3027768/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:561243/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:885057/62‑2 (MQ=255) gcagcGTGCGCCGTTGCCGCACTGCTCGACCTCGGAGCCGTCGGCGTTGAAGATCCGGTAgc                                             <  1:934786/62‑2 (MQ=255)                                            gCATTGAAAATGCGATAGTGAAAATCCAGTTCAGGATCATACGGCGGCTCAACCACcagcag  >  1:1028472/1‑62 (MQ=255)                                              |                                                            GCAGCGCGCACCGTTGCCGCACTGCGCCACTTCACTGCCATCAGCATTGAAAATGCGATAGTGAAAATCCAGTTCAGGATCATACGGCGGCTCAACCACCAGCAG  >  W3110S.gb/3641683‑3641787

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