Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4589492	4589496	5	7 [0]	[1] 5	hsdR	endonuclease R

GCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGA  >  W3110S.gb/4589421‑4589491                                                                       | gcTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGTGGa  <  1:834397/71‑1 (MQ=255) gcTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGa  <  1:1188779/71‑1 (MQ=255) gcTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGa  <  1:1419381/71‑1 (MQ=255) gcTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGa  <  1:653786/71‑1 (MQ=255) gcTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGa  <  1:898512/71‑1 (MQ=255)  ctctACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGa  <  1:2663846/70‑1 (MQ=255)  ctctACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGa  <  1:3003234/70‑1 (MQ=255)                                                                       | GCTCTACCTGCTCGCCTTTGGAGAGATAAACCCCCTCCTGCGCGTTGCGGGTGATGATCTGAATAGGCGGA  >  W3110S.gb/4589421‑4589491

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