Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4589285	4589289	5	10 [0]	[0] 71	hsdR	endonuclease R

CCACCATATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGG  >  W3110S.gb/4589223‑4589284                                                              | ccaccaTATCGGCATGGGCATTGGTGACGCTGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:142219/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:1050529/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:1177730/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:1373433/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:1376373/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:485455/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:524317/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:62154/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:697765/62‑1 (MQ=255) ccaccaTATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCggg  <  1:881582/62‑1 (MQ=255)                                                              | CCACCATATCGGCATGGGCATTGGTGACGCAGAAGACCAGCGTTTTTTGCGATCCGGTCGGG  >  W3110S.gb/4589223‑4589284

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