Question: Add INFO tag to one VCF if record present in a second VCF
gravatar for Sean Davis
5.7 years ago by
Sean Davis26k
National Institutes of Health, Bethesda, MD
Sean Davis26k wrote:

I have a VCF file to which I would like to add an INFO flag for every record that is present in another (LARGE) VCF file.  There are a number of tools to do filtering and another set of tools that will merge INFO fields from one VCF to another, but this isn't exactly what I need (though the latter could be adapted, I suppose).  

The use case is that I am trying to do a simple panel-of-normals filter to a set of somatic variants.  Should be simple, I know!

vcf • 2.1k views
ADD COMMENTlink modified 5.6 years ago by Pierre Lindenbaum130k • written 5.7 years ago by Sean Davis26k
gravatar for Pierre Lindenbaum
5.6 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum130k wrote:

I wrote a tool like this but it's undocumented and I should update this old pre-htsjdk code.

The compilation should be:

ant vcfvcf

It works with a tabix indexed file. A usage would be:

java  -jar dist/vcfvcf.jar -H

java -jar  dist/vcfvcf.jar ACF=CONFLICTALT  INFO=TAG1 INFO=TAG2 \
  TABIX=tabix.vcf.gz \
  IN=in.vcf > out.vcf


ADD COMMENTlink written 5.6 years ago by Pierre Lindenbaum130k
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: 1381 users visited in the last hour