brentp / duphold

don't get DUP'ed or DEL'ed by your putative SVs.
MIT License
100 stars 9 forks source link

missing insertion #46

Open Danisov opened 2 years ago

Danisov commented 2 years ago

Hi @brentp I have a question about missing insertion variants. I used sample HG002 for variant calling with smoove using the command below. But it seems strange to me that I can't find the variants of Insertions? Do you have any idea why? does it have a notation instead of INS? Thank you in advance for your answer. Note: I also used a WGS of my samples and I still haven't found the Insertions.

$ singularity run smoove.sif smoove call \
-x \
--name my_project \
--exclude ${BED} \
--fasta ${hg38} \
-p 1 \
--genotype \
 ${input} 
brentp commented 2 years ago

what does this mean:

I can't find the variants of Insertions

? Do you mean duphold isn't annotating your insertions? Can you show a variant VCF line output from duphold that doesn't show what you expect?

Danisov commented 2 years ago

My apologies I wanted to write this issue in smoove not on duphold. Do you want me to delete it and re-write it in the smoove issues?

brentp commented 2 years ago

Here is fine, I just don't have enough information to help you.

Danisov commented 2 years ago

I just noticed that the insertions generated with smoove or lumpy are in BND format. Thanks @brentp