Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,295,693	C→T	A18T (GCG→ACG)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,295,693	0	C	T	100.0%	113.2 / NA	38	A18T (GCG→ACG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/29);  total (9/29)

TTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCA  >  W3110S.gb/2295649‑2295719                                             |                           ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:295081/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:1206080/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:1264029/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:837353/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:2310324/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:396245/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:1538833/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:183418/1‑64 (MQ=255) ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc         >  1:1888285/1‑64 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:505118/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2525176/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:319850/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:32423/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:461336/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:497049/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1075543/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:57291/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:654638/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:708962/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:773112/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:86049/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:988300/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2489630/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2482087/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2300873/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2130609/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2043987/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1899445/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:177237/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1537043/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1510896/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1459466/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1437338/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1400019/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:117129/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:1054199/64‑1 (MQ=255)        tgtgCTTCATCTTTAAACTGCAACACATCGATGGGCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2379504/64‑1 (MQ=255)                                     gATGGTCGTCAGCGCTGAGCCGGAGATcatcagca  <  1:2440418/35‑1 (MQ=255)                                             |                           TTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCA  >  W3110S.gb/2295649‑2295719

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