Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,655,558	A→G	Y79Y (TAT→TAC)	yifK ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,655,558	0	A	G	94.2%	262.4 / 3.7	87	Y79Y (TAT→TAC)	yifK	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (5/0);  new base G (0/82);  total (5/82)
Fisher's exact test for biased strand distribution p-value = 2.71e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GGTGAGATAGCCAAAGAACGGGCTCATATAACGATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGA  >  W3110S.gb/3655519‑3655589                                        |                                ggTGAGATAACCAAAGAATGGGCTCATATAAGGATGCGCGTAAAc                            <  1:592135/45‑1 (MQ=255) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2680503/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:435269/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:342773/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:319862/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2925586/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2886521/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2877793/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2840612/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2771183/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2713013/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:495239/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2669345/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2622365/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2562712/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2457087/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2396614/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2387214/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2278303/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1011859/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2243903/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:762992/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:992186/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:983581/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:964116/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:96228/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:886043/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:873578/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:87301/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:842017/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:837263/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:488806/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:737600/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:708856/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:705820/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:687690/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:651837/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:606930/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:601844/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:552150/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2252004/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1284654/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1586292/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2212297/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:15631/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1543212/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1499530/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1495102/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1405540/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1374699/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1351454/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1326749/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1583985/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1273162/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1228530/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1215861/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1100293/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:109777/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1073496/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1058293/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:104276/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:102548/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1918231/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:213461/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:209812/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2058033/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2033321/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2021377/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:2017885/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1995677/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1979708/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1590250/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1897352/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1895433/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1888458/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1771645/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1715892/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1702697/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1675691/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:1653356/45‑1 (MQ=38) ggTGAGATAACCAAAGAATGCGCTCATATAACGATGCGCGTAAAc                            <  1:1634536/45‑1 (MQ=25)  gTGAGATAACCAAAGAATGGGCTCATATAACGATGCGCGTAAAc                            <  1:898457/44‑1 (MQ=37)           ccAAAGAACGGGCTCATATAACTATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGa  >  1:75946/1‑61 (MQ=255)           ccAAAGAACGGGCTCATATAACTATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGa  >  1:2426906/1‑61 (MQ=255)           ccAAAGAACGGGCTCATATAACGATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGa  >  1:1374946/1‑61 (MQ=255)           ccAAAGAACGGGCTCATATAACGATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGa  >  1:2605027/1‑61 (MQ=255)           ccAAAGAACGGGCTCATATAACGATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGa  >  1:578201/1‑61 (MQ=255)                                        |                                GGTGAGATAGCCAAAGAACGGGCTCATATAACGATGCGCATAAACGGCGAACGAACCGGTAACCGGTTCGA  >  W3110S.gb/3655519‑3655589

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