Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,947,335	A→G	100%	L158P (CTG→CCG)	dnaT ←	DNA biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,947,335	0	A	G	100.0%	259.1 / NA	68	L158P (CTG→CCG)	dnaT	DNA biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (68/0);  total (68/0)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAC  >  minE/2947321‑2947391               |                                                         caTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt       >  1:219852/1‑66 (MQ=255) caTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt       >  1:2069457/1‑66 (MQ=255) caTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt       >  1:913412/1‑66 (MQ=255) caTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt       >  1:89232/1‑66 (MQ=255) caTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt       >  1:131986/1‑66 (MQ=255) caTCTCGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgt        >  1:2110701/1‑65 (MQ=255) caTCTCGATTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTttgtt       >  1:2082709/1‑66 (MQ=255)     tcGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:988338/1‑67 (MQ=255)     tcGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:160279/1‑67 (MQ=255)     tcGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:2050430/1‑67 (MQ=255)      cggttCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1937098/4‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGCCCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1289572/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGCCCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:395203/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCt                        >  1:2122796/1‑44 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:293152/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:2103491/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:421167/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:42455/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:292986/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:449996/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:512947/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:555321/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:637321/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:647295/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:66910/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:669362/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:681261/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:720305/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:769804/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:778428/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:792620/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:801140/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:846515/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:849406/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:884838/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:897519/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:944259/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:953800/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:973987/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:156783/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1077367/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1084539/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1086897/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1136025/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1214419/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:129910/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1367599/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1400376/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1403055/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1473956/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1496301/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1516078/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:273296/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1628678/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1649590/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:165890/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1667139/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1699655/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1716460/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1799284/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1809916/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1816502/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1947683/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:2000611/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:20480/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1041708/1‑66 (MQ=255)      cGTTTCGGCGGTCCGCCGTTGCTGGCACGACAGATTTGCAGGCTgc                      >  1:871507/1‑46 (MQ=38)      agTTTCGGCGGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAc  >  1:1127101/2‑66 (MQ=255)               |                                                         CATCTCGTTTCGGCAGTCCGCCGTTGCTGGCACGACCGATTTGCAGGCTGCGCGCCAGTTTTTGTTGCCAC  >  minE/2947321‑2947391

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