Question: Snpsift to filter a number of genes
gravatar for
4.2 years ago by
United States
krithika.gu260 wrote:


I have a VCF file, and would like to pull out all variants in various genes of interest (I have about 200 genes)

I created a txt file with all my genes of interest.

cat variants.vcf | java -jar SnpSift.jar filter --set input.txt "EFF[*].GENE in SET[0]" > filtered.vcf

This seemed to filter out all the variants, and I got an empty VCF File without any variants.

  • Is the SET function not meant to be used like this ?
  • Do I have to use a loop instead ?

Could someone help with this ?

Thanks, K

ADD COMMENTlink written 4.2 years ago by krithika.gu260

old questioun, but I need exactly the same and get it to work using ANN[*].GENE

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by Jimbou730
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: 1833 users visited in the last hour