Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4586752	4586757	6	4 [0]	[0] 76	hsdM	DNA methylase M

GAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGA  >  W3110S.gb/4586758‑4586819 |                                                              gaATGTCGTTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:885721/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:356744/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:559024/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:552666/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:552608/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:494229/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:469268/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:405690/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:396565/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:384372/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:370720/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:561250/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:345026/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:326900/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:321320/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:313172/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:312156/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:305897/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:291608/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:276983/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:863004/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:997243/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:950581/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:943723/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:923083/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:918252/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:908810/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:895606/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:892110/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:275483/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:845752/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:828726/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:803744/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:801814/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:782044/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:749802/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:676136/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:60718/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1191566/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1359693/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1332284/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1330016/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1274752/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1252592/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1252337/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1222680/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1222452/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1215574/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1361477/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1179756/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1178030/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1148532/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1141143/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1093332/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1073338/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1043583/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:104093/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:157180/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:275120/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:247967/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:243250/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:234455/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:22443/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:192147/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:190188/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:183133/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:18158/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1025809/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:147724/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1465362/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1463272/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1444639/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1409407/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1407422/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:137859/62‑1 (MQ=255) gaATGTCGGTGCCTTTGCCGCCATCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGa  <  1:1406347/62‑1 (MQ=255) |                                                              GAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGA  >  W3110S.gb/4586758‑4586819

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