Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4585699	4585709	11	7 [0]	[0] 5	hsdS	specificity determinant for hsdM and hsdR

AGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCC  >  W3110S.gb/4585638‑4585698                                                             | aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:1034951/61‑1 (MQ=255) aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:1346605/61‑1 (MQ=255) aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:189229/61‑1 (MQ=255) aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:2451149/61‑1 (MQ=255) aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:2754390/61‑1 (MQ=255) aGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:336972/61‑1 (MQ=255)  gCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGcc  <  1:1399858/60‑1 (MQ=255)                                                             | AGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCC  >  W3110S.gb/4585638‑4585698

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