Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,609,944	G→A	R134R (CGC→CGT)	fhuF ←	ferric iron reductase involved in ferric hydroximate transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,609,944	0	G	A	94.0%	148.0 / 0.2	50	R134R (CGC→CGT)	fhuF	ferric iron reductase involved in ferric hydroximate transport
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (0/47);  total (0/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

AACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGCCTTTTC  >  W3110S.gb/4609934‑4610003           |                                                            aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:428837/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1040069/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2442919/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2520714/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2703374/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2704456/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2795287/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2801336/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2849072/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2879907/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:363064/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:400393/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2384443/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:524951/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:530124/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:58047/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:640845/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:643723/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:751625/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:763209/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:766542/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:815008/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:878732/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:970164/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2040131/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1113153/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1441275/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1473614/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1504338/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1579323/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1593146/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1610750/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1733394/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1861170/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:1918022/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2012799/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2128622/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2139792/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2151823/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2196374/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2207881/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2211229/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2228096/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2269239/62‑1 (MQ=255) aaCAGGCGGCACGGCCGGTTTCGTGAAACTCAGCAAGGAAATGTTCCGGCGACACATCTAAt          <  1:2155202/62‑1 (MQ=255) aaCAGGCAGCACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:2339813/62‑1 (MQ=255)     ggcggcACGGCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAAt          <  1:369008/58‑1 (MQ=255)         aCGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGCCTTTTc  <  1:2072468/62‑1 (MQ=255)         aCGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGCCTTTTc  <  1:1967357/62‑1 (MQ=255)         aCGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGCCTTTTc  <  1:2603898/62‑1 (MQ=255)           |                                                            AACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGCCTTTTC  >  W3110S.gb/4609934‑4610003

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