Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4582364	4582496	133	43 [0]	[0] 30	mcrC	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

TAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGA  >  W3110S.gb/4582302‑4582363                                                              | tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:681850/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1006688/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:344467/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:369391/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:38136/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:446652/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:455353/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:504403/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:515742/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:524920/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:533618/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:2928/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:726647/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:744745/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:77804/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:817290/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:872272/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:877377/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:916759/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:973207/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:999680/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1459244/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1123543/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1124002/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1134159/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1187403/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1194663/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1324678/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1380661/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1453397/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:272646/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1521063/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1547683/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1570110/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1576768/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1591690/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:162851/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:1632239/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:169788/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:226192/62‑1 (MQ=255) tAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:239452/62‑1 (MQ=255)   aTGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:886302/60‑1 (MQ=255)                           ccAGCGTGTCTTCATTAAGCATATCAAAAGTACTGa  <  1:822490/36‑1 (MQ=255)                                                              | TAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGA  >  W3110S.gb/4582302‑4582363

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