Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	36,728	Δ1 bp	coding (66/1218 nt)	caiB ←	crotonobetainyl CoA:carnitine CoA transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	36,728	0	G	.	100.0%	243.6 / NA	55	coding (66/1218 nt)	caiB	crotonobetainyl CoA:carnitine CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (0/55);  total (0/55)

ACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCG  >  minE/36712‑36774                 |                                               aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCATTCCGGCCAACGGCCCGAATTTCg  <  1:159324/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:744702/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:511262/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:527257/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:550710/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:558962/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:577600/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:595284/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:604111/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:61002/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:632380/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:669030/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:684217/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:708869/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:716052/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:502425/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:757852/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:761336/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:764115/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:773670/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:846995/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:884079/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:895718/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:897367/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:913199/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:923627/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:951204/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:9996/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:287556/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:1009816/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:1022473/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:104968/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:141948/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:158783/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:159399/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:165640/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:194588/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:208404/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:252546/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:267676/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:285774/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:1008623/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:348311/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:359200/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:381332/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:41413/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:417563/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:441993/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:470194/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:475064/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:479746/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:485451/62‑1 (MQ=255) aCGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCAGGCCAACGGCCCGAATTTCg  <  1:601736/62‑1 (MQ=255)  cGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:408699/61‑1 (MQ=255)  cGGTCCGGCGATTTC‑ATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCg  <  1:12465/61‑1 (MQ=255)                 |                                               ACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCG  >  minE/36712‑36774

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