Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,276,155	T→C	H159H (CAT→CAC)	yodB →	predicted cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,276,155	0	T	C	100.0%	384.4 / NA	100	H159H (CAT→CAC)	yodB	predicted cytochrome
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/54);  total (46/54)

CTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACAT  >  minE/1276098‑1276163                                                          |         cTGTGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:961966/62‑1 (MQ=255) cTGGTCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1650555/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCTCAGCACTCTTTCACCACt      <  1:2665535/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:3000964/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2125817/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2132282/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2163077/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:222460/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2271496/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2333759/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2488503/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2591896/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2596010/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1000517/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1207106/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2746241/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2764727/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2854252/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2908101/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2675794/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:319389/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:331225/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:3325224/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:3519718/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:3606343/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:3920/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:465252/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:472419/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:500072/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:723822/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:773455/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:843526/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:901534/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:113928/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1807534/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1165704/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1570037/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1738570/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1770469/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1556661/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1804656/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1551895/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:2017078/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1498624/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1346469/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1853073/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1249470/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1911608/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1243623/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1975617/62‑1 (MQ=255) cTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1983098/62‑1 (MQ=255)   gggCAAATATAGGCTACTTTTTAATTGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1703178/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCGCGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3330615/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTa     >  1:3206788/1‑61 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTa     >  1:300721/1‑61 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTa     >  1:161692/1‑61 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2344988/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3515103/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1476785/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3635092/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:380042/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1418587/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:422968/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:459223/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:464090/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1348014/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:134734/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2708431/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:511812/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:596693/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:606800/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1258755/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:769092/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1256507/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2284146/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1939218/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1239796/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2555513/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2767166/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:286365/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2886507/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2697477/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2917993/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2956083/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1785581/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2563903/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3480020/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2467405/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3226776/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3247774/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:1548127/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2460523/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:2387985/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3388671/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACAACTAc    >  1:539688/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTTTAATAGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:310335/1‑62 (MQ=255)   gggCAAATATAGGCTACTTTCTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTAc    >  1:3329368/1‑62 (MQ=255)    ggCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:852947/59‑1 (MQ=255)      cAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACTACAt  <  1:335980/61‑1 (MQ=255)          tataGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCACCACt      <  1:1485135/53‑1 (MQ=255)                                                          |         CTGGGCAAATATAGGCTACTTTTTAATCGCAGCTCACGCTGGCGCAGCACTCTTTCATCACTACAT  >  minE/1276098‑1276163

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