Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4585242	4585614	373	14 [0]	[0] 48	hsdS	specificity determinant for hsdM and hsdR

GAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAAC  >  W3110S.gb/4585181‑4585241                                                             | gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:1090295/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:1090561/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:1292092/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:1308071/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:1438246/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:206038/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:293648/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:441733/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:511188/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:714623/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:802946/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAac  >  1:953198/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCAACAAac  >  1:69118/1‑61 (MQ=255) gAATAAGTTTATCAGGATATAGCAAATTTTCATGTTGTAATTTTTTCAATAACCCACAAac  >  1:364579/1‑61 (MQ=255)                                                             | GAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAAC  >  W3110S.gb/4585181‑4585241

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