Does The Vcffilter Tool Work On String Type Info Tags?
0
0
Entering edit mode
10.6 years ago
dwaggott ▴ 290

Hey,

vcflib looks great! To get started I need a bit of direction on filtering using vcffilter. Specifically, does it work on string type info tags? Also, can one group expressions similar to test5 below? I tried it on dbSNP138 clinical and ran into some problems.

Best, Daryl

### dbSNP138
#ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/

wget ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/clinvar_20130808.vcf.gz 
wget ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606/VCF/clinvar_20130808.vcf.gz.tbi

### test1 and test2 work but not test3-5
vcffilter -f "PM" clinvar_20130808.vcf.gz > test1.vcf  # PM is a flag
vcffilter -f "COMMON = 1" clinvar_20130808.vcf.gz > test2.vcf  # COMMON is a integer
vcffilter -f "CLNSIG = 4" clinvar_20130808.vcf.gz > test3.vcf #  CLINSIG and CLNORIGIN are strings
vcffilter -f "CLNORIGIN = 4" clinvar_20130808.vcf.gz > test4.vcf 
vcffilter -f "(COMMON = 0) & (MUT | PM | CLNSIG = 3 | CLNSIG = 4 | CLNORGIN = 4 | CLNORGIN = 32)"   00-All.vcf.gz > test5.vcf
vcf • 2.4k views
ADD COMMENT
0
Entering edit mode

So the problem is that you have an integer in your filter, but it's a String type in the file...

ADD REPLY

Login before adding your answer.

Traffic: 2917 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