Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	minE	2408075	2408079	5	17 [0]	[0] 76	hemC	hydroxymethylbilane synthase

GAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGAT  >  minE/2408039‑2408074                                    | gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:2541234/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:3484631/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:3436417/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:3368816/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:324150/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:2994287/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:2970525/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:2915810/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:2819607/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:1107730/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:244768/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:2407458/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:1900113/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:1788618/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:1724689/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:1655152/1‑36 (MQ=255) gagCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGat  >  1:1604744/1‑36 (MQ=255)                                    | GAGCTGGAAGTCGCGCTCCTCGAAAATCGCGCCGAT  >  minE/2408039‑2408074

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: ‑