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

Filter HOM_REF records without GQ information #1805

Open
biozzq opened this issue May 14, 2022 · 0 comments
Open

Filter HOM_REF records without GQ information #1805

biozzq opened this issue May 14, 2022 · 0 comments

Comments

@biozzq
Copy link

biozzq commented May 14, 2022

Instructions

Feature request

When filtering the variants, I found that genotypes without GQ information have been filtered. I wonder that if can add any parameter to change this default behavior. For example, if genotypes do not have GQ or any other tags, we can assume them to be PASS.

Tool(s) involved

Tool name(s), special parameters?
FilterVcf
java -Xmx20g -Djava.io.tmpdir=${tmp} -jar picard.jar FilterVcf --MIN_DP 3 --MIN_GQ 20 --INPUT ${chr}.raw.vcf.gz --OUTPUT ${chr}.flag.vcf.gz --CREATE_INDEX true

Description

Specify whether you want a modification of an existing behavior or addition of a new capability.
Provide examples, screenshots, where appropriate.
1 23454 rs1109686627 C . 132.87 PASS AN=32;DP=568 GT:AD:DP:FT 0/0:33:33:LowGQ 0/0:38:38:LowGQ 0/0:33:33:LowGQ 0/0:35:35:LowGQ 0/0:36:36:LowGQ 0/0:36:36:LowGQ 0/0:35:35:LowGQ 0/0:45:45:LowGQ 0/0:36:36:LowGQ 0/0:35:35:LowGQ 0/0:33:33:LowGQ 0/0:34:34:LowGQ 0/0:34:34:LowGQ 0/0:35:35:LowGQ 0/0:36:36:LowGQ 0/0:34:34:LowGQ

I want to keep these non-variant sites, thus I ask for your help here. Thank you in advance.

Best wishes,
Zheng zhuqing

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