Question: extract polyformisms non synonymous using SnpSift
1
gravatar for zion22
19 months ago by
zion2250
merilla
zion2250 wrote:

Hi, sorry to bother you. I ask in this forum why maybe someone can help me, some of you know how I could extract polyformisms non synonymous using SnpSift? and if you have an example of the command I appreciated it very much. Thank you

snp sequence next-gen genome • 342 views
ADD COMMENTlink modified 19 months ago by Arup Ghosh2.7k • written 19 months ago by zion2250
4
gravatar for Arup Ghosh
19 months ago by
Arup Ghosh2.7k
India
Arup Ghosh2.7k wrote:

The SnpSift filter option can be used to extract variants based on annotation and other parameters.

java -jar SnpSift.jar filter "( EFF[0].EFFECT = 'NON_SYNONYMOUS_CODING' )"  annotated.vcf > non_syn.vcf

Ref: http://snpeff.sourceforge.net/SnpSift.html#filter

ADD COMMENTlink modified 19 months ago • written 19 months ago by Arup Ghosh2.7k

Hi. Thank you so much for helping me. but I have a problem, when I run the next command:

java -jar SnpSift.jar filter "( EFF[0].EFFECT = 'NON_SYNONYMOUS_CODING' )"  my_original_file.vcf > NON_SYNONYMOUS_file.vcf

when I compare the two files, original and filtered, are exactly the same, even when I change the filter parameter NON_SYNYMOUS_CODING to SYNONYMOUS_CODING the output file is the same. if you know what I could do it would be great since I have to do it with 10 other vcf files and I have the same problem

ADD REPLYlink modified 18 months ago • written 18 months ago by zion2250

Can you share a snippet of the vcf file including the header?

ADD REPLYlink written 18 months ago by Arup Ghosh2.7k
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 2.3.0
Traffic: 1075 users visited in the last hour