Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,745,614	(A)5→6	coding (791/1176 nt)	mdtL ←	multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,745,609	1	.	A	100.0%	38.5 / NA	16	K266* (AAG→TAG)	mdtL	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (0/16);  total (0/16)

TGCGAGGTGATCATCAACGTACGTGGCTT‑AAAAATTCCCAGCGCAAATGGCGTGGAGAATG  >  W3110S.gb/3745581‑3745641                              |                                 tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:1210146/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:1252394/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:1410045/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:181013/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:37958/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:560211/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:604088/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:613611/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:616739/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:633587/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:75897/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:782135/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:834623/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:91064/62‑1 (MQ=255) tGCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:940498/62‑1 (MQ=255)  gCGAGGTGATCATCAACGTACGTGGCTTAAAAAATTCCCAGCGCAAATGGCGTGGAGAATg  <  1:450921/61‑1 (MQ=255)                              |                                 TGCGAGGTGATCATCAACGTACGTGGCTT‑AAAAATTCCCAGCGCAAATGGCGTGGAGAATG  >  W3110S.gb/3745581‑3745641

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