前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >1 wes相关软件下载与安装并添加环境变量

1 wes相关软件下载与安装并添加环境变量

作者头像
Y大宽
发布2019-06-03 10:52:43
1.4K0
发布2019-06-03 10:52:43
举报
文章被收录于专栏:Y大宽

1创建wes环境

先下载安装conda,因为它是.sh文件,直接bash即可,然后

代码语言:javascript
复制
conda create -n wes python=2 bwa
conda info --envs
source activate wes#激活wes环境

2 所有需要的软件安装在wes环境下

代码语言:javascript
复制
conda install sra-tools
conda install samtools
conda install bcftools
conda install vep
conda install snpEFF
conda install multiqc
conda install qualimap

可以一次性安装

代码语言:javascript
复制
conda install -y vep bcftools multiqc

下载gatk4,用迅雷下载比较快。GATK4下载地址 或者直接wget下载(我的速度慢)

代码语言:javascript
复制
wget https://github.com/broadinstitute/gatk/releases/download/4.1.2.0/gatk-4.1.2.0.zip

下载完成或迅雷下载复制到gatk4文件夹,解压

代码语言:javascript
复制
$ unzip gatk-4.1.2.0.zip
$ cd gatk-4.1.2.0
$ ./gatk

注意,需要下载100G的必备数据才能使用GATK。

3 gatk加入环境变量

代码语言:javascript
复制
$ echo 'export PATH=/home/kelly/biosoft/gatk/gatk-4.1.2.0/:$PATH' >>~/.bashrc
$ source ~/.bashrc
$ gatk --help

4 熟悉参考基因组及必备的数据库

4.1 参考基因组下载

我是从服务器上下载下来放本地电脑了 下载方式1:直接去gatk官网下载,下载链接为ftp://ftp.broadinstitute.org/bundle/

image.png

下载方式2:也是官网,但通过ftp匿名登录下载

代码语言:javascript
复制
location: ftp.broadinstitute.org/bundle
username: gsapubftp-anonymous
password:

下载后的hg38的bwa_index文件夹内有以下文件:

代码语言:javascript
复制
kelly@DESKTOP-MRA1M1F:/mnt/f/kelly/bioTree/server/wesproject/hg38$ tree -h
.
├── [1.8G]  1000G_phase1.snps.high_confidence.hg38.vcf.gz
├── [2.0M]  1000G_phase1.snps.high_confidence.hg38.vcf.gz.tbi
├── [4.0K]  bwa_index
│   ├── [ 20K]  gatk_hg38.amb
│   ├── [445K]  gatk_hg38.ann
│   ├── [3.0G]  gatk_hg38.bwt
│   ├── [767M]  gatk_hg38.pac
│   ├── [1.5G]  gatk_hg38.sa
│   ├── [6.2K]  hg38.bwa_index.log
│   ├── [   0]  index.129248.err
│   ├── [   0]  index.129248.out
│   └── [ 566]  run.sh
├── [4.0K]  cnv
│   ├── [5.8M]  hg38.chr.bed
│   ├── [6.7M]  list.interval_list
│   ├── [ 843]  readme.txt
│   └── [5.7M]  targets.preprocessed.interval.list
├── [3.2G]  dbsnp_146.hg38.vcf.gz
├── [3.0M]  dbsnp_146.hg38.vcf.gz.tbi
├── [ 59M]  hapmap_3.3.hg38.vcf.gz
├── [1.5M]  hapmap_3.3.hg38.vcf.gz.tbi
├── [568K]  Homo_sapiens_assembly38.dict
├── [3.0G]  Homo_sapiens_assembly38.fasta
├── [157K]  Homo_sapiens_assembly38.fasta.fai
├── [ 20M]  Mills_and_1000G_gold_standard.indels.hg38.vcf.gz
├── [1.4M]  Mills_and_1000G_gold_standard.indels.hg38.vcf.gz.tbi
└── [1018]  run.sh
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2019.05.30 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 1创建wes环境
  • 2 所有需要的软件安装在wes环境下
  • 3 gatk加入环境变量
  • 4 熟悉参考基因组及必备的数据库
    • 4.1 参考基因组下载
    领券
    问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档