Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	CP000730	459969	460046	78	2 [1]	[0] 2	hsdS1	type I site‑specific deoxyribonuclease specificity subunit

GGTGATTTCTTTAAAAAAATGGATATATTGATAAGTAAACAGAAAATGAAAATTGAAATATTAGAAAAAGAGAAACAATCCTTTTTACAAAAAATGTTCTTATAACTTTGATAAATACATAGATTGCATAAGAATAAAATTTGTATAATTT  >  CP000730/460047‑460197 |                                                                                                                                                       ggtggtTTCTTTAAAAAAATGGGTATATTGATAAGTAAACAGAAAATGAAAATTGAAATATTAGAAAAAGAGAAACAATCCTTTTTACAAAAAATGTTCTTATAACTTTGATAAATACATAGATTGCATAAGAATAAAATTTGTATAAttt  >  2:215004/1‑151 (MQ=255) ggTGATTTCTTTAAAAAAATGGATATATTGATAAGTAAACAGAAAATGAAAATTGAAATATTAGAAAAAGAGAAACAATCCTTTTTACAAAAAATGTTCTTATAACTTTGATAAATACATAGATTGCATAAGAATAAAATTTGTATAAttt  >  1:104409/1‑151 (MQ=255) |                                                                                                                                                       GGTGATTTCTTTAAAAAAATGGATATATTGATAAGTAAACAGAAAATGAAAATTGAAATATTAGAAAAAGAGAAACAATCCTTTTTACAAAAAATGTTCTTATAACTTTGATAAATACATAGATTGCATAAGAATAAAATTTGTATAATTT  >  CP000730/460047‑460197

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