added ability to set max x & y-axis values in chartdirector output; 02-vcf-analysis script loads ngs.locations db table for use in start position exclusions
extended download function to sample data, split existing vep & non-cds downloads into separate files; moved update_vep_table() function into new save_data() method
simplified args passed in NGS.pm method calls - src_data & data_file are equivalent so replaced all refs to former with latter; simplified config block in _parse()
wrapped call to NGS::VEP in try/catch to handle die via 500.tt - works for tests as well; test for invalid filename; display list of valid filenames in upload page
11 years ago
README.txt
http://www.ensembl.org/info/docs/variation/vep/vep_script.html
Setup/update VEP - delete scripts/variant_effect_predictor dir
extract variant_effect_predictor.tar.gz into scripts dir
run scripts/variant_effect_predictor/INSTALL.pl - skip cache.
Cache manually installed/updated (needs to match VEP version) from:
ftp://ftp.ensembl.org/pub/release-<xx>/variation/VEP/
see README.txt in NGS/script dir for more details
Cache directory structure:
vep dir containing core & refseq dirs, each with homo_sapiens/xx cache
* on dev box - symlink /home/raj/.vep -> /media/sf_WIN_DRIVE/vep
* on deployment - /home/raj/.vep
to run from command-line:
perl script/variant_effect_predictor.pl -config=script/vep.ini \
-i=t/data/myeloid_variants.vep -o=output.txt --polyphen=b --sift=b \
--check_existing --coding_only --regulatory --dir=refseq/core
or:
perl script/variant_effect_predictor.pl -config=/home/raj/apps/NGS/script/vep.ini \
-o=output.txt --dir=/home/raj/.vep/refseq --polyphen=b --sift=b --check_existing \
--coding_only --regulatory -i=/tmp/20_06_14.vep
PERFORMANCE:
a) 4055-row vep input with unsorted chromosome order on 163.160.171.48:
fork=0 255 sec
fork=2 580 sec
fork=4 620 sec
b) 1028-row vep input with alpha-numerically sorted chromosome order on 163.160.171.48:
fork=0 95 sec
fork=2 65 sec (optimal)
fork=3 68 sec
fork=4 70 sec
fork=5 75 sec
fork=6 75 sec
performance much WORSE on dev server with fork enabled
splitting a 1110-row vep input into separate chromosome input files slightly
REDUCED performance (77sec vs 66sec) on dev
