Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,827,842	Δ1 bp	coding (1085/2562 nt)	mutS →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,827,840	0	C	.	100.0%	85.3 / NA	19	coding (1083/2562 nt)	mutS	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (19/0);  total (19/0)

GTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCCACGCGATCTGGCCCGTATG  >  minE/1827790‑1827861                                                   |                      gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGcca                     >  1:91301/1‑52 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:439925/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:935044/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:817732/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:802251/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:676634/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:66165/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:588354/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:558466/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:541724/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1094085/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:434261/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:241769/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:223621/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1531251/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1469557/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1427238/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1374760/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1344417/1‑71 (MQ=255) gTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCG‑CCACGCGATCTGGCCCGTATg  >  1:1108620/1‑71 (MQ=255)                                                   |                      GTCGGCGACCTGGAACGTATTCTGGCACGTCTGGCTTTACGAACTGCTCGCCCACGCGATCTGGCCCGTATG  >  minE/1827790‑1827861

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑