Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	1408894	1408916	23	2 [0]	[0] 3	ydaM	predicted diguanylate cyclase, GGDEF domain signalling protein

CCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTC  >  W3110S.gb/1408917‑1408978 |                                                              ccGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTc  <  1:1005457/62‑1 (MQ=255) ccGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTc  <  1:1149938/62‑1 (MQ=255) ccGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTc  <  1:842910/62‑1 (MQ=255) |                                                              CCGGGTAAATGCGAGATTTCATGCATGATAGGCATGACGCGACGCCCGAGCATATTTATTTC  >  W3110S.gb/1408917‑1408978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑