How to filter VEP results of "missing values"
1
0
Entering edit mode
8 weeks ago
Yiming ▴ 10

Hello,

I have a VEP annotated VCF files with CADD and LOFTEE annotations. I need to filter in results with either CADD score >= 20 OR LOFTEE HC without any flag

My code is as following

filter_vep -i VEP_annotated.vcf -o VEP_filtered.vcf --format vcf --force_overwrite --filter "CADD_PHRED >= 20 or (LoF = HC and LoF_filter = . and LoF_flags = .)"

I guess the problem is that "LoF_filter" and "LoF_flags" should be blank (no value) but I am not sure how I should express this in the code (I tried = "" but still not working)?

Thank you very much

filter VEP • 279 views
ADD COMMENT
0
Entering edit mode

Have you tried with a space between the " "?

ADD REPLY
0
Entering edit mode

The VEP dose not seem to recognize " ". It returned "MSG: ERROR: Error parsing filter string - incomplete parentheses sets?"

ADD REPLY
1
Entering edit mode
8 weeks ago
Yiming ▴ 10

Thank you so much!

I found an error with my code. And this works:

filter_vep -i VEP_annotated.vcf -o VEP_filtered.vcf --format vcf --force_overwrite --filter "(CADD_PHRED >= 20 and LoF = "") or (LoF = HC and LoF_filter = "" and LoF_flags = "")"
ADD COMMENT

Login before adding your answer.

Traffic: 1880 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6