Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ minE 153475 153541 67 3 [0] [0] 25 hemL glutamate‑1‑semialdehyde aminotransferase

GGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCG  >  minE/153408‑153474
                                                                  |
ggTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTcggcg  >  1:1569469/1‑67 (MQ=255)
ggTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTcggcg  >  1:444568/1‑67 (MQ=255)
ggTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTcggcg  >  1:782664/1‑67 (MQ=255)
                                                                  |
GGTCACCAGTTGCGCCATTTTCACTTCCATTTCGGTTGGTGCACCAAAGCTTAAACCACGCTCGGCG  >  minE/153408‑153474

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