Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	2,970,677	(T)5→6	coding (403/645 nt)	ytjB ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,970,672	1	.	T	100.0%	39.1 / NA	13	N136K (AAC→AAA)	ytjB	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (13/0);  total (13/0)

CTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑TTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCC  >  minE/2970617‑2970719                                                         |                                                cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGTTAGCCGAGCGGTCCG‑ttttt                                            <  1:446694/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:752820/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:366600/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:368883/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:377220/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:423009/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:4802/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:514925/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:51514/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:656692/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:65772/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:702310/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:454682/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:777026/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:78028/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:796680/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:911377/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:918083/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:960/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:976814/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:980976/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:989369/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:994004/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1180696/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:101349/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1014780/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1042094/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1061166/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1064220/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1071749/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1080246/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1107536/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1155593/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1171972/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:1190215/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:169731/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:302147/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:244401/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:175288/61‑1 (MQ=255) cTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:207244/61‑1 (MQ=255)  ttCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:785532/60‑1 (MQ=255)  ttCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:271664/60‑1 (MQ=255)  ttCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:463323/60‑1 (MQ=255)  ttCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:775964/60‑1 (MQ=255)  ttCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCAGAGCGGTCCG‑ttttt                                            <  1:1078737/60‑1 (MQ=255)   tCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:712942/59‑1 (MQ=255)      gTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGTTAGCCGAGCGGTCCG‑ttttt                                            <  1:921011/56‑1 (MQ=255)            aGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑ttttt                                            <  1:980333/50‑1 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:465628/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:131793/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:1009977/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:891091/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:553256/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:26288/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:648552/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgcc  >  1:404068/1‑62 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgc   >  1:548736/1‑61 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgc   >  1:478096/1‑61 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgc   >  1:947829/1‑61 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgc   >  1:21707/1‑61 (MQ=255)                                           aGCCGAGCGGTCCGTTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAgccgc   >  1:339422/1‑61 (MQ=255)                                                         |                                                CTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCG‑TTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCC  >  minE/2970617‑2970719

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