Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	minE	1039306	1039328	23	5 [0]	[0] 35	nth	DNA glycosylase and apyrimidinic (AP) lyase

TCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGA  >  minE/1039246‑1039305                                                            | tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa  <  1:1659201/60‑1 (MQ=255) tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa  <  1:2140029/60‑1 (MQ=255) tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa  <  1:2526972/60‑1 (MQ=255) tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGa  <  1:2708438/60‑1 (MQ=255) tCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTAGCCTTTTGa  <  1:2900312/60‑1 (MQ=255)                                                            | TCGCCTGCGTGAGAACAATCCTCATCCCACCACCGAGCTTAATTTCAGTTCGCCTTTTGA  >  minE/1039246‑1039305

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