Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,128,917	A→G	A137A (GCT→GCC)	fcl ←	bifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase and GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,128,917	0	A	G	81.8%	70.0 / 9.7	33	A137A (GCT→GCC)	fcl	bifunctional GDP‑fucose synthetase: GDP‑4‑dehydro‑6‑deoxy‑D‑mannose epimerase and GDP‑4‑dehydro‑6‑L‑deoxygalactose reductase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (27/0);  total (29/4)
Fisher's exact test for biased strand distribution p-value = 3.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GTCCGTACTGGCGGTTGTATGATTCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATAAGGCTCGTTAGTCGGCTCCAGCGT  >  W3110S.gb/2128862‑2128946                                                        |                              gTCCGTACTGGCGGTTGTATGATTCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATaa                         <  1:1165521/62‑1 (MQ=255) gTCCGTACTGGCGGTTGTATGATTCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATaa                         <  1:1835807/62‑1 (MQ=255) gTCCGTACTGGCGGTTGTATGATTCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATaa                         <  1:1544492/62‑1 (MQ=255) gTCCGTACTGGCGGTTGTATGATTCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATaa                         <  1:1301091/62‑1 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGc                             >  1:1568920/1‑56 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGc                             >  1:432053/1‑56 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:564508/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:555994/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:680222/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:24209/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:482668/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:683182/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:715995/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:337547/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:3092943/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:3070864/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:3016393/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:2934862/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:2911243/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:243583/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:1162299/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:2186703/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:217206/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:195711/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:1832938/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:1783185/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:145051/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:142096/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:1228747/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:1225420/1‑57 (MQ=255)   ccGTACTGGCGGTTGTAAGATTCACACAGTTTGATCCCGGCGATTTTGGCAATGGCa                            >  1:1210372/1‑57 (MQ=255)                        tCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATAAGGCTCGTTAGTCGGCTCCAGCGt  >  1:355541/1‑62 (MQ=255)                        tCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATAAGGCTCGTTAGTCGGCTCCAGCGt  >  1:815904/1‑62 (MQ=255)                                                        |                              GTCCGTACTGGCGGTTGTATGATTCGCACAGTTTGATCCCGGCGATTTTGGCAATAGCATAAGGCTCGTTAGTCGGCTCCAGCGT  >  W3110S.gb/2128862‑2128946

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