Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	569,274	A→G	Y414H (TAC→CAC)	hcp ←	hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	569,274	0	A	G	100.0%	99.5 / NA	29	Y414H (TAC→CAC)	hcp	hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/22);  total (7/22)

GAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCC  >  minE/569221‑569304                                                      |                               gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:841752/1‑62 (MQ=255) gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:698780/1‑62 (MQ=255) gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:111990/1‑62 (MQ=255) gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:1190787/1‑62 (MQ=255) gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:606164/1‑62 (MQ=255) gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:440233/1‑62 (MQ=255) gAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCt                        >  1:1018571/1‑62 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:356104/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:919636/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:901644/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:714829/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:641256/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:584068/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:528978/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:505525/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:389572/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:315231/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:291534/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:220673/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:206818/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:1131678/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:1105002/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:1080317/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:1064512/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:1031880/60‑1 (MQ=255)                         cGGCACGCTGGTGGCGAAATCGGTGAAGTGGGGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:572859/60‑1 (MQ=255)                          ggCACGCTGGTGGCGAAATCGGTGAAGTGGTGGCGCGCGCCGCGTGCGCCGTCACAGcc  <  1:1022100/59‑1 (MQ=255)                          ggCACGCTGGGGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:469243/59‑1 (MQ=39)                                cTGGTGGCGAAATCGGTGAAGTGGTGGCGCTCGCCGCGTGCGCCGTCACAGcc  <  1:147683/53‑1 (MQ=255)                                                      |                               GAGGGTCAGGATCAGGCAGTCATCCGGCACGCTGGTGGCGAAATCGGTGAAGTAGTGGCGCTCGCCGCGTGCGCCGTCACAGCC  >  minE/569221‑569304

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