Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	minE	1702921	1702971	51	16 [0]	[0] 24	hcaE	3‑phenylpropionate dioxygenase, large (alpha) subunit

ACTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCT  >  minE/1702972‑1703026 |                                                       aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGc   >  1:594984/1‑54 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGc   >  1:228767/1‑54 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1969304/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:848301/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:439374/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:337194/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:3279396/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:3266297/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:2655363/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:2642549/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:246794/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:2254099/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1086669/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1759475/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1707638/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1655305/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1610680/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1554968/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1357632/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1204173/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1133264/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1132765/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1132539/1‑55 (MQ=255) aCTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCt  >  1:1118807/1‑55 (MQ=255) |                                                       ACTGCCGCTCGTGGAATGTACCAACGCTGGGCCGATCTTCTGAGTAGCGAAAGCT  >  minE/1702972‑1703026

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