Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4586138	4586149	12	6 [3]	[0] 21	[hsdS]–[hsdM]	[hsdS],[hsdM]

GATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCA  >  W3110S.gb/4586079‑4586143                                                           |       gATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCg        >  1:1348554/1‑59 (MQ=255) gATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCg        >  1:1774172/1‑59 (MQ=255) gATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCg        >  1:932235/1‑59 (MQ=255) gATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTc   >  1:1596941/1‑64 (MQ=255) gATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTc   >  1:1949169/1‑64 (MQ=255) gATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCa  >  1:613191/1‑65 (MQ=255)                                                           |       GATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCA  >  W3110S.gb/4586079‑4586143

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