Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4583191	4583213	23	25 [0]	[0] 27	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

TCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAA  >  W3110S.gb/4583127‑4583190                                                                | tacATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:199562/3‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:2096255/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:977743/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:882316/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:659657/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:552600/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:548438/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:471730/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:415163/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:32454/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:306909/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:298760/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:2108640/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1130378/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:2049169/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:2007927/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1881700/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1779744/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1505499/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1503226/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1463572/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1333205/1‑64 (MQ=255) tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1274119/1‑64 (MQ=255) tCCATTAACATCATCACTTCCCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:1781956/1‑64 (MQ=255) tCCATTAACATAATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTataa  >  1:562375/1‑64 (MQ=255)                                                                | TCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAA  >  W3110S.gb/4583127‑4583190

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