Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	39974	40049	76	7 [0]	[0] 15	caiA	crotonobetaine reductase subunit II, FAD‑binding

TGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTC  >  W3110S.gb/39912‑39973                                                              | tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:1030978/1‑62 (MQ=255) tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:1624238/1‑62 (MQ=255) tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:394420/1‑62 (MQ=255) tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:409736/1‑62 (MQ=255) tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:506739/1‑62 (MQ=255) tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:846090/1‑62 (MQ=255) tgctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc  >  1:913609/1‑62 (MQ=255)                                                              | TGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTC  >  W3110S.gb/39912‑39973

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