Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4583009	4583258	250	2 [0]	[0] 43	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

CGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAT  >  W3110S.gb/4583259‑4583322 |                                                                cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGtat                            >  1:1549159/1‑38 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaaa                        >  1:196890/1‑42 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaa                          >  1:790026/1‑40 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCa                  >  1:836661/1‑48 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGa   >  1:700756/1‑63 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:966833/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:995298/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:405063/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:408781/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:476293/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:476552/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:610050/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:614276/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:657288/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:873772/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:675789/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1864609/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:719728/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:726409/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:940440/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:893677/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:855224/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1390636/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1104967/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1120212/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1158821/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:125540/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1293917/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1296327/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1326686/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1337334/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1342212/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:193871/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1467824/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1568856/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1589327/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1592248/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1712645/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1769140/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1820113/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1837864/1‑64 (MQ=255) cGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1073795/1‑64 (MQ=255) cGTCAGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAt  >  1:1201176/1‑64 (MQ=255) |                                                                CGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGAT  >  W3110S.gb/4583259‑4583322

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