Skip to content

Commit

Permalink
Merge pull request #512 from vasiliosz/patch-1
Browse files Browse the repository at this point in the history
Fix workflow example in docs
  • Loading branch information
arq5x committed Jul 24, 2015
2 parents bef96ce + 9381610 commit 25effdd
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions docs/index.rst
Original file line number Diff line number Diff line change
Expand Up @@ -61,20 +61,20 @@ A workflow for the above steps is given below.
# setup
VCF=/path/to/my.vcf
NORMVCF=/path/to/my.norm.vcf
NORMVCF=/path/to/my.norm.vcf.gz
REF=/path/to/human.b37.fasta
SNPEFFJAR=/path/to/snpEff.jar
# decompose, normalize and annotate VCF with snpEff.
# NOTE: can also swap snpEff with VEP
#NOTE: -classic and -formatEff flags needed with snpEff >= v4.1
zless $VCF \
| sed 's/ID=AD,Number=./ID=AD,Number=R/'
| sed 's/ID=AD,Number=./ID=AD,Number=R/' \
| vt decompose -s - \
| vt normalize -r $REF - \
| java -Xmx4G -jar $SNPEFFJAR -formatEff -classic GRCh37.75 \
| java -Xmx4G -jar $SNPEFFJAR GRCh37.75 -formatEff -classic \
| bgzip -c > $NORMVCF
tabix $NORMVCF
tabix -p vcf $NORMVCF
# load the pre-processed VCF into GEMINI
gemini load --cores 3 -t snpEff -v $NORMVCF $db
Expand Down

0 comments on commit 25effdd

Please sign in to comment.