Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,581,155	G→A	G93G (GGC→GGT)	ypfG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,581,155	0	G	A	94.0%	243.0 / 2.3	100	G93G (GGC→GGT)	ypfG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (0/94);  total (6/94)
Fisher's exact test for biased strand distribution p-value = 8.39e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCCTCCGGCGACTT  >  W3110S.gb/2581116‑2581208                                        |                                                      tGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGgc                                  >  1:1381881/1‑61 (MQ=255) tGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGgc                                  >  1:353098/1‑61 (MQ=255) tGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGgc                                  >  1:984267/1‑61 (MQ=255) tGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGgc                                  >  1:590273/1‑61 (MQ=255) tGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGgc                                  >  1:907395/1‑61 (MQ=255) tGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGgc                                  >  1:575774/1‑61 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGTCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1365488/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCTCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:809544/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:515850/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:625911/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:622100/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:582573/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:549809/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:544339/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:540613/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:516524/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:400817/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:507140/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:491325/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:464611/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:453398/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:421364/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:416364/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:400894/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1068202/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:890103/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:996383/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:991451/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:980642/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:969331/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:924867/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:915061/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:910621/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:905964/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:905113/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:678057/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:889247/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:879634/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:866852/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:797372/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:795685/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:778177/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:70131/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:696048/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:69339/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1166967/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1414088/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1386038/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1382485/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:13286/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:129402/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1276426/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1226381/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1214954/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1209039/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:36195/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1157415/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1133007/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1130158/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1130100/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1115254/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1094879/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1086170/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1019102/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1427732/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:387291/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:332990/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:319078/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:264477/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:263759/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:238762/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:233193/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1464777/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGAAGCATCCGGCGACtt  <  1:221139/62‑1 (MQ=255)                                aaCGGCTCACCATCTAACAGCAGCAGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:938738/62‑1 (MQ=255)                                aaCGGCGCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:710881/62‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1470758/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:602496/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:987617/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:986774/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:322950/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:674173/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:446701/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:937248/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:308753/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:489899/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:30251/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1372299/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1355516/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1230390/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1265214/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:812395/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:236665/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1331098/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1294268/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:760962/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:749932/61‑1 (MQ=255)                                 aCGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:715915/61‑1 (MQ=255)                                  cGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:1215879/60‑1 (MQ=255)                                  cGGCTCACCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCATCCGGCGACtt  <  1:303177/60‑1 (MQ=255)                                        |                                                      TGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAGCCGTGGCGCTATCTCCCCTTCTGACGCCTCCGGCGACTT  >  W3110S.gb/2581116‑2581208

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