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]

ATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAAT  >  NC_000913/4577753‑4578241                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               aTGATCTCGGTGTTAGGATTGTAATCAAGCACAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGTCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTTCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAg                                                                                                                                                                                                                                                <  2:67834/251‑1 (MQ=255)                                                                                                                                                                             aGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCTCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTc                                                                     >  1:13903/1‑250 (MQ=255)                                                                                                                                                                              gATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCa                                                                    >  1:66446/1‑250 (MQ=255)                                                                                                                                                                                                                                               cggttgttGATAGGGGTTAGAAAAAAAATATTCGTCGAGTAAAGGGGTGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGGTCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGATTCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCTTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTGTTATTCAGTAAAAAAt  <  2:13903/249‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               ATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGGGGTCATCAAAGAAATATTCTTCGAGTAAAGGGGCGATATCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAAT  >  NC_000913/4577753‑4578241

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