Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4221317	4221361	45	22 [0]	[0] 14	aceA	isocitrate lyase

ACGTGACGGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATC  >  W3110S.gb/4221362‑4221423 |                                                              acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGAt        >  1:159130/1‑56 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:116359/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:1404804/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:1458320/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:1559168/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:181494/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:1848190/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:2410114/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:2476383/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:2873269/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:29332/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:964971/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATc  >  1:966334/1‑62 (MQ=255) acgtgacgGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCAGATc  >  1:2571861/1‑62 (MQ=255) |                                                              ACGTGACGGGCGTTCCAACCCTGCTGGTTGCCCGTACCGATGCTGATGCGGCGGATCTGATC  >  W3110S.gb/4221362‑4221423

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