Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,511,988	C→A	E146* (GAG→TAG)	aroK ←	shikimate kinase AroK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,511,988	0	C	A	100.0%	24.3 / NA	8	E146* (GAG→TAG)	aroK	shikimate kinase AroK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (6/2);  total (6/2)

TTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAA  >  NZ_CP009273/3511911‑3512060                                                                              |                                                                         ttCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCgg                                                              >  2:74715/1‑90 (MQ=255)       cATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCg                                                        <  2:370050/90‑1 (MQ=255)                    tGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAAc                                           >  2:66862/1‑90 (MQ=255)                      gTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGc                                         >  2:114227/1‑90 (MQ=255)                                           ttGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGcggcgg                    >  2:39514/1‑90 (MQ=255)                                           ttGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGcggcgg                    >  2:58904/1‑90 (MQ=255)                                                   caGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCaa            <  2:343186/90‑1 (MQ=255)                                                             ggTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCaa  >  1:380941/1‑90 (MQ=255)                                                                              |                                                                         TTCCAGCATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAA  >  NZ_CP009273/3511911‑3512060

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