Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,951,954	G→A	R902R (CGC→CGT)	recB ←	exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,951,954	0	G	A	100.0%	138.3 / NA	51	R902R (CGC→CGT)	recB	exonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (51/0);  total (51/0)

AATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGCTGGTGACGCGCCAGTTATCG  >  W3110S.gb/2951902‑2951966                                                     |             gagccTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:19769/4‑62 (MQ=255) gaTCCTGAGCGATACCGTGCCCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1210447/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAACCCGGAGTAGCTGGTGACACGCCAGtt      >  1:1394955/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:414570/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:1325500/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:1445788/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:62678/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:653441/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:1255338/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:678614/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:1068955/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGtt      >  1:913593/2‑61 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:399594/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:408117/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:385576/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:357002/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:347694/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:996546/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:334251/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:468736/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:5669/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:614032/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:669674/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:760767/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:799706/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:819994/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:893631/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:91366/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1005695/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1304514/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1006620/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1034528/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1043288/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:10970/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1131341/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:114201/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1203513/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1217204/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1256012/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:268455/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1309290/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1364512/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1384427/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:1386147/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:156431/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:161020/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:176838/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:238195/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACGCCAGTTa     >  1:251096/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGCTGCTGCAAACCGGAGTAGCTGGTGACACCCCAGTTa     >  1:1226342/2‑62 (MQ=255) gaTCCTGAGCGATACCGTGACCACGC‑‑‑TGCAAACCGGAGTAGCTGGTGAAACGCCAGTTATCg  >  1:413921/2‑62 (MQ=255)                                                     |             AATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTAGCTGGTGACGCGCCAGTTATCG  >  W3110S.gb/2951902‑2951966

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