Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	281537	281816	280	2 [1]	[1] 2	yagA	integrase core domain‑containing protein YagA

CGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCG  >  NC_000913/281812‑281958      |                                                                                                                                              cggcggCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGGGCGATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTGCTGCGGGGCGAACAAAACAAACGCAGTACGTAATGACATGGTAtct      >  1:172857/1‑143 (MQ=255)     agcACCTTCCTGAGCCCAGCGCTGAAGCCACTTGTAGCCGGTGGCAGGCGAAATGCCGAAACGACGGCAGAGGGAACGGATGTTCGCCCCGGCCTGCGAGGCGAACAAAACAAACTCAGTACGTCATGACATGGTATCTcgcg  <  2:172856/142‑1 (MQ=255)      |                                                                                                                                              CGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCG  >  NC_000913/281812‑281958

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