Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4590510	4590822	313	15 [0]	[0] 62	hsdR	endonuclease R

GAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCT  >  W3110S.gb/4590472‑4590509                                      | gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1029163/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1146528/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1307240/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1369299/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1423040/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1542397/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1666287/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1676967/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1699483/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:432210/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:496115/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:567714/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:634838/38‑1 (MQ=255) gAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:713924/38‑1 (MQ=255)  aGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct  <  1:1732266/37‑1 (MQ=255)                                      | GAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCT  >  W3110S.gb/4590472‑4590509

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