Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,298,568	T→C	E49E (GAA→GAG)	ccmE ←	periplasmic heme chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,298,568	0	T	C	100.0%	152.4 / NA	41	E49E (GAA→GAG)	ccmE	periplasmic heme chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/18);  total (23/18)

CAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGAGCGCAG  >  W3110S.gb/2298528‑2298627                                         |                                                            cAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTcc                                     >  1:801189/1‑65 (MQ=255) cAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTcc                                     >  1:1135543/1‑65 (MQ=255) cAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTcc                                     >  1:1397860/1‑65 (MQ=255) cAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTc                                      >  1:1669022/1‑64 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:2055806/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:988701/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:788617/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:510087/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:48859/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:2144377/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1036576/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1887728/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1878121/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1593227/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1037957/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1142582/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1463435/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1493605/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1659020/65‑1 (MQ=255)                   cGACTTCCGGCATTTGCCGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:1690366/65‑1 (MQ=255)                    gACTTCCGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:686660/64‑1 (MQ=255)                         ccGGCATTTGCTGAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAgaga                   <  1:312990/59‑1 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:447977/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:97801/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:958130/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:861861/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1077441/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1096015/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:647876/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1400972/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:464093/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1572322/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:317604/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:309616/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:218948/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:2185441/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:2154637/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1723643/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:209468/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1814922/1‑64 (MQ=255)                                     gAGTCTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGagcgcag  >  1:1531913/1‑64 (MQ=255)                                         |                                                            CAACGCGCAGACGCTGACCGACTTCCGGCATTTGCTGAGTTTCACGCTTGCCGTAGAGAATTTCCCCCGGCGTATAAAAGAGATCGATATTCGAGCGCAG  >  W3110S.gb/2298528‑2298627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑