This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Does The Vcffilter Tool Work On String Type Info Tags?

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

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

0 answers

No answers yet.

Log in to answer this question.