Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,472,315	+C	intergenic (‑594/‑102)	yfaL ← / → nrdA	adhesin/ribonucleoside diphosphate reductase 1, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,472,312	1	.	C	100.0%	145.0 / NA	40	intergenic (‑591/‑105)	yfaL/nrdA	adhesin/ribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (0/40);  total (0/40)

AAGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAAT‑CCCAAACATACC  >  minE/1472264‑1472324                                                  |             aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:60854/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:30054/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1049102/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:31526/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:359046/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:409935/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:411559/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:432801/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:575276/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:607212/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:313807/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:64470/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:658790/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:753141/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:807047/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:864246/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:867852/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:899890/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:971941/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:227888/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1085155/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1124144/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1125884/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1157068/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1164737/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:117004/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1175300/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1185101/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1190462/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:1198928/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:146426/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:157979/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:172424/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:180736/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:219094/62‑1 (MQ=255) aaGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:275371/62‑1 (MQ=255)  aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:632443/61‑1 (MQ=255)  aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:789427/61‑1 (MQ=255)  aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:20906/61‑1 (MQ=255)  aGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAATCCCCAAACATAcc  <  1:932360/61‑1 (MQ=255)                                                  |             AAGTTATGCACTTGCAAGAGGGTCATTTTCACACTATCTTGCAGTGAAT‑CCCAAACATACC  >  minE/1472264‑1472324

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