Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	204,556	A→G	intergenic (+156/‑28)	alaV → / → rrlH	tRNA‑Ala/23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	204,556	0	A	G	100.0%	120.6 / NA	40	intergenic (+156/‑28)	alaV/rrlH	tRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/35);  total (5/35)

GAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCG  >  minE/204525‑204600                                |                                             gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:614498/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:3183857/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:3240934/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:3476415/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:3574660/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:3631370/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:4427/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:532783/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:580203/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:60916/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:118589/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:664481/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:719976/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:739217/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:764471/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:814486/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:82901/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:839199/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:863504/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:956465/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:1003814/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2956807/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2926726/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2867118/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2805153/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:1403977/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2306121/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:228776/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2178321/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2143688/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:1456175/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:1704164/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:1810877/61‑1 (MQ=25) gAGTCTCTCAAATTTTCGCAACACGATGATGGATCACAAGAAACATCTTCGGGTTGTGAgg                 <  1:2973494/61‑1 (MQ=255)     ctctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAgg                 <  1:2780583/57‑1 (MQ=255)               ttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg  >  1:1883292/1‑62 (MQ=25)               ttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg  >  1:1958717/1‑62 (MQ=25)               ttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg  >  1:1416740/1‑62 (MQ=25)               ttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg  >  1:3079/1‑62 (MQ=25)               ttCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg  >  1:992479/1‑62 (MQ=25)                                |                                             GAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCG  >  minE/204525‑204600

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