Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,588,674	A→G	E808E (GAA→GAG)	acrD →	aminoglycoside/multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,588,674	0	A	G	100.0%	75.4 / NA	42	E808E (GAA→GAG)	acrD	aminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/0);  total (42/0)

TTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAT  >  W3110S.gb/2588654‑2588710                     |                                     ttCGCGACCTCACGCTGGGAGCCAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:68443/1‑57 (MQ=38) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTgg                   >  1:57428/1‑40 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCtat           >  1:24702/1‑48 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:74281/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:49673/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:51978/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:52571/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:55047/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:57149/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:60191/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:6099/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:71489/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:4855/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:74669/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:7472/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:76129/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:77945/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:81253/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:84068/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:84379/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:9625/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:9727/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:25070/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:12577/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:13396/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:13556/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:14817/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:15556/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:16522/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:1867/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:22304/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:23565/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:11034/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:25593/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:28499/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:3193/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:36930/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:39426/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:40807/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:42012/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:42866/1‑57 (MQ=255) ttCGCGACCTCACGCTGGGAGACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAt  >  1:43714/1‑57 (MQ=255)                     |                                     TTCGCGACCTCACGCTGGGAAACAGGCTCGCCGCGTCTGGAACGCTATAACGGTTAT  >  W3110S.gb/2588654‑2588710

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