接昨天pvactools 安裝蜒秤，今天測(cè)試的目的是軟件跑通捍岳。
VEP install 可以直接下載cache files , FASTA files and plugins . 昨天由于下載過(guò)慢，轉(zhuǎn)為手動(dòng)下載。
注意：cache 文件 分為VEP 和 VEP_indexed ,官方解釋，這兩個(gè)文件都可以使用悴侵，indexed cache 包注釋速度非常快拭嫁。 #需要重新下載indexed版本
#indexed_Vep
curl -O https://ftp.ensembl.org/pub/release-110/variation/indexed_vep_cache/homo_sapiens_vep_110_GRCh38.tar.gz
#vep
wget -c https://ftp.ensembl.org/pub/release-110/variation/vep/homo_sapiens_merged_vep_110_GRCh38.tar.gz
wget -c https://ftp.ensembl.org/pub/release-110/variation/vep/homo_sapiens_refseq_vep_110_GRCh38.tar.gz
wget -c https://ftp.ensembl.org/pub/release-110/variation/vep/homo_sapiens_vep_110_GRCh38.tar.gz
image.png
今天分析的時(shí)候遇到了新問(wèn)題可免。
#啟動(dòng)griffithlab/pvactools 容器
docker run -it -v /share/home/test/00test/neoantigens_test/pvacseq_test:/test_data griffithlab/pvactools
#下載vep plugin
#pvactools 鏡像中已經(jīng)包含VEP plugins
pvacseq install_vep_plugin /test_data/vep_plugin
###vep 測(cè)試 ok
docker run -it \
-v /share/home/test/00test/neoantigens_test/pvacseq_test:/test_data \
-v /share/home//scripts/datas/ensembl-vep:/ensembl-vep \
ensemblorg/ensembl-vep vep \
--input_file /test_data/sample1_gDNA.vardict.vcf \
--output_file /test_data/sample1_gDNA.vep3.vcf \
--format vcf --vcf --symbol --terms SO --tsl --biotype --hgvs \
--fasta /ensembl-vep/hg38/hg38.fa \
--offline --cache /ensembl-vep/ --plugin Frameshift --plugin Wildtype \
--dir_plugins /test_data/vep_plugin \
--pick --transcript_version
##
### pvacseq 測(cè)試可以跑通
docker run -it \
-v /share/home/dujl/test/00test/neoantigens_test/pvacseq_test:/test_data \
griffithlab/pvactools pvacseq run \
/test_data/sample1_gDNA.vep3.vcf \
sample1_gDNA \
HLA-A*29:02,HLA-B*45:01,DRB1*04:05 \
all \
/test_data/outer2 \
-e1 8,9,10 \
-e2 15 \
--iedb-install-directory /opt/iedb
位點(diǎn)多，分析速度較慢做粤；
多線程后浇借，速度超級(jí)加倍~
pvacseq 使用過(guò)濾后的vcf ，HLA分型結(jié)果對(duì)計(jì)算結(jié)果影響怕品？
