Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4578004	4580102	2099	4 [3]	[2] 4	[mcrB]–[hsdS]	[mcrB],[mcrB],symE,symR,[hsdS]

GATTAGCGATGCGGGTGTGTTGGCGCACATCCGCACCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGATTGCGTCAGGTTGTTGACGCGGGCTAAGGCGTTGTTGACCTGTTTTTCTATGGTGTCGGCGTAGGCGAAGAGTTGCTCGACGCGGCGAACGATTTCGGCTTGTTCTTTTACTGGAGGTAATAA  >  NC_000913/4579870‑4580352                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          gagtagtgATGCGGGTGGGTTGGCGAACATCCGCACCGCGCTAAATTCCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGGTAAAAATAGAATAAAATTGGTAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAATAGCGCACCGGAAAGGTGCGCCTGAAAATAATGTTCAGGACTTTTTACGTGAGGCTTTTTTACCCCCGCTATCTGTGAGTTCAGCt                                                                                                                                                                                                                                           <  2:156092/243‑1 (MQ=255)                                                                                 ccAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAAGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCg                                                                                                                                                                >  1:123451/1‑245 (MQ=255)                                                                                                                                                                                                                                          tAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGGTCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGATTGCGTCAGGTTGTTGACGCGGGCTAAGGCGTTGTTGACCTGTTTTTCTATGGTGTCGGCGTAGGCGAAGAGTTGCTCGACGCGGCGAACGATTTCGGCTTGTTCTTTTACTGGAGGtaataa  >  1:314967/1‑250 (MQ=255)                                                                                                                                                                                                                                          tAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGATTGCGTCAGGTTGTTGACGCGGGCTAAGGCGTTGTTGACCTGTTTTTCTATGGTGTCGGCGTAGGCGAAGAGTTGCTCGACGCGGCGAACGATTTCGGCTTGTTCTTTTACTGGAGGtaataa  >  1:180437/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          GATTAGCGATGCGGGTGTGTTGGCGCACATCCGCACCGCGCTAAATACCTGTATATATCATCAGTAAATATGGGGAAAGTCCAGCTAAAAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGATTGCGTCAGGTTGTTGACGCGGGCTAAGGCGTTGTTGACCTGTTTTTCTATGGTGTCGGCGTAGGCGAAGAGTTGCTCGACGCGGCGAACGATTTCGGCTTGTTCTTTTACTGGAGGTAATAA  >  NC_000913/4579870‑4580352

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