Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,286,501	C→A	60.0%	A261S (GCG→TCG)	yejK ←	nucleotide associated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,286,501	0	C	A	60.0%	‑5.5 / 27.0	20	A261S (GCG→TCG)	yejK	nucleotide associated protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (3/9);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTC  >  W3110S.gb/2286450‑2286556                                                    |                                                        tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACt                                      >  1:203280/1‑71 (MQ=255) tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACt                                      >  1:1657917/1‑71 (MQ=255) tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACt                                      >  1:1520722/1‑71 (MQ=255) tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACt                                      >  1:1355261/1‑71 (MQ=255)        gTCTGCCGGGTAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                 <  1:262545/69‑1 (MQ=255)        gTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                 <  1:1361286/69‑1 (MQ=255)                          ttCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTc                                    <  1:1059183/48‑1 (MQ=255)                          ttCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTc                                    <  1:1647575/48‑1 (MQ=255)                            ccAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTctc                                  >  1:987366/1‑46 (MQ=255)                            ccAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGc                                  >  1:12022/1‑48 (MQ=255)                            ccAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGc                                  >  1:115399/1‑48 (MQ=255)                               gTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc              <  1:1304099/65‑1 (MQ=255)                               gTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc              <  1:1408202/65‑1 (MQ=255)                               gTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc              <  1:1236366/65‑1 (MQ=255)                               gTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc              <  1:1198251/65‑1 (MQ=255)                               gTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc              <  1:616340/65‑1 (MQ=255)                                      gCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTc  <  1:615128/70‑1 (MQ=255)                                      gCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTc  <  1:1010655/70‑1 (MQ=255)                                      gCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTc  <  1:1128027/70‑1 (MQ=255)                                      gCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTc  <  1:977840/70‑1 (MQ=255)                                                    |                                                        TGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTC  >  W3110S.gb/2286450‑2286556

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