Question: problem to filter freq from vcf produced with varscan
0
gravatar for fabbri.marco
2.9 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 • 889 views
ADD COMMENTlink written 2.9 years ago by fabbri.marco0
1
gravatar for Pierre Lindenbaum
2.9 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum131k wrote:

from this old file https://github.com/blachlylab/mucor/blob/master/examples_package/data/NA12878.VarScan.indels.vcf (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 2.9 years ago by Pierre Lindenbaum131k
Please log in to add an answer.

Help
Access

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