Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ minE 1717647 1717671 25 3 [1] [0] 3 glyA serine hydroxymethyltransferase

GTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAA  >  minE/1717585‑1717663
                                                             |                 
gTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGccc                   <  1:28415/62‑1 (MQ=255)
gTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGccc                   <  1:499165/62‑1 (MQ=255)
                 aCAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGaa  >  1:41702/1‑62 (MQ=255)
                                                             |                 
GTGCGCCATATCAACGAACAGGTAAGCACCGATGCTGTCAGCGATTTCACGCATTTTCGCCCAGTCCACCACGCCGGAA  >  minE/1717585‑1717663

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