Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	2146545	2146563	19	31 [0]	[0] 6	yegE	predicted diguanylate cyclase, GGDEF domain signalling protein

AACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCTACTA  >  W3110S.gb/2146484‑2146544                                                             | aaCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:551207/61‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1037650/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:941449/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:932237/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:8186/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:7526/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:570523/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:479319/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:438500/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:382543/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:372513/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:352867/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2913924/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2896640/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2853182/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2800540/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2509215/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:243872/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2385027/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1801577/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1658963/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1581926/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1453669/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1430289/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1430186/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1290516/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1161326/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1059325/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:1004835/60‑1 (MQ=255)  aCAGGTTGAAAAGCTGATAAGCCGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:261142/60‑1 (MQ=255)                       cGGTGAAATAAACACCTATTCAATGGAAAAACGctacta  <  1:2025360/39‑1 (MQ=255)                                                             | AACAGGTTGAAAAGCTGATAAGCGGTGAAATAAACACCTATTCAATGGAAAAACGCTACTA  >  W3110S.gb/2146484‑2146544

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