Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,371,684	G→A	A249V (GCG→GTG)	nanT ←	sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,371,684	0	G	A	76.2%	33.4 / 8.4	21	A249V (GCG→GTG)	nanT	sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (16/0);  total (19/2)
Fisher's exact test for biased strand distribution p-value = 4.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTC  >  W3110S.gb/3371654‑3371743                               |                                                            gtgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:1679906/3‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:2729476/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:972411/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:9400/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:706324/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:615238/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:388635/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:373975/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:2877830/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:1205003/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:2643313/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:2413230/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:2068418/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:1976654/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:1472278/1‑61 (MQ=255) gcgcATAACAGCCCAAGAACAGCGACGATCACGGCATTTTGCAGGTTACCGGCGAAGCAGa                               >  1:1468587/1‑61 (MQ=255)               aaGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAg                >  1:1982568/1‑62 (MQ=255)               aaGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAg                >  1:928051/1‑62 (MQ=255)               aaGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAg                >  1:1404077/1‑62 (MQ=255)                             tCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTc  <  1:2985904/62‑1 (MQ=255)                             tCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTc  <  1:761761/62‑1 (MQ=255)                               |                                                            GCGCATAACAGCCCAAGAACAGCGACGATCGCGGCATTTTGCAGGTTACCGGCGAAGCAGAACCACAGCGCAGTAGCCGCCGCCAGTGTC  >  W3110S.gb/3371654‑3371743

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