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%	166.4 / NA	51	A18T (GCG→ACG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/28);  total (23/28)

AGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCAT  >  W3110S.gb/2295631‑2295714                                                               |                      aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:435937/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:2583628/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:2474318/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:264736/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:2803931/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:2147894/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:2879373/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:1973227/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:1562358/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:1191175/1‑65 (MQ=255) aGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:942613/1‑65 (MQ=255) aGTGAGCAGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                     >  1:82534/1‑65 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGct           <  1:318058/64‑2 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:1523072/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:1092138/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:999162/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:943148/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:2807229/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:1431323/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:3088558/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:888897/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:376827/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:1748962/64‑1 (MQ=255)            ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGcc           <  1:1575861/64‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:1190087/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:2659859/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:438987/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:503215/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:2375335/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:1261734/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:1261914/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:1612216/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:175371/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:1899684/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:2250862/65‑1 (MQ=255)               aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:2344457/65‑1 (MQ=255)               aCTGTTGGTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:636593/65‑1 (MQ=255)                cTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:2129972/64‑1 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:532714/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:2466633/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:962291/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:2726995/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:1329976/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:1526617/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:2730144/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:1640147/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:1711742/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatc    >  1:1717042/1‑64 (MQ=255)                   ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGCTCGTCAGCGCTGAGCCGGAGata    >  1:2001367/1‑63 (MQ=255)                     gttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGatcat  <  1:500847/64‑1 (MQ=255)                        ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGgag       <  1:659094/56‑1 (MQ=255)                                                               |                      AGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCAT  >  W3110S.gb/2295631‑2295714

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