Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,620,788	A→G	62.5%	E312E (GAA→GAG)	purM →	phosphoribosylaminoimidazole synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,620,788	0	A	G	62.5%	9.1 / 11.7	16	E312E (GAA→GAG)	purM	phosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (7/3);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATC  >  W3110S.gb/2620728‑2620848                                                             |                                                             tGAAATGTATCGCACCTTCGACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAg                                                             >  1:2019567/1‑62 (MQ=255)   aaaTGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGc                                                   <  1:622414/70‑1 (MQ=255)         aTCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGcc                                            >  1:974167/1‑71 (MQ=255)                   cAACTGCGGCGTCGGGATGATTATTGCCTTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATg                                  >  1:1326264/1‑71 (MQ=255)                   cAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCt                                                >  1:2388387/1‑57 (MQ=255)                      cTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGcc                                >  1:2016226/1‑70 (MQ=255)                              tCGGGATGATTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGaa                       <  1:2347131/71‑1 (MQ=255)                                    tgatTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGc                                        >  1:2485083/1‑48 (MQ=255)                                    tgatTATTGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGaaaa                     >  1:2577296/1‑67 (MQ=255)                                       ttattGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTc                                      <  1:784432/47‑1 (MQ=255)                                        tattGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTgg              >  1:1446413/1‑70 (MQ=255)                                         attGCCCTGCCTGCTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGAGGaa            >  1:1257701/1‑71 (MQ=255)                                             cccTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATc        <  1:2265277/71‑1 (MQ=255)                                               ctgcctgcTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCgg      >  1:879507/1‑71 (MQ=255)                                                  cctgcTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTat   >  1:682423/1‑71 (MQ=255)                                                   ctgcTCCGGAGGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTatc  <  1:1414417/71‑1 (MQ=255)                                                             |                                                             TGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATC  >  W3110S.gb/2620728‑2620848

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