Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4589352	4589426	75	73 [0]	[0] 22	hsdR	endonuclease R

ATAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAA  >  W3110S.gb/4589290‑4589351                                                              | aTAATTTGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1635331/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCGGGGATCACCAGGCCACGGTTaa  <  1:933438/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:502398/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:707359/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:647000/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:598767/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:579762/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:573851/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:566463/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:546675/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:5275/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:754218/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:4906/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:456116/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:443484/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:43709/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:417458/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:364112/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:355736/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:860838/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:984729/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:980562/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:973729/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:954815/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:907591/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:90220/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:899603/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:873453/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:324394/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:853612/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:82153/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:814967/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:799590/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:790477/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:784986/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:76303/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1280181/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1515365/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1458435/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1429639/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1390585/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1353803/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1343246/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1317755/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1316412/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1522367/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1277631/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1255416/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1252909/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:121999/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1188643/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1179557/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1133403/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1125076/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1632071/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:31150/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:30110/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:296592/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:283254/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:267861/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:253692/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:21925/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:202398/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:174850/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1103598/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1598651/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1559772/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1556299/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1542659/62‑1 (MQ=255) aTAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:1525357/62‑1 (MQ=255) aTAATTGATGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:553242/62‑1 (MQ=255)    aTTGGTGAGTTCTTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:557824/59‑1 (MQ=255)                gTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTaa  <  1:142493/47‑1 (MQ=255)                                                              | ATAATTGGTGAGTTCGTTACAGACGGCGCGGTTAAACGCCGGGATCACCAGGCCACGGTTAA  >  W3110S.gb/4589290‑4589351

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