Question: SNPsift filtering problem
0
gravatar for sara.mehrabi
13 months ago by
sara.mehrabi40
Sweden
sara.mehrabi40 wrote:

Dear All.

I would like to filter my VCF file with SNPsift for DP, Quality, SRP, SAP and Effect (remove all the intergenic_region(MODIFIER||||||||||1).

The first problem is that I can not filter for all of them together, I followed the instruction, but its not working. I can just filter for one of them at the time.

the second problem is I can not filter for Effect, I tried several things, but its not working.

I appreciate your thought or experience.

Best, Sara

snpsift filtering • 441 views
ADD COMMENTlink written 13 months ago by sara.mehrabi40
0
gravatar for Guy Reeves
13 months ago by
Guy Reeves1.0k
Germany
Guy Reeves1.0k wrote:

First step, are all the variables yo are filtering with present in your vcf. Have you annotated your .vcf e.g.. with SNPeff? When working with SNPsift it is best to put in simple expressions first, check the output to make sure it makes sense and only then to attempt to combine expressions. This is especially useful for less experienced users, the syntax of very complex expressions is often not trivial. Go stepwise. Cheers Guy

ADD COMMENTlink written 13 months ago by Guy Reeves1.0k
1

First step: yes they are available. yes I did annotate with SnpEff. I am able to filter for DP, Quality, SRP, SAP one by one, but not all together. exp: I filter for DP, then filter for QUAL on the DP filtered file and so on and it worked. but still I am not able to filter for Effect. Thanks.

ADD REPLYlink written 13 months ago by sara.mehrabi40

First part are you placing each filter within (). unless specified this makes each filter an AND. So for (DP>23) (QUAL>200) only sites which pass both filters will be retained. If you want to use OR then you need to look at "|". look at the the instructions on http://snpeff.sourceforge.net/SnpSift.html

they are vey good.

About filtering the 'ANN' ' field where I guess your 'intergenic_region' annotation is have you read the section in the above link?

ADD REPLYlink written 13 months ago by Guy Reeves1.0k
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 16.09
Traffic: 157 users visited in the last hour