Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1623602	1623748	147	2 [1]	[1] 2	dgcZ	diguanylate cyclase DgcZ

TATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCTTTTTCACAGT  >  NC_000913/1623722‑1623870                            |                                                                                                                          tATCTCTGGATTACTGGTCTCGCTCGCGACCACGCTGTGAACCTTACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCa                       <  1:599942/128‑1 (MQ=255)                      tctccggACCCCGCTGTGAAAACTTCCCCCCAGGCACTGGTAGTTGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCTTTTTCACAGt  <  1:234654/122‑1 (MQ=255)                            |                                                                                                                          TATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATTTAACAAGATGGCATCAATTTCCGTTGTCTTCTTGATCATTGCCACTCCTTTTTCACAGT  >  NC_000913/1623722‑1623870

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