Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Using PacBio/pbsv VCF for delly force genotyping (-v)? #382

Open
bbimber opened this issue May 21, 2024 · 0 comments
Open

Using PacBio/pbsv VCF for delly force genotyping (-v)? #382

bbimber opened this issue May 21, 2024 · 0 comments

Comments

@bbimber
Copy link

bbimber commented May 21, 2024

Hello,

We have a dataset of SVs created from PacBio data using pbsv. We'd like to take those SVs and use delly to genotype a dataset of Illumina WGS samples at these specific sites. We made a sites-only VCF from the PacBio data and provided this VCF as "-v" for delly call. This command dies with the error "Delly genotyping requires local SV assembly (INFO/CONSENSUS) and breakpoint (INFO/CONSBP) introduced in delly v1.1.7!", which maps to this code:

std::cerr << "Error: Delly genotyping requires local SV assembly (INFO/CONSENSUS) and breakpoint (INFO/CONSBP) introduced in delly v1.1.7!" << std::endl;

I assume the problem is that the VCF lack from INFO annotations delly expects.

My question is: is there a way to use a VCF generated by another caller (pbsv in particular) as the whitelist of delly's forced genotyping?

Thanks for any help.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant