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

For joint variant calling, produce a final VCF file with filtered variants (those that PASS filters) #1102

Closed
amizeranschi opened this issue Jun 13, 2023 · 0 comments
Assignees
Labels
bug Something isn't working

Comments

@amizeranschi
Copy link
Contributor

Description of feature

Following this discussion on Slack, after running joint calling and VQSR, it would be helpful to produce a final VCF file with filtered variants, i.e. those that have PASS on their filter column, and not . or anything else. Something like the following bcftools command:

bcftools view -f PASS joint_germline_recalibrated.vcf.gz > joint_germline_recalibrated_final.vcf.gz

The current behavior is to have a mix of both calibrated and uncalibrated variants in joint_germline_recalibrated.vcf.gz, which means that there are two records for each variant.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
Archived in project
Development

No branches or pull requests

3 participants