Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4591473	4591476	4	22 [0]	[0] 6	hsdR	endonuclease R

AATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCA  >  W3110S.gb/4591477‑4591541 |                                                                 aaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTca  <  1:1301100/65‑1 (MQ=255) aaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTca  <  1:2237468/65‑1 (MQ=255) aaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTca  <  1:2619023/65‑1 (MQ=255) aaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTca  <  1:297162/65‑1 (MQ=255) aaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTca  <  1:509733/65‑1 (MQ=255) aaTATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTca  <  1:55305/65‑1 (MQ=255) |                                                                 AATATTTAAGGCCCATAACATCTCATCTTAGCTTTCTGTACCTTTCCGGGCAATGACCACGGTCA  >  W3110S.gb/4591477‑4591541

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