Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	2588947	2588999	53	24 [0]	[0] 10	acrD	aminoglycoside/multidrug efflux system

GCGGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGC  >  W3110S.gb/2589000‑2589060 |                                                             gcgGGCTGGAAAACGACGTTTACTTCCAGGTGTGCCTGTTAACGGTCATTGGTTTATCgg   >  1:1333296/1‑60 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCgg   >  1:112694/1‑60 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:259062/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:314014/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:327246/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:371695/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:406574/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:582309/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:903569/1‑61 (MQ=255) gcgGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGc  >  1:952277/1‑61 (MQ=255) |                                                             GCGGGCTGGAAAACGACGTTTACTTCCAGGTGGGCCTGTTAACGGTCATTGGTTTATCGGC  >  W3110S.gb/2589000‑2589060

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