Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,125,199	+G	coding (1479/1500 nt)	glcD ←	glycolate oxidase subunit, FAD‑linked

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,125,196	1	.	G	82.7%	146.5 / 24.9	52	coding (1482/1500 nt)	glcD	glycolate oxidase subunit, FAD‑linked
Reads supporting (aligned to +/- strand):  ref base . (0/9);  new base G (43/0);  total (43/9)
Fisher's exact test for biased strand distribution p-value = 2.72e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GCGTAGCATCAGAAACGCTCCAGTTCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTA  >  W3110S.gb/3125170‑3125254                            |                                                           gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCAcatgcatg                           >  1:2725962/1‑61 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCAcatgcatg                           >  1:2742352/1‑61 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:556149/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2413997/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1027592/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2454229/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2480021/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2578080/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2584068/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2687511/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2825358/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:391912/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:485331/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2449875/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:651792/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:657945/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:807157/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:850777/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:854062/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:883921/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:888617/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:89965/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:907124/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:986948/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:999387/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1802386/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1327385/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1340580/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1346277/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1430389/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1446595/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1234710/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1586013/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1605663/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1789435/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1202103/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1157006/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:115386/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:1972617/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:11450/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2063190/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGTGATGCACATGCATgg                          >  1:2361867/1‑62 (MQ=255) gcgTAGCATCAGAAACGCTCCAGTTCAGGGGAAAGGTAAATGACCGAGATGCAcatgcatg                           >  1:262991/1‑61 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:1084762/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:918683/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:1560657/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:1827032/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:18481/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:1904895/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:2566984/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:1997654/62‑1 (MQ=255)                        ttCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTa  <  1:2088093/62‑1 (MQ=255)                            |                                                           GCGTAGCATCAGAAACGCTCCAGTTCA‑GGGAAAGGTAAATGACCGTGATGCACATGCATGGCACCAAATTCAGCACAGCGGTGTA  >  W3110S.gb/3125170‑3125254

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