Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
High performance computing & Tophat & Cufflinks paolo.kunder Bioinformatics 1 01-24-2012 09:02 AM
snpEff and extracting EFF field kasthuri Bioinformatics 2 12-30-2011 01:32 PM
recommendations needed for filtering the variant file kjaja Bioinformatics 0 12-09-2011 01:47 PM variant filtering vidhya General 0 06-07-2011 01:54 AM
MAQ: Obtain quality & read depth info non-variant (referent) sites erichpowell Bioinformatics 1 12-18-2009 01:35 PM

Thread Tools
Old 11-15-2011, 01:00 AM   #1
Location: Planet Earth

Join Date: Jan 2011
Posts: 26
Default snpEff & snpSIFT variant filtering


I want to filter a .vcf file with snpSIFT and have the problem that I can't use the filter expressions described here:

I figured out that instead of:
cat variants.vcf | java -jar SnpSift.jar " ( QUAL >= 30 )" > filtered.vcf

I have to use:
java -jar SnpSift.jar filter -f variants.vcf " ( QUAL >= 30 )" > filtered.vcf

to get this done. But I dont' have a clue how to target the INFO field, where all the important information like DP, AF or in case of an snpEFF annotated file the EFF is located. When I use a expression like ( DP > 10 ) I get an error No such field 'DP'...

Anyone ?
doc.ramses is offline   Reply With Quote
Old 03-05-2014, 07:55 AM   #2
Junior Member
Location: Concepcion, Chile

Join Date: Dec 2012
Posts: 9

I have the same question:
I need to filter my vcf file by allele frequencies in order to keep only those samples that have differences between each other (polymorphic SNVs), not between them and the reference.
If you got the answer, please share it =)

Danielbenitezr is offline   Reply With Quote
Old 10-21-2014, 01:29 PM   #3
Junior Member
Location: houston

Join Date: Oct 2014
Posts: 2

Here you go!

cat VCF_file.vcf | java -jar SnpSift.jar filter "(( QUAL >= 30) & (na ID ) & (DP >= 200) & (AO >= 20) & (SAF >= 5) & (SAR >= 5))" | > Output.vcf

you can do similar commands with different portions of the INFO field. The "na ID" field is for me to output the non-SNPs.

If anyone knows how to combine the ("na ID" or "(exists ID & Eur_AF < .02)) fields, I'd appreciate the knowledge. Right now i have to do two different vcf files and cat them.
jreuther is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:42 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO