Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,764,766	Δ1 bp	coding (193/1074 nt)	yfjN →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,764,763	0	G	.	77.3%	40.5 / 18.6	22	coding (190/1074 nt)	yfjN	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base . (0/17);  total (5/17)
Fisher's exact test for biased strand distribution p-value = 3.80e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGG  >  W3110S.gb/2764718‑2764815                                              |                                                     accagaccagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACaa                                    >  1:1168828/1‑64 (MQ=255) accagaccagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACaa                                    >  1:1049742/1‑64 (MQ=255) accagaccagGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACaa                                    >  1:13174/1‑64 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:688328/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:508610/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:418905/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:416023/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:399786/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:331152/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:2514763/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:248666/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:244031/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1891501/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1484190/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1478983/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1398902/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1391937/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1368157/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTCGACTACAATTCAata                             <  1:1265507/65‑1 (MQ=255)      accagGGGTTAACTTGGCAACTGTGGACATTTTTTATAAC‑GGGACGGTTAGACTACAATTCAata                             <  1:270633/65‑1 (MQ=255)                                   ttttttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTggg  >  1:1507481/1‑64 (MQ=255)                                   ttttttATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTggg  >  1:1442604/1‑64 (MQ=255)                                              |                                                     ACCAGACCAGGGGTTAACTTGGCAACTGTGGACATTTTTTATAACGGGGACGGTTCGACTACAATTCAATATCTCACTGGAGCCAATCGTTCTCTGGG  >  W3110S.gb/2764718‑2764815

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