Question: How to use snpEff filter for getting particluar column in annotated file
0
gravatar for bishtamit078
2.8 years ago by
bishtamit0780 wrote:

I am trying to add column such as Transcript, Nucleotide Change and PopFreqMax in my output annotated file. Here are the Steps that i have performed.

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.

After successfully completing these steps i got an annotated file, but it does not have any of the field mention above. I also tried "cat cancer.vcf | java -jar SnpSift.jar filter "(exists SNP)"" but got error "Exception in thread "main" java.lang.RuntimeException: INFO field 'SNP' not found in VCF header".

Any command or script that i need to run ?

please help me out in solving this problem.

snp genome gene • 1.7k views
ADD COMMENTlink modified 2.8 years ago by Medhat8.2k • written 2.8 years ago by bishtamit0780
1
gravatar for Medhat
2.8 years ago by
Medhat8.2k
Texas
Medhat8.2k wrote:

This fields "Transcript, Nucleotide Change and PopFreqMax" not added by snpEff.

Only these fields under ANN:

ALLELE (alias GENOTYPE) EFFECT (alias ANNOTATION): Effect in Sequence ontology terms (e.g. 'missense_variant', 'synonymous_variant', 'stop_gained', etc.) IMPACT: { HIGH, MODERATE, LOW, MODIFIER } etc...

For more please see the documentation http://snpeff.sourceforge.net/SnpEff_manual.html#input

There is No filed called SNP in the header so you can not use

java -jar SnpSift.jar filter "(exists SNP)""

ADD COMMENTlink modified 26 days ago • written 2.8 years ago by Medhat8.2k
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: 872 users visited in the last hour