Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4585777	4585782	6	2 [0]	[0] 13	hsdS	specificity determinant for hsdM and hsdR

AAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATT  >  W3110S.gb/4585783‑4585843 |                                                             aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATgc                >  1:2003963/1‑47 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:1017449/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:1175609/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:1347521/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:1662295/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:2013977/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:2162087/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:2265768/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:2533471/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:2777356/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:295634/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:403782/1‑61 (MQ=255) aaaCCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACAtt  >  1:608565/1‑61 (MQ=255) |                                                             AAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATT  >  W3110S.gb/4585783‑4585843

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: ‑