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%	46.4 / NA	14	E146* (GAG→TAG)	aroK	shikimate kinase AroK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (8/6);  total (8/6)

ATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGC  >  NZ_CP009273/3511918‑3512065                                                                       |                                                                              aTGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGt                                                            >  2:270247/1‑90 (MQ=255)         aatCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTg                                                               <  1:53747/79‑1 (MQ=255)         aatCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTg                                                               >  2:53747/1‑79 (MQ=255)              ggTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACg                                               <  1:23342/90‑1 (MQ=255)              ggTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACg                                               <  1:84706/90‑1 (MQ=255)              ggTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACg                                               <  2:67843/90‑1 (MQ=255)                        cACTTTAGCGCTTTGATCATCAGTACGAATGGACACGTCGGCAATCTATTCATACAGAGGATTGCGCACATTGGCCAACGCTACCAGCAc                                     >  1:239465/1‑90 (MQ=255)                                cgcTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCAcgcg                             <  2:303019/90‑1 (MQ=255)                                                   aaTGGTCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTgcag          >  1:384180/1‑90 (MQ=255)                                                        tCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGcggcgg                         <  1:239491/70‑1 (MQ=255)                                                        tCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGcggcgg                         >  2:239491/1‑70 (MQ=255)                                                        tCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACg     >  1:372329/1‑90 (MQ=255)                                                        tCACGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACg     >  2:287234/1‑90 (MQ=255)                                                           cGTCGGCAATCTATTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGc  >  1:167974/1‑90 (MQ=255)                                                                       |                                                                              ATGTGAATAATCTGGTTTGCAACCACTTTAGCGCTTTGATCATCAGTACGAATGGTCACGTCGGCAATCTCTTCATACAGCGGATTGCGTTCATTGGCCAACGCTTCCAGAACTTCACGCGGCGGTGTTTCAACGTGCAGCAACGGGC  >  NZ_CP009273/3511918‑3512065

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