生信技能树 WES分析教程学习(1)conda安装软件,配置环境
1.首先添加conda的镜像
本文后续的分析主要参考的是这篇文章,主要是想用E.coil的数据进行后续的分析
https://zhuanlan.zhihu.com/p/33891718
但是再下载SRA数据的时候一直没下载下来
因为后续下载到了GATK上bundel上的100G的数据,所以直接俄按照jimmy后续的分析操作进行了。
conda的软件源有很多,其中专门收录了生物信息学软件的软件源,亦即bioconda。所以bioconda仅仅是conda的软件源之一,与bioconductor之于CRAN,bioperl之于CPAN是类似的。所以掌握bioconda,事实上也就是掌握conda,反之亦然。
常用的镜像有,清华大学的镜像,中科大的镜像和兰州大学的镜像等。
安装Anaconda或者miniconda,如若没有root权限,建议安装后者,后者包轻便,前者包含自带的一个python,安装过程如下:
#下载包
wget "https://repo.continuum.io/archive/Anaconda3-5.0.1-Linux-x86_64.sh"
#安装程序
bash Anaconda3-5.0.1-Linux-x86_64.sh
#添加环境变量
vim /etc/bashrc
export PATH=/tools/anaconda3/bin:$PATH ##是你安装Anaconda的路径
配置bioconda
如果root身份运行则添加到全局,如果是用户身份运行则会修改用户主目录下的.condarc文件,与bashrc, bashprofile等类似
conda config --add channels conda-forge
conda config --add channels defaults
conda config --add channels r
conda config --add channels bioconda
以上,无论是conda默认的软件源还是bioconda软件源都是国外的,速度非常慢,
所以需要增加国内软件源,同时bioconda已经有清华,中科大两个国内镜像,也添加进去。
conda config --add channels https://mirrors.tuna.tsinghua.edu.cn/anaconda/pkgs/free/
conda config --set show_channel_urls yes
conda config --add channels https://mirrors.tuna.tsinghua.edu.cn/anaconda/cloud/conda-forge/
conda config --add channels https://mirrors.tuna.tsinghua.edu.cn/anaconda/cloud/msys2/
#bioconda
conda config --add channels https://mirrors.tuna.tsinghua.edu.cn/anaconda/pkgs/free
conda config --set show_channel_urls yes
#在广州应该使用中科大的bioconda源(清华的源可能远水救不了近火...)
conda config --add channels https://mirrors.ustc.edu.cn/anaconda/cloud/bioconda/
conda config --set show_channel_urls yes
#查看目前conda软件源情况
conda info
gatk数据下载
先在目录下新建文件夹 input output 和bin
然后在input数据里下载E.coli K12的参考基因组序列,进行数据测试
可以在NCBI的FTP上进行下载
ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/005/845/GCF_000005845.2_ASM584v2/GCF_000005845.2_ASM584v2_genomic.fna.gz wget ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/005/845/GCF_000005845.2_ASM584v2/GCF_000005845.2_ASM584v2_genomic.fna.gz
得到的时Ecoil K12的全基因组序列
如下图所示
然后解压缩,并命名为 E.coli_K12_MG1655.fa
gzip -dc GCF_000005845.2_ASM584v2_genomic.fna.gz > E.coli_K12_MG1655.fa
head E.coli_K12_MG1655.fa
作者接着使用samtools进行构建参考基因组的索引
因此在Biosoft文件夹下载samtools
mkdir ~/Biosoft/samtools
cd /root/Biosoft/samtools
生信技能树 WES分析教程学习(1)conda安装软件,配置环境相关推荐
- R:生信技能树学习笔记一
生信技能树小破站:R应该这样学1-4 1.查看已经安装的包的地址 .libPaths() 2.怎么查看函数用法 #在RStudio的右下角窗口的help可以看到 ?函数名 3.三个有用的函数 1.he ...
- R:生信技能树学习笔记二
生信技能树小破站:R应该这样学5-7 1.热图 rm(list=ls()) library(pheatmap) a1=rnorm(100) dim(a1)=c(5,20) #设置维度 pheatmap ...
- CancerSubtypes包的介绍(根据生信技能树Jimmy老师分享的乳腺癌分子分型包资料整理)
CancerSubtypes包的介绍(根据生信技能树Jimmy老师分享的乳腺癌分子分型包资料整理,感谢Jimmy老师!) 1. 引言 2. 数据处理 2.1 基本处理 2.1.1 通过检查数据分布来分 ...
- ProTICS包的介绍(根据生信技能树Jimmy老师分享的乳腺癌分子分型包资料整理)
ProTICS包的介绍(根据生信技能树Jimmy老师分享的乳腺癌分子分型包资料整理,感谢Jimmy老师!) 1.设置环境 2.Part1的结果 3.Part2的结果 4.Part3的结果 5.相关函数 ...
- python基因差异分析_R语言之生信②差异基因分析2
目录 R语言之生信②差异基因分析2 样品的无监督聚类 检查基因表达分析最重要的探索性策略之一是多维定标(MDS)图或类似的图.该图以无监督的方式显示了样本之间的相似性和不相似性,以便人们可以了解在进行 ...
- ionic3 教程(一)安装和配置 1
链接: ionic3教程(一)安装和配置 ionic3教程(二)登录页制作 ionic3教程(三)设置页制作 ionic3教程(四)安卓硬件返回键处理 ionic3 教程(五)基本的网络请求 学习本教 ...
- ZED相机学习笔记1——安装与配置(Win10 + Python)
系列文章目录 ZED相机学习笔记1--安装与配置(Win10 + Python) 文章目录 系列文章目录 前言 一.ZED2 相机 二.配置ZED相机环境 1.安装CUDA 2.安装ZED-SDK 3 ...
- Xamarin.Forms教程下载安装JDK配置环境变量
Xamarin.Forms教程下载安装JDK配置环境变量 Xamarin.Form环境配置下载安装JDK JDK是编程Java程序必须的软件.也许有人会问我们用的C#为什么还有Java呢?这是因为我们 ...
- mesos 学习笔记-- mesos安装和配置
2019独角兽企业重金招聘Python工程师标准>>> mesos 学习笔记-- mesos安装和配置 博客分类: 架构 mesos 参考资料: 官方文档:http://mesos. ...
- PIC单片机入门教程(二)—— 安装集成开发环境(MPLAB X IDE)
PIC单片机入门教程(二)-- 安装集成开发环境(MPLAB X IDE) 教程中使用的电脑运行Windows 10 专业版 64位系统 1.下载 MPLAB X IDE v4.05 历史版本:htt ...
最新文章
- Survey | 生物医学文本挖掘最新进展
- R语言使用psych包的fa函数对指定数据集进行因子分析(输入数据为相关性矩阵)、使用rotate参数指定进行斜交旋转提取因子、使用factor.plot函数可视化斜交旋转因子分析、并解读可视化图形
- 备忘录同步android ios9,IOS备忘录
- 希尔排序是一种稳定的排序算法_全面解析十大排序算法之四:希尔排序
- 回调函数到底是怎么一回事呢
- 关于Oracle undostat中的2012和ORA-01555问题的自我解答
- 粉丝回馈 | 用30天换来了这辈子独一无二的纪念
- pageable设置size_总结SpringBoot使用Pageable实现分页源码
- 利用DllPlugin分割你的第三方库
- 维护LINQ to SQL多对多表间关系[转]
- 第15章-使用远程服务
- word中插入对号、错号符号
- 信息安全毕设论文选题推荐
- Office更新了那么多代,哪个版本好看呢?
- 走向.NET架构设计—第三章—分层设计,初涉架构(后篇)
- java 随机金额_java_微信随机生成红包金额算法java版,最近几年玩得最疯狂的应该是 - phpStudy...
- 不用代码免费采集网站内容_看一遍就能上手的火车头采集图文全流程教程
- 强跟踪ukf matlab,一种改进的强跟踪UKF算法及其在SINS大方位失准角初始对准中的应用...
- 证明:因为矩阵A非0,所以AA'也非零
- 消息传递框架MPNN: Neural Message Passing for Quantum Chemistry