IGV基因组浏览器可视化高通量测序数据
IGV是本地浏览测序数据功能最为强大的基因组浏览器,支持多种不同类型的输入格式和不同的显示方式,如峰图、线图、柱状图、Sashimi-plot。同时还可以配合bedtools使用。
新版的Windows IGV分发包中包含了Java运行环境,点击链接http://data.broadinstitute.org/igv/projects/downloads/IGV_Win_2.3.93.zip即可下载,解压后,点击igv.bat文件即可启动。
若启动失败,使用记事本打开并编辑igv.bat文件,在文件的最后新起一行输入pause,保存后,再尝试打开,就可以在Windows下的命令行界面(黑色背景的小框中)看到错误信息,根据信息提示去解决问题。
IGV具体使用如下:
如果基因组未存在于当前列表中,则选择更多从Broad IGV获取对应的基因组信息。若所研究的物种没有被Broad IGV收录,则需要自己构建基因组信息。具体参见前面发过的文章 (在公众号后台输入 IGV基因组索引 即可获取)。
IGV展示所用的bw文件可以直接从bam文件得来。下面列出几个工具和相应的代码辅助转换。考虑到不同样品的测序深度不同,需要根据样品的测序量做一个标准化处理,建议使用deepTools的bamCoverage命令。
IGV是一个界面可操作软件,学习的方式就是每个菜单点开看看,点选下看下效果,比如修改下Track的颜色,Track的高度,纵轴的坐标范围Auto scale, Log scale 或 Data range等方便比较。
如果想比较多个基因的表达峰图,可以先在IGV中添加相应的基因列表,然后针对特定的列表选择view,就可以在一个屏幕显示所有加入的基因。
如果我们要查询的区域很多,可以配合bedtools igv -i input.bed,生成一个脚本,批量捕获多个区域的比对信息或coverage信息。使用Tools - Run Batch Script - 选择bedtools igv的输出结果,即可以看到动画般的截屏界面。完成后,就可以一张张的浏览图片或者传给没有安装IGV的老师或朋友查看。
构建自己的基因组索引
- 在IGV的菜单栏点击
Genomes
>Create .genome File
,会弹出一个窗口填写需要的信息。 - 输入一个ID和基因组的描述信息,这个ID等同于常见的mm10,会作为构建的基因组索引的名字。
- 使用
samtools faidx genome.fasta
索引基因组文件,然后输入基因组FASTA格式文件的路径。 - 提供基因注释的
GTF
文件,这个文件可以是拼装基因组预测的基因,也可以是cufflinks
等拼装的转录组的GTF文件。 - 其它的选项可以忽略。
- 点击
Save
,选择基因组索引文件*.genome
的存储路径。
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- Linux中修改mysql默认编码的方法步骤
- ubuntu16.04搭建nfs服务的方法
- Ubuntu16.04搭建php5.6Web服务器环境
- Linux上通过SSH挂载远程文件系统方法详解
- Linux上创建、列出、删除Docker容器方法总结
- Centos 7.4服务器时间同步配置方法【基于NTP服务】
- PowerBI 通用万能日历模板,想干嘛就干嘛
- Linux持续集成自动化安装Maven的方法
- 在CentOS上安装搭建PHP+Apache+Mysql的服务器环境
- Centos7.4服务器安装apache及安装过程出现的问题解决方法
- ubuntu系统修改时区和时间的方法
- CentOS7 LNMP+phpmyadmin环境搭建 第一篇虚拟机及centos7安装
- flutter 输入框组件TextField的实现代码
- Mac通过不同终端SSH连接远程服务器的讲解
- Linux关于透明大页的使用与禁用介绍