Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4585527	4585585	59	4 [0]	[0] 6	hsdS	specificity determinant for hsdM and hsdR

CAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTT  >  W3110S.gb/4585586‑4585650 |                                                                 cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGAttttt  >  1:1528245/1‑65 (MQ=255) cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGAttttt  >  1:40435/1‑65 (MQ=255) cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGAttttt  >  1:437177/1‑65 (MQ=255) cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGAttttt  >  1:898336/1‑65 (MQ=255) cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGAtttt   >  1:1087919/1‑64 (MQ=255) cAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTGTTCAGCGATGAttttt  >  1:700651/1‑65 (MQ=255) |                                                                 CAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTT  >  W3110S.gb/4585586‑4585650

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