Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4585893	4585908	16	11 [0]	[0] 25	hsdS	specificity determinant for hsdM and hsdR

ACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACC  >  W3110S.gb/4585909‑4585968 |                                                            aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1902592/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:937389/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:809759/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:668189/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:468499/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:332123/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2597384/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2590037/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2576548/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2369524/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2350923/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2331398/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:2029436/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1220360/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1766297/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1738437/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1687854/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1624820/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1616945/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1564538/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1544013/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:147082/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1445463/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1406933/60‑1 (MQ=255) aCAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAAcc  <  1:1267869/60‑1 (MQ=255) |                                                            ACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGGAACAAAAACC  >  W3110S.gb/4585909‑4585968

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