search

HKU-BAL / Clair3

Clair3 - Symphonizing pileup and full-alignment for high-performance long-read variant calling
246 stars 27 forks source link

Variants below min_coverage still being called #262

Closed mproberts99 closed 8 months ago

mproberts99 commented 10 months ago

Hi,

I have set the min_coverage parameter to 10, but am still observing calls with read counts below that number in the final VCF. When I up the parameter to 20, then I see those have disappeared. Is there a range around the min_coverage that Clair3 still allows that would explain the presence of these calls? I have attached the VCF and the log. run_clair3.log merge_output.vcf.gz

zhengzhenxian commented 10 months ago

Hi, @mproberts99,

Thanks for the log. Seem all low-coverage variants were reported by the full-alignment model. We will check and get back to you soon.

aquaskyline commented 9 months ago

Hi, @mproberts99, sorry for the late response, we have checked and confirmed it is a problem in our code. The problem doesn't affect the correctness of each variant but allows some variants below the minimum coverage to be outputted. We will solve the problem in the next update, but before that, a hard filter on DP after calling could also solve the problem.

aquaskyline commented 8 months ago

Fixed in v1.0.6