Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4581244	4581977	734	4 [3]	[3] 4	[hsdS]–[hsdM]	[hsdS],[hsdM]

CGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAAC  >  NC_000913/4581752‑4582182                                                                                                                                                                                                                                   |                                                                                                                                                                                                             cGCTGTAGACAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCAAACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCAGAAACTCGGCATATTGGTACGCAGGTCATACACCAACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGGCACCGTCCCTTTGGTAAAGAACAGAAAGTTGGTCTTcgcg                                                                                                                                                                                      >  2:162217/1‑248 (MQ=255)                                                                                                                    ggTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCg                                                                   >  1:151937/1‑251 (MQ=255)                                                                                                                                                                                     cGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAACGCACAAc  >  1:224545/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaaGAACAGCACGTTGGTCTTCACGCCCTGAGCGTCATAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGCTAATAtg                  >  1:110569/1‑189 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                             CGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAAC  >  NC_000913/4581752‑4582182

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