ensembl-vep (docker:下載使用plugins、fasta宣谈、caches)

教程http://asia.ensembl.org/info/docs/tools/vep/script/vep_download.html

docker pull ensemblorg/ensembl-vep?

docker run -ti ensemblorg/ensembl-vep ./vep (試驗是否下載可運行)

mkdir /home/vep_data

chmod a+rwx /home/vep_data? ??

docker run -tiv /home/vep_data:/opt/vep/ .vep ensemblorg/ensembl-vep perl INSTALL.pl (按提示下載cache等,啰里啰嗦键科,放棄闻丑!)

借此機會,開個新終端勋颖,查看下image的結(jié)構(gòu)

docker ps (查看下container_id)·? ??

docker exec -it container_id bash (進(jìn)入image內(nèi)部嗦嗡,看下examples下內(nèi)容吧)

docker stop container_id(關(guān)閉)

docker run -t -i -v /home/vep_data:/opt/vep/ .vep ensemblorg/ensembl-vep perl INSTALL.pl -a cfp -s homo_sapiens -y GRCh37 -g all (自動下載所有東西,但是實在太慢了胺沽帷=募馈!茄厘!放棄)

(a?(or?--AUTO) with the parameters?c?(cache),?f?(fasta) and?p?(plugins) . If you wish to include the?VEP plugins, add the 'p' value to the?-a?flag and the?--PLUGINS?(or?-g) flag as well:)

本地文件夾結(jié)構(gòu):與容器共享數(shù)據(jù)路徑下矮冬,創(chuàng)建input,output次哈,Plugins文件夾(名稱胎署、路徑只是建議,隨便啦)


1)下載caches,?ftp://ftp.ensembl.org/pub/release-96/variation/VEP/homo_sapiens_vep_96_GRCh38.fa.gz

在/home/vep_data下解壓窑滞,獲得homo_sapiens文件夾/home/vep_data/homo_sapiens/96_GRCh38

2)下載fasta數(shù)據(jù)

ftp://ftp.ensembl.org/pub/release-96/fasta/homo_sapiens/dna/Homo_sapiens.GRCh38.dna.primary_assembly.fa.gz

在/home/vep_data/fasta下解壓

3)去github下載plugin插件琼牧,https://github.com/Ensembl/VEP_plugins,例如G2P.pm,

https://asia.ensembl.org/info/docs/tools/vep/script/vep_plugins.html#g2p參看下載G2P.csv?https://www.ebi.ac.uk/gene2phenotype/downloads?

將G2P.pm和2個G2P.csv放入/home/vep_data/Plugins中哀卫,

試一下

docker run -u root -t -i -v /home/vep_data:/opt/vep/.vep ensemblorg/ensembl-vep \

./vep --cache --offline --format vcf --force_overwrite \--dir_cache /opt/vep/.vep/ \

--hgvs --fasta /opt/vep/.vep/fasta/Homo_sapiens.GRCh38.dna.primary_assembly.fa \

--input_file /opt/vep/.vep/input/input.vcf \

--output_file /opt/vep/.vep/output/output.vcf \

--dir_plugins /opt/vep/.vep/Plugins \

--plugin G2P,file=/opt/vep/.vep/Plugins/CancerG2P.csv,html_report=/opt/vep/.vep/output/report.html,txt_report=/opt/vep/.vep/output/report.txt

--assembly GRCh38

(--vcf, Writes output in?VCF format.?)

report.txt長這樣:

report.html樣子長這樣:


最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末巨坊,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子此改,更是在濱河造成了極大的恐慌抱究,老刑警劉巖,帶你破解...
    沈念sama閱讀 219,427評論 6 508
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件带斑,死亡現(xiàn)場離奇詭異鼓寺,居然都是意外死亡勋拟,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 93,551評論 3 395
  • 文/潘曉璐 我一進(jìn)店門妈候,熙熙樓的掌柜王于貴愁眉苦臉地迎上來敢靡,“玉大人,你說我怎么就攤上這事苦银⌒ル剩” “怎么了?”我有些...
    開封第一講書人閱讀 165,747評論 0 356
  • 文/不壞的土叔 我叫張陵幔虏,是天一觀的道長纺念。 經(jīng)常有香客問我,道長想括,這世上最難降的妖魔是什么陷谱? 我笑而不...
    開封第一講書人閱讀 58,939評論 1 295
  • 正文 為了忘掉前任,我火速辦了婚禮瑟蜈,結(jié)果婚禮上烟逊,老公的妹妹穿的比我還像新娘。我一直安慰自己铺根,他們只是感情好宪躯,可當(dāng)我...
    茶點故事閱讀 67,955評論 6 392
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著位迂,像睡著了一般访雪。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上掂林,一...
    開封第一講書人閱讀 51,737評論 1 305
  • 那天冬阳,我揣著相機與錄音,去河邊找鬼党饮。 笑死肝陪,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的刑顺。 我是一名探鬼主播氯窍,決...
    沈念sama閱讀 40,448評論 3 420
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼蹲堂!你這毒婦竟也來了狼讨?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,352評論 0 276
  • 序言:老撾萬榮一對情侶失蹤柒竞,失蹤者是張志新(化名)和其女友劉穎政供,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,834評論 1 317
  • 正文 獨居荒郊野嶺守林人離奇死亡布隔,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,992評論 3 338
  • 正文 我和宋清朗相戀三年离陶,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片衅檀。...
    茶點故事閱讀 40,133評論 1 351
  • 序言:一個原本活蹦亂跳的男人離奇死亡招刨,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出哀军,到底是詐尸還是另有隱情沉眶,我是刑警寧澤,帶...
    沈念sama閱讀 35,815評論 5 346
  • 正文 年R本政府宣布杉适,位于F島的核電站谎倔,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏猿推。R本人自食惡果不足惜片习,卻給世界環(huán)境...
    茶點故事閱讀 41,477評論 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望彤守。 院中可真熱鬧,春花似錦哭靖、人聲如沸具垫。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,022評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽筝蚕。三九已至,卻和暖如春铺坞,著一層夾襖步出監(jiān)牢的瞬間起宽,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,147評論 1 272
  • 我被黑心中介騙來泰國打工济榨, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留坯沪,地道東北人。 一個月前我還...
    沈念sama閱讀 48,398評論 3 373
  • 正文 我出身青樓擒滑,卻偏偏與公主長得像腐晾,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子丐一,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 45,077評論 2 355

推薦閱讀更多精彩內(nèi)容