Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1407487	1407511	25	2 [1]	[1] 3	dgcM	diguanylate cyclase DgcM

TCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGA  >  NC_000913/1407506‑1407639       |                                                                                                                                tctgtCAGCTTTCGCAAAACAAGCCGACAGCTTAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCtttt          >  1:585067/1‑126 (MQ=255)       aGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTGGATAAAACCGTTTTGGTTTTGCGTAT‑GGGGAGGGATTTTGGGAGGTTTTTTtttggtgg  >  1:226532/1‑123 (MQ=255)       aGCTTTCGCAAAAAAAGCCGATAGCCTAATGCTGTTTCCTCTTTTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTgagga  <  1:602114/128‑1 (MQ=255)       |                                                                                                                                TCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGA  >  NC_000913/1407506‑1407639

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