Question: Udpate the FILTER field in VCF to PASS if the variants is represented as true in another set of VCF?
0
gravatar for kirannbishwa01
3.1 years ago by
kirannbishwa011.2k
United States
kirannbishwa011.2k wrote:

I have a VCF files called from several samples (Say Vcf01), and also have a VCF files with sites that are true variants (Say Vcf02). So, what I want to do is update the FILTER field for the sites in Vcf01 to PASS if the variants is represented in Vcf02.

I have checked several VCF tools and also GATK but this functionality is lacking.

Any ideas?

Thanks,

pass snp vcf • 1.2k views
ADD COMMENTlink modified 3.1 years ago by Pierre Lindenbaum133k • written 3.1 years ago by kirannbishwa011.2k
1
gravatar for Pierre Lindenbaum
3.1 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum133k wrote:

gatk variant annotator https://software.broadinstitute.org/gatk/gatkdocs/3.8-0/org_broadinstitute_gatk_tools_walkers_annotator_VariantAnnotator.php

and peek any INFO from the second vcf:

   --resource:foo resource.vcf \
   -E foo.AF \

followed by

gatk variant filtration: https://software.broadinstitute.org/gatk/documentation/tooldocs/current/org_broadinstitute_gatk_tools_walkers_filters_VariantFiltration.php

where you test for the presence of foo.AF.

On my side, I've written http://lindenb.github.io/jvarkit/VcfIn.html

with options:

-A, --allalt
  ALL user ALT must be found in VCF-database ALT
  Default: false
-fi, --filterin
  Do not discard variant but add this FILTER if the variant is found in 
  the database
  Default: <empty string>
-fo, --filterout
  Do not discard variant but add this FILTER if the variant is NOT found 
  in the database
  Default: <empty string>
ADD COMMENTlink written 3.1 years ago by Pierre Lindenbaum133k
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: 2234 users visited in the last hour
_