Using SnpEff to extract relevant fields
2
0
Entering edit mode
7.8 years ago

Hi all , I am new to SnpEff tool, i want to produce an annotated vcf file that contains relevant columns such as Gene, Transcript and Nucleotide Change and PopFreqMax. I think it is possible using SnpSift filters but don't know how to approach it, Is there any command or script... if yes plz comment.

SNP genome gene • 4.1k views
ADD COMMENT
1
Entering edit mode
7.8 years ago
morovatunc ▴ 550

If that measure is on the VCF header. We can do that. In order to filter some of the parameters you should first determine their measure in the vcf header. For example, in the snp+indel vcf, you may filter snp by basically saying cat vcf | java -jar $snpSiftjar filter "(exists SNP)"

ADD COMMENT
0
Entering edit mode

Thank you for the reply.. Here are steps what all i have done so far, please correct me where i am going wrong.

Step 1. Downloaded SnpEff zar file, Extracted into my home directory. Step 2. Downloaded GRCh37.75 from SnpEff official website and extract it into same location Step 3. Ran command java -Xmx4g -jar snpEff.jar -v -stats ex1.html GRCh37.75 example/file.vcf > eaxmples/file.ann.vcf Step 4. Got three files as output i.e. html,annotated and txt file.

Now i am trying to run the command you shared, but getting below error "Exception in thread "main" java.lang.RuntimeException: INFO field 'SNP' not found in VCF header"

Please let me know where I am going wrong..

ADD REPLY
0
Entering edit mode
7.7 years ago
sbombin ▴ 10

Hi, I had the same problem. Actually to filter your annotated VCF file based on impact like "High" impact for example, you need to run:

 java -jar SnpSift.jar filter "ANN[0].IMPACT has 'HIGH'" input_annot.vcf > output_annot_filtered.vcf
ADD COMMENT

Login before adding your answer.

Traffic: 2536 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6