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%	223.7 / NA	67	Q102* (CAA→TAA)	bcp	thiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/38);  total (29/38)

AGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAG  >  minE/1646999‑1647113                                                   |                                                                 agCTGCTTAACTTTACGCTCCTGTCTGATGAGGAGCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1807451/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:266487/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3851158/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3850453/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3733714/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3709716/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:370007/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3689639/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:3655974/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2964720/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:4309843/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2684316/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2234873/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:2120425/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:756449/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:115157/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1404284/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:182079/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1519543/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1711092/65‑1 (MQ=255) agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:1696761/65‑1 (MQ=255)  gctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgg                                                    <  1:4080540/64‑1 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:1968860/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:3723766/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:1566976/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:4077779/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:3476754/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:3443179/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:1654293/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:4344681/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                  >  1:582602/1‑65 (MQ=255)   ctgctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTAGGCGTCTgggg                                                  >  1:1418398/1‑65 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:4088827/70‑1 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:4323257/70‑1 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:1619927/70‑1 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:2383450/70‑1 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:812331/70‑1 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:140293/70‑1 (MQ=255)          aCTTTACGCTCCTGTCTGATGAGGACCACCAGGGGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCtt                                      <  1:2900827/70‑1 (MQ=255)                     cTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAAc                          <  1:4234048/71‑1 (MQ=255)                     cTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAAc                          <  1:3437399/71‑1 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:812615/1‑68 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:576513/1‑68 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:3947111/1‑68 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:2775967/1‑68 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:1239491/1‑68 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                 >  1:1829055/1‑68 (MQ=255)                                 gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACCATgg                 >  1:2991157/1‑68 (MQ=255)                                    caccaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAAc                          >  1:2911917/1‑56 (MQ=255)                                     accaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCAt          <  1:3956866/71‑1 (MQ=255)                                     accaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCAt          <  1:1897528/71‑1 (MQ=255)                                     accaGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCAt          <  1:2433377/71‑1 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3999313/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:653749/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:641857/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:4377201/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:2548188/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3924825/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:4283407/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3525863/1‑65 (MQ=255)                                        aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt             >  1:3985278/1‑65 (MQ=255)                                             tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:455203/71‑1 (MQ=255)                                             tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:204531/71‑1 (MQ=255)                                             tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  >  1:3077129/1‑71 (MQ=255)                                             tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  >  1:720182/1‑71 (MQ=255)                                             tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:378254/71‑1 (MQ=255)                                              gCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAg  <  1:3139324/70‑1 (MQ=255)                                                   |                                                                 AGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAG  >  minE/1646999‑1647113

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