Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,647,049	C→T	Q102* (CAA→TAA)	bcp →	thiol peroxidase, thioredoxin‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,647,049	0	C	T	100.0%	45.4 / NA	15	Q102* (CAA→TAA)	bcp	thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/8);  total (7/8)

CGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAT  >  minE/1646990‑1647101                                                            |                                                     cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCGCCAGGTGTGCGAATAATTCGGCGTc                                           <  1:2221740/71‑1 (MQ=255) cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:1778738/71‑1 (MQ=255) cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:2060994/71‑1 (MQ=255) cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:2587296/71‑1 (MQ=255) cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:288702/71‑1 (MQ=255) cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:2991846/71‑1 (MQ=255)          agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:1729039/62‑1 (MQ=255)             tgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTc                                           <  1:1696076/59‑1 (MQ=255)                                 tCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTa           >  1:1858150/1‑71 (MQ=255)                                 tCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTa           >  1:2526359/1‑71 (MQ=255)                                 tCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTa           >  1:2820635/1‑71 (MQ=255)                                 tCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTa           >  1:2840906/1‑71 (MQ=255)                                          gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAt  >  1:1447475/1‑71 (MQ=255)                                          gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAt  >  1:2561006/1‑71 (MQ=255)                                          gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAt  >  1:657334/1‑71 (MQ=255)                                                            |                                                     CGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAT  >  minE/1646990‑1647101

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