Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ minE 1039568 1039606 39 25 [1] [0] 32 nth DNA glycosylase and apyrimidinic (AP) lyase

AGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTC  >  minE/1039506‑1039568
                                                             | 
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:374233/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:998405/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:808989/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:772152/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:736348/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:73534/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:71842/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:694895/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:683829/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:678037/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:621141/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:612919/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:1007170/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:345045/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:303438/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:287911/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:184175/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:16412/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:116769/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:115579/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:114330/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:113101/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:108715/1‑62 (MQ=255)
agcagcATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt   >  1:1009200/1‑62 (MQ=255)
agcagcATAATG‑‑GAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc  >  1:989720/1‑61 (MQ=255)
                                                             | 
AGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTC  >  minE/1039506‑1039568

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