Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,789,595	+A	coding (958/1269 nt)	ygaF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,789,591	1	.	A	82.3%	40.7 / 5.6	17	coding (954/1269 nt)	ygaF	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (0/14);  total (1/16)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CAGGACTGGGCGAGATGAAAAACTCGCTGTGC‑AAAAGCGGCTATCTGCGGCTGGTGCAAAAGTATTGTCCCCGGCTTTCGTTAAGCGA  >  W3110S.gb/2789560‑2789647                                 |                                                         cAGGACTGGGCGAGATGAAAAACTCGCTGTGC‑AAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:2080936/62‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:865888/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1114703/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:840222/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:819736/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:2685966/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:2537615/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:2534255/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1863433/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1584076/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1558727/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1471391/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1356999/60‑1 (MQ=255)    gACTGGGCGAGATGAAAAACTCGCTGTGCAAAAAGCGGCTATCTGCGGCTGGTGCAAAAg                            <  1:1125744/60‑1 (MQ=255)                  aaaaaCTCGCTGTGC‑AAAAGCGGCTATCTGCGGCTGGTGCAAAAGTATTGTCCCCGGCTTTc           <  1:2899028/62‑1 (MQ=255)                   aaaaCTCGCTGTGCAAAAAGCGTCTATCTGCGGCTGGTGCATAAg                            <  1:2684848/45‑1 (MQ=255)                           cTGTGC‑AAAAGCGGCTATCTGCGGCTGGTGCAAAAGTATTGTCCCCGGCTTTCGTTAAGCGa  >  1:13829/1‑62 (MQ=255)                                 |                                                         CAGGACTGGGCGAGATGAAAAACTCGCTGTGC‑AAAAGCGGCTATCTGCGGCTGGTGCAAAAGTATTGTCCCCGGCTTTCGTTAAGCGA  >  W3110S.gb/2789560‑2789647

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