Question: problem to filter freq from vcf produced with varscan
gravatar for fabbri.marco
3.2 years ago by
fabbri.marco0 wrote:

I have a VCF file obtained from the Varscan2 tool. I want to filter the variants based on "FREQ". I used SnpSift (SnpSiftVersion="SnpSift 4.2 (build 2015-12-05), by Pablo Cingolani"). When I used the command: java -Xmx64g -jar ${SNPSIFT_JAR} extractFields $svcf "GEN[1].FREQ" I got the right field of my vcf file: 0.52% 0% 0.72% 0% 0.31% 0.21% 0.31% 0.31%

But when I tried to filter variant with a freq greter than 0.1, no variant pass the filter. I ran this command: java -Xmx64g -jar ${SNPSIFT_JAR} filter " ( GEN[0].FREQ > 0.1 ) " $svcf

I got the feeling that freq is not an integer.

Any suggstion?

thank you Marco

snpsift varscan2 filter vcf • 935 views
ADD COMMENTlink written 3.2 years ago by fabbri.marco0
gravatar for Pierre Lindenbaum
3.2 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum133k wrote:

from this old file (2015)

FREQ is a string

##FORMAT=<ID=FREQ,Number=1,Type=String,Description="Variant allele frequency">

check this in your VCF header.

ADD COMMENTlink written 3.2 years ago by Pierre Lindenbaum133k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 992 users visited in the last hour