Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	61,620	A→T	W549R (TGG→AGG)	hepA ←	RNA polymerase‑associated helicase protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	61,620	0	A	T	100.0%	80.3 / NA	30	W549R (TGG→AGG)	hepA	RNA polymerase‑associated helicase protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (0/30);  total (0/30)

CAGCAGTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCAGGCGGCAGCGCGGTCACGTTCGA  >  W3110S.gb/61582‑61643                                       |                        cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTTGCGGCAGCGCGGTCACGTTCGa  <  1:1382065/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:429082/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:913039/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:747283/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:742402/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:719665/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:698606/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:683542/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:677300/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:613336/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:588049/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:557832/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:451606/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:434794/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:431030/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:100573/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:358341/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:244286/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:214893/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:200166/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:1348775/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:134173/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:1016351/62‑1 (MQ=255) cagcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:1014296/62‑1 (MQ=255)  agcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:275125/61‑1 (MQ=255)  agcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:1384797/61‑1 (MQ=255)  agcagTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:114327/61‑1 (MQ=255)  agcagTACCTGTGAGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:1279954/61‑1 (MQ=255)  agcagTACCCGTGCGCCGGTGTCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:714817/61‑1 (MQ=255)                      gtCTTCTTCGGCAAACCTGGCGGCAGCGCGGTCACGTTCGa  <  1:1589938/41‑1 (MQ=255)                                       |                        CAGCAGTACCTGTGCGCCGGTGTCTTCTTCGGCAAACCAGGCGGCAGCGCGGTCACGTTCGA  >  W3110S.gb/61582‑61643

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