SEQanswers

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 08:02 AM
snpEff and extracting EFF field kasthuri Bioinformatics 2 12-30-2011 12:32 PM
recommendations needed for filtering the variant file kjaja Bioinformatics 0 12-09-2011 12:47 PM
vcfutils.pl variant filtering vidhya General 0 06-07-2011 12:54 AM
MAQ: Obtain quality & read depth info non-variant (referent) sites erichpowell Bioinformatics 1 12-18-2009 12:35 PM

Reply
 
Thread Tools
Old 11-15-2011, 12:00 AM   #1
doc.ramses
Member
 
Location: Planet Earth

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

Hi,

I want to filter a .vcf file with snpSIFT and have the problem that I can't use the filter expressions described here: http://snpeff.sourceforge.net/SnpSift.html

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, 06:55 AM   #2
Danielbenitezr
Junior Member
 
Location: Concepcion, Chile

Join Date: Dec 2012
Posts: 9
Default

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 =)

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

Join Date: Oct 2014
Posts: 2
Default

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
Reply

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 11:46 AM.


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