Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4587009	4587062	54	15 [0]	[0] 29	hsdM	DNA methylase M

CGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGT  >  W3110S.gb/4587063‑4587124 |                                                              cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:2364953/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:773542/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:761033/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:686817/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:626172/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:599174/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:525576/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:497474/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:393540/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:312614/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:2888401/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:277895/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:270229/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:2506749/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1166657/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:2361881/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:20915/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1951257/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1924594/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1868297/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1789182/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1712400/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1691457/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1686327/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1584297/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1396617/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1353889/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:1294291/62‑1 (MQ=255) cgacgGGTGCCGGGCACCAGTTCGACGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgt  <  1:2636871/62‑1 (MQ=255) |                                                              CGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGT  >  W3110S.gb/4587063‑4587124

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