Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NZ_CP009273	1206069	1206091	23	3 [1]	[1] 2	mcrA	5‑methylcytosine‑specific endonuclease McrA

GTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAAT  >  NZ_CP009273/1205979‑1206074                                                                                          |       gtagtaTATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGAtt        <  1:18314/90‑1 (MQ=255) gtagtaTATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGAtt        <  1:84382/90‑1 (MQ=255)       tatCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAAt  >  1:137308/1‑90 (MQ=255)                                                                                          |       GTAGTATATCTGGCGTCATCATCAGTCAGAAAACATATGCATTCTTTGGATGAAAGAAAAATCCATCCTGGTGAATATTTTACTTTGATTGGTAAT  >  NZ_CP009273/1205979‑1206074

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