Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4588361	4588377	17	19 [0]	[0] 35	hsdR	endonuclease R

GGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATA  >  W3110S.gb/4588296‑4588360                                                                 | ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1726502/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:995670/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:765874/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:2275854/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:2275332/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:2147858/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1962981/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1908176/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1844560/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1260518/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1708582/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1696324/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1682629/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1672297/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1668793/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1626208/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1594000/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATtaata  >  1:1399193/1‑65 (MQ=255) ggCGGTCAAACCACTCCTGTAGCACGACCAGCCCTTTACGGGTGAGATAGCGCGGGCGATtaata  >  1:537581/1‑65 (MQ=255)                                                                 | GGCGGTCAAACCACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATA  >  W3110S.gb/4588296‑4588360

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