This codebase is branched off of official SnpEff version 3.6
The snpEff make script assumes a specific Eclipsey directory structure. Let's not muck with that.
mkdir ~/workspace ~/snpEff
cd ~/workspace
git clone git@github.com:CBMi-BiG/snpEff.git SnpEff
cd SnpEff
Maven/Ant/Ivy have a hard time getting Picard and Samtools jars, but they are in the Picard distribution. Alter the pom.xml versions as necessary.
wget http://downloads.sourceforge.net/project/picard/picard-tools/1.84/picard-tools-1.84.zip
unzip picard-tools-1.84.zip
mvn install:install-file -DgroupId=net.sf.samtools -DartifactId=Sam -Dversion=1.84 \
-Dpackaging=jar -Dfile=$PWD/picard-tools-1.84/sam-1.84.jar
mvn install:install-file -DgroupId=net.sf.picard -DartifactId=Picard -Dversion=1.84 \
-Dpackaging=jar -Dfile=$PWD/picard-tools-1.84/picard-1.84.jar
- 3.1h2 Introduces proper HGVS for insertions and deletions, including those which require "walking and rolling" to identify the correct indel frame.
- 3.6 Merges official SnpEff 3.6
To build snpEff-3.6-jar-with-dependencies.jar in the ~/workspace/SnpEff/target directory:
./scripts/make.sh
This is copied to ~/workplace/SnpEff/snpEff.jar per the config file
To get GRCh37 annotations:
mkdir -p data/GRCh37.64
java -jar snpEff.jar download GRCh37.64
To test (assuming you have downloaded all the test cases):
java -cp snpEff.jar \
ca.mcgill.mcb.pcingola.snpEffect.testCases.TestSuiteAll
java -Xmx4g -jar snpEff.jar hg19 -i vcf -o vcf tests/hgvs_test_in.vcf
This fork of SnpEff introduces the following annotations to VCF files
Segment(previously Exon) - a verbose description of exons (e.g.NM_152486.2.ex.3) or introns (NR_024540.1_intron_7) for all applicable transcript hitsHGVS_DNA_nomenclatureHGVS coding DNA nomenclature for exonic SNPs e.g.NM_001005484.1:c.655C>T. See http://www.hgvs.org/mutnomen/recs-DNA.html for more info.
##SnpEffVersion="SnpEff_cbmi 3.6h"
##SnpEffCmd="SnpEff hg19 -i vcf -o vcf tests/hgvs_test_in.vcf "
##INFO=<ID=EFF,Number=.,Type=String,Description="Predicted effects for this variant.Format: 'Effect ( Effect_Impact | Functional_Class | Codon_Change | Amino_Acid_change| Amino_Acid_length | Gene_Name | Gene_BioType | Coding | Transcript | Rank | Segment | HGVS_DNA_nomenclature | HGVS_protein_nomenclature [ | ERRORS | WARNINGS ])' ">
#CHROME POS ID REF ALT QUAL FILTER INFO
1 16856 . A G 100.0 PASS DP=100;EFF=DOWNSTREAM(MODIFIER|||||DDX11L1||NON_CODING|NR_046018.2||||),INTRON(MODIFIER|||||WASH7P||NON_CODING|NR_024540.1|7|NR_024540.1_intron_7||),SPLICE_SITE_DONOR(HIGH|||||WASH7P||NON_CODING|NR_024540.1|7|NR_024540.1.ex.5||)
1 69745 . C T 100.0 PASS DP=100;EFF=STOP_GAINED(HIGH|NONSENSE|Cag/Tag|Q219*|305|OR4F5||CODING|NM_001005484.1|1|NM_001005484.1.ex.1|NM_001005484.1:c.655C>T|p.Q219*)
1 865697 . A G 100.0 PASS DP=100;EFF=NON_SYNONYMOUS_CODING(MODERATE|MISSENSE|Atc/Gtc|I79V|681|SAMD11||CODING|NM_152486.2|3|NM_152486.2.ex.3|NM_152486.2:c.235A>G|p.I79V)
1 1114650 . C T 100.0 PASS DP=100;EFF=STOP_GAINED(HIGH|NONSENSE|Cga/Tga|R19*|673|TTLL10||CODING|NM_001130045.1|4|NM_001130045.1.ex.4|NM_001130045.1:c.55C>T|p.R19*),UPSTREAM(MODIFIER||||404|TTLL10||CODING|NM_153254.2||||)