Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,897,032	+A	coding (653/999 nt)	yiaK ←	2,3‑diketo‑L‑gulonate dehydrogenase, NADH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,897,028	1	.	A	100.0%	56.5 / NA	22	L219F (TTG→TTT)	yiaK	2,3‑diketo‑L‑gulonate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (0/22);  total (0/22)

TTTCCAGTAGCCCATCGGC‑AAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCA  >  W3110S.gb/3897010‑3897071                    |                                            tttCCAGTAG‑CCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:843628/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:261424/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:95180/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:894470/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:892525/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:766624/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:519178/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:498118/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:402268/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:395982/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:348602/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1003456/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1444269/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1376224/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1351618/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1199376/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1173875/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1150978/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1094443/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1081028/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1065842/62‑1 (MQ=255)  ttCCAGTAGCCCATCGGCAAAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCa  <  1:1006628/62‑1 (MQ=255)                    |                                            TTTCCAGTAGCCCATCGGC‑AAAATGCGGCGATTCTTCTCGATAACGCCAGGTTCTTTGGTCA  >  W3110S.gb/3897010‑3897071

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