Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,933,009	+A	coding (415/1179 nt)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,933,006	1	.	A	100.0%	85.9 / NA	25	E138K (GAA→AAA)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (0/25);  total (0/25)

TGCCCACTTCCACTTATCGTTTTGCCGATAGCG‑AAAGCCTTTTCCGCGAGGCGGTTGCTGACAT  >  W3110S.gb/1932974‑1933037                                  |                                tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1026795/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:954998/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:917298/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:770819/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:508984/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:421152/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:381115/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:2492444/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:2268854/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:2203944/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:2194780/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:2141859/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1987339/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1945016/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1797860/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1626530/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:142961/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1341391/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1340143/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1306134/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1064831/65‑1 (MQ=255) tgcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1014947/65‑1 (MQ=255) tgcCCACTGCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:2051319/65‑1 (MQ=255)  gcCCACTTCCACTTATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:1524767/64‑1 (MQ=255)    ccacttccacttATCGTTTTGCCGATAGCGAAAAGCCTTTTCCGCGAGGCGGTTGCTGACAt  <  1:718894/62‑1 (MQ=255)                                  |                                TGCCCACTTCCACTTATCGTTTTGCCGATAGCG‑AAAGCCTTTTCCGCGAGGCGGTTGCTGACAT  >  W3110S.gb/1932974‑1933037

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