Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4581792	4581794	3	6 [0]	[0] 42	mcrC	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

CAGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTC  >  W3110S.gb/4581795‑4581855 |                                                             caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:794406/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:257393/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:266593/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:282327/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:303406/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:337427/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:353361/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:67421/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:731446/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:760560/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:766241/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:198773/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:80170/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:813737/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:890178/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:905751/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:908765/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:918761/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:919681/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:935680/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:942336/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1248297/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:103240/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1037697/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1079983/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1116057/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1124476/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1132631/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1154017/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1163662/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1185628/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1212683/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1022961/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1252041/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1302570/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1333063/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1390649/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1392984/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1404394/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1416873/61‑1 (MQ=255) caGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:152494/61‑1 (MQ=255) caGTTGATAAAGATTTTGAGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGctc  <  1:1360592/61‑1 (MQ=255) |                                                             CAGTTGATAAAGATTTTGCGAATGAAATTTTTCTGTTCCCATTCGTCGTGAAAAAATGCTC  >  W3110S.gb/4581795‑4581855

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