Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,181,222	C→G	R192P (CGC→CCC)	thiH ←	2‑iminoacetate synthase ThiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,181,222	0	C	G	100.0%	26.6 / NA	10	R192P (CGC→CCC)	thiH	2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (5/5);  total (5/5)

ATCAATCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGCGGGCATAAGTCGCCTCGTGATATGTCTCCTGATAAACCATCACGCCATCCAGACCAAGTTGCTTTAACTCG  >  NZ_CP009273/4181143‑4181293                                                                                |                                                                        aTCAATCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGt                                                               <  2:64032/90‑1 (MQ=255) aTCAACCCCCCCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTG‑TTTCTGCCTTTCAGATGATGGGGGGCATAAGTc                                                              <  2:323342/90‑1 (MQ=255)      tCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCt                                                          >  1:342305/1‑90 (MQ=255)      tCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCt                                                          >  1:402862/1‑90 (MQ=255)      tCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCt                                                          >  1:435682/1‑90 (MQ=255)      tCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCt                                                          >  2:312183/1‑90 (MQ=255)                              gCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCTCGTGATATGTCTCCTGATAAACCa                                  >  1:279330/1‑90 (MQ=255)                                    ccAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCTCGTGATATGTCTCCTGATAAACCATCACGc                            <  1:160294/90‑1 (MQ=255)                                    ccAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGGGGGCATAAGTCGCCTCGTGATATGTCTCCTGATAAACCATCACGc                            <  2:460534/90‑1 (MQ=255)                                                              ttGCCTTTCAGATGATGGGGGGCATAAGTCGCCTCGTGATATGTCTCCTGATAACCCATCACGCCATCCAGACCAAGTTGCTTTAACTCg  <  2:228338/90‑1 (MQ=255)                                                                                |                                                                        ATCAATCCCCGCACGCCCCAGCCGATCCGGCGTTTCCAGCCGCCAGAAGAAGTCCTGTTTTTTGCCTTTCAGATGATGGCGGGCATAAGTCGCCTCGTGATATGTCTCCTGATAAACCATCACGCCATCCAGACCAAGTTGCTTTAACTCG  >  NZ_CP009273/4181143‑4181293

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