search

HKU-BAL / Clair3

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

Correct variant found by pileup is changed after full alignment #316

Closed f-ferraro closed 1 month ago

f-ferraro commented 2 months ago

Hi, Thanks for your work with clair3!

We are running clair3 v1.0.9 via docker and we encountered an expected result that we would like your opinion about.

We have a previously known and validated variant in EHMT1.

We recently sequenced the same sample with ONT (SUP basecalling, dna_r10.4.1_e8.2_400bps_sup@v5.0.0) and we didn't find this variant in the clair3 merge_output.vcf. We initially suspected something went wrong with the sequencing itself but the variant is clearly visible in the bam. We however realized that the variant is found in the pileup vcf but it's changed into another variant after the full alignment. Below you can see the screenshot of the variant and zooming out a bit more of the region. Do you perhaps have any suggestions? Thanks in advance! Federico

image

image

aquaskyline commented 2 months ago

Weird, would you be able to send us a minibam (only the reads spanning the variant of interest) so we could look into the details.