Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	1408663	1408664	2	27 [0]	[0] 28	ydaM	predicted diguanylate cyclase, GGDEF domain signalling protein

GCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAA  >  W3110S.gb/1408665‑1408729 |                                                                 gCGATGGCTAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1227308/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTTGCCGATTCAGTaa  >  1:654828/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:222216/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:948873/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:815061/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:714393/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:641511/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:608279/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:561130/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:414870/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:311758/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:2513929/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:2479642/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:2437909/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:2428074/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:2205034/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:2064157/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1968982/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1961276/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1650205/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1581383/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1486231/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1434456/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1423193/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1382478/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:135075/1‑65 (MQ=255) gCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTa   >  1:1549407/1‑64 (MQ=255) gCGATGCCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTaa  >  1:1971619/1‑65 (MQ=255) |                                                                 GCGATGGCGAGATGCTGCATTTGGCCTGGTTCCGTAATGTGATAAAACTGTCGCCGATTCAGTAA  >  W3110S.gb/1408665‑1408729

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