Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	1248726	1248732	7	4 [0]	[1] 4	treA	periplasmic trehalase

AGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGC  >  W3110S.gb/1248655‑1248725                                                                       | aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc  <  1:1289467/71‑1 (MQ=255) aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc  <  1:1636716/71‑1 (MQ=255) aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc  <  1:843216/71‑1 (MQ=255) aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc  <  1:946808/71‑1 (MQ=255)                                                                       | AGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGC  >  W3110S.gb/1248655‑1248725

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