原文地址:https://www.jianshu.com/p/f101c9238fc5
https://www.jianshu.com/p/1b9426d0f9f4
下载参考基因组
登陆网站http://genome.ucsc.edu/index.html下载hg19参考基因组
image.png
image.png
image.png
image.png
hg下载前是900M 左右,解压后是3G大小
这里需要在终端解压一下,普通解压方法不好用
tar -zxvf chromFa.tar.gz
#解压后可以发现,参考序列是按照染色体号分开列出的,我们还需要把所有的序列写入到一个文件中。
cat *.fa > hg19.fa
#最后删除其他无用的文件
rm chr*.fa
下载注释文件
官网:http://www.gencodegenes.org
image.png
image.png
# GTF格式主要是用来描述基因的注释axel ftp://ftp.sanger.ac.uk/pub/gencode/Gencode_human/release_27/GRCh37_mapping/gencode.v27lift37.annotation.gtf.gz# GFF文件是一种用来描述基因组特征的文件,现在我们所使用的大部分都是第三版)(GFF3)axel ftp://ftp.sanger.ac.uk/pub/gencode/Gencode_human/release_27/GRCh37_mapping/gencode.v27lift37.annotation.gff3.gz# 解压并删除原来的文件gzip -d gencode.v27lift37.annotation.gtf.gzgzip -d gencode.v27lift37.annotation.gff3.gz
image.png
IGV软件的下载和安装
IGV软件全称:Intergrative Genomics Viewer 是一个高效的查看基因数据的可视化软件。
官网:http://software.broadinstitute.org/software/igv/home
安装igv 之前需要先安装java 8 以上版本
官网:https://java.com/en/download/mac_download.jsp
IGV软件的使用
窗口
主窗口布局:
tool bar(工具栏),menu bar(菜单栏),pop-up menus(弹出式菜单)
染色体上的红色盒子表示显示这部分染色体,显示完整染色体是红框会消失
尺度显示了染色体的可见部分,刻度线显示了染色体的位置,跨度列表显示了当前显示的碱基的数量
IGV在水平行显示的数据称为tracks。通常,每个tracks代表一个样本或实验。这个例子展示了甲基化、基因表达、拷贝数,LOH和突变数据
IGV也显示某些特性,比如在tracks中的基因。默认情况下,IGV在一个面板显示数据,在另一个面板显示数据特性。拖放一个track名称,将一个track从一个面板移动到另一个地方
Track名称列在最左边面板。名字的易读性取决于 tracks的高度,例如,track越小,它的名字的可读性越小
属性名称被列在顶部的属性面板。彩色块代表属性值,每个独特的值被都有一个独特的颜色。鼠标放在一个颜色块的附近来查看其属性值
导入参考基因组及注释信息,查看感兴趣基因的结构
导入前面的 hg19.fa
接着导入gtf文件,需要先sort,才能导入
把count 改成 sort,在input里面输入gff3文件(先解压)
在output里面就会生成sorted文件然后在file-load from files-输入 sorted 的文件
-现在可以导入sort 后的gtf 文件了。
高通量测序图解见下章分析
参考文献
http://www.jianshu.com/p/48b5a0972301(GTF/GFF文件的差异及其相互转换)
http://www.jianshu.com/p/3e545b9a3c68(hoptop)
http://fbb84b26.wiz03.com/share/s/3XK4IC0cm4CL22pU-r1HPcQQ1lZQRc2nKQhn2SthRW24I8CZ(greenhillman)