Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	minE	501672	501704	33	3 [0]	[0] 35	modA	molybdate transporter subunit

GGGAAAATCACGGTGTTCGCCGCCGCATCACTGACTAACGCAATGCAGGA  >  minE/501622‑501671                                                  | gggAAAATCACGGTGTTCGCCGCCGCATCACTGACTAACGCAATGCAGGa  >  1:2003040/1‑50 (MQ=255) gggAAAATCACGGTGTTCGCCGCCGCATCACTGACTAACGCAATGCAGGa  >  1:2432169/1‑50 (MQ=255) gggAAAATCACGGTGTTAGCCGCCGCATCACTGACTAACGCAATGCAGGa  >  1:817495/1‑50 (MQ=255)                                                  | GGGAAAATCACGGTGTTCGCCGCCGCATCACTGACTAACGCAATGCAGGA  >  minE/501622‑501671

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