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

GTTTTGCGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCT  >  minE/1646984‑1647115                                                                  |                                                                   gTTTTGCGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTc                                                               <  1:892884/71‑1 (MQ=255) gTTTTGCGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTc                                                               >  1:2657375/1‑71 (MQ=255)    ttGCGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGc                                                            >  1:1903205/1‑71 (MQ=255)       cGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGt                                                          <  1:1078555/70‑1 (MQ=255)         gAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTg                                                       <  1:2626155/71‑1 (MQ=255)            aaaGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTgggg                                                    >  1:1855256/1‑71 (MQ=255)              agagCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGt                                                   >  1:2331874/1‑70 (MQ=255)                agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGaa                                                <  1:241193/71‑1 (MQ=255)                agCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGa                                                 >  1:2672433/1‑70 (MQ=255)                     ctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCTAATAATTCGGCGTCTGGGGTGAAAAg                                             <  1:2479346/69‑1 (MQ=255)                     ctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAg                                             <  1:37350/69‑1 (MQ=255)                     ctTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGt                                            >  1:1145955/1‑70 (MQ=255)                        aaCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGt                                                   >  1:733981/1‑60 (MQ=255)                        aaCTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTcc                                          <  1:537659/69‑1 (MQ=255)                             tACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCAt                                     <  1:89684/69‑1 (MQ=255)                             tACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATgg                                   >  1:2064129/1‑71 (MQ=255)                             tACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATg                                    <  1:1181508/70‑1 (MQ=255)                                gCTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTg                                                  <  1:1725937/53‑1 (MQ=255)                                 cTCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCaa                               <  1:1425411/71‑1 (MQ=255)                                  tCCTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCaa                               >  1:1137242/1‑70 (MQ=255)                                   ccTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCaaaa                             >  1:3019818/1‑71 (MQ=255)                                    cTGTCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCaaa                              <  1:1101752/69‑1 (MQ=255)                                       tCTGATGAGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCaa                               <  1:27570/65‑1 (MQ=255)                                            tgaGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATg                    >  1:1075964/1‑71 (MQ=255)                                             gaGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                   >  1:743449/1‑71 (MQ=255)                                              aGGACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATgg                   <  1:2635074/70‑1 (MQ=255)                                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAt                <  1:1751890/71‑1 (MQ=255)                                                gACCACCAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAt                >  1:380895/1‑71 (MQ=255)                                                 accaccAGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCAtt               >  1:2270130/1‑71 (MQ=255)                                                     ccaGGTGTGCGAATAATACGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATc           >  1:1233984/1‑71 (MQ=255)                                                       aGGTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCg          >  1:1670893/1‑70 (MQ=255)                                                        ggTGTGCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCa        <  1:1116017/71‑1 (MQ=255)                                                            tgCGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCa     <  1:1372171/70‑1 (MQ=255)                                                              cGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCt                          >  1:333982/1‑47 (MQ=255)                                                              cGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTAc                        >  1:1850844/1‑49 (MQ=255)                                                              cGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTAc                        >  1:1563854/1‑49 (MQ=255)                                                              cGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGc   <  1:2251254/70‑1 (MQ=255)                                                              cGAATAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCt  >  1:1043142/1‑71 (MQ=255)                                                                  |                                                                   GTTTTGCGGAAAAAGAGCTGCTTAACTTTACGCTCCTGTCTGATGAGGACCACCAGGTGTGCGAACAATTCGGCGTCTGGGGTGAAAAGTCCTTCATGGGCAAAACCTACGATGGCATTCATCGCATCAGCT  >  minE/1646984‑1647115

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