你会用到的网站:

写在前面:

之前mac不小心升级了一下java,然后igv就不能用了,要写教程必须降级java首先,看官方说明,需要安装Java -8,9以上版本不支持。我的mac不知道什么时候更新到了java 10,按说可以向下兼容,但是事与愿违,igv不能正常使用了。

需要降级Java,mac用户可以直接参考,windows可以试下直接下载安装IGV:先删除原来的javaterminal打开终端,复制粘贴一下三条命令:sudo rm -fr /Library/Internet\ Plug-Ins/JavaAppletPlugin.plugin

sudo rm -fr /Library/PreferencesPanes/JavaControlPanel.prefPane

sudo rm -fr ~/Library/Application\ Support/Java⚠️:不要通过/usr/bin 删除 Java 工具来卸载 Java。此目录是系统软件的一部分,下次对操作系统执行更新时,Apple 会重置所有更改。

finder中进入 /Library/Java/JavaVirtualMachines,然后删除之前的jdk.版本号

【Windows用户下载,解压后,点击igv.bat文件即可启动;如果启动失败,用记事本打开并编辑igv.bat文件,在文件的最后新起一行输入pause,保存后,再尝试打开,就可以在Windows下的命令行界面(cmd命令提示符)看到错误信息,再根据信息提示去解决问题;不过一般问题不大】

IGV安装

正文开始:

好啦,问题解决啦,开始正式IGV介绍!

什么是IGV

它是一款本地的探索基因组数据的可视化浏览器,有多个系统版本,支持多种不同类型的输入格式,包括芯片测序、二代测序、基因组注释文件等。推荐使用BAM与SAM格式,主要格式见下表数据来源文件格式序列比对SAM/BAM

显示覆盖率TDF

拷贝数SNP、CN

基因表达GCT、RES

基因注释GFF3/GTF、BED

突变数据MUT

追踪参考基因组覆盖度、测序深度(UCSC)WIG、BW

一睹IGV

每次打开会自动加载hg19.fa文件,也就是人类基因组,一会进入主界面

主界面

自己构建基因组信息

这里我会举一个昆虫中一种——棉铃虫,这个基因组是17年2月更新在NCBI伤的,属于小众物种,IGV并没有收录。正好拿来练手,当然如果你研究的领域也有基因组被测出来,也可以试一试【注意:在提交基因组文件到IGV之前,要先构建索引】

这些工作都可以在本地进行,只需要打开你本地的git_bash或者putty/xshell或者terminal,解压缩基因组文件=》下载samtools(推荐用conda管理)=>构建索引samtools faidx genome.fasta=>IGV中 输入fasta文件路径=》提供注释文件(可以是组装基因组预测的基因注释文件,也可以是拼接转录组用的gtf文件)=〉其他选项可以忽略=》点击OK推弹出一个框让你输入存储路径

自定义基因组

查看注释文件这里以人类基因组注释文件为例,下载gtf到电脑

下载完不要急着导入,需要先构建索引

导入注释文件

然后会生成gff3.idx或者gtf.idx文件,说明构建了索引,接着导入File -> Load from file,选择sorted的注释文件

查找基因

查看bam文件

我这里准备的bam文件大小是2.8G,是由人类转录组测序数据得到的,准备的参考基因组是hg19,注释文件是gencode.v28lift37.annotation.sorted.gff3bam文件在导入前,要先使用samtools进行sort和index,samtools sort test.bam test.sort``samtools index test.sort.bam,生成一个后缀为“.fai”的文件,它根据文件名自动和.bam关联, 另外这两个文件要在一个文件夹下,最后将bam导入IGV中

载入bam后,默认会出现两个track(翻译的话,可以理解为不同的轨道,显示不同的信息)Coverage track和Alignment track。

载入bam后

另外基因组信息也可以有collapsed、expanded、squished三种展示形式

基因组信息查看Coverage track

它的意思是显示比对文件的覆盖率和测序深度。横坐标是基因组上的位置,纵坐标是该位置的测序深度。【鼠标放在每一个位点都会显示一个小方框,其中的的内容就是显示总共有多少reads在这个位置,每个碱基各是什么】

点右上角+放大reads可视化窗口后,track会以灰色的条形图来显示每个位点的测序深度。如果某一个核苷酸与参考序列相比,有超过20%的reads是不同的,条形图会显示不同的颜色

Coverage track

显示特定位点变异

关于上图中的右键菜单,解释如下功能含义Rename Track更改track名

Change Track Color更改背景色,比如把Coverage Track灰色变紫色

Change Track Height改变每一个track的高度

Change Font Size改变IGV最左侧字体大小

Set Data Range覆盖深度的范围设置

Log scale用对数尺度作图

AutoScale是否自动缩放比较多个基因

比较多个基因

这里看到有许多颜色,这些颜色是根据定义不同比对类型而不同,

不同比对类型颜色也不同查找结构变异

灰色:与参看基因组可以比对的reads

紫色I:插入(鼠标查看插入的碱基信息)

黑色横线:——缺失

欢迎关注我们的公众号~_~

我们是两个农转生信的小硕,打造生信星球,想让它成为一个不拽术语、通俗易懂的生信知识平台。需要帮助或提出意见请后台留言或发送邮件到Bioplanet520@outlook.com

Welcome to our bioinfoplanet!

打开 igv java_必备可视化Integrative Genomic Viewer(IGV)相关推荐

  1. JSON必备工具之Json Viewer

    玩json的朋友有时候看那一些json数据很郁闷,很杂乱,要是有一款工具能像XML一样那样的格式显示就好了,现在向大家推荐Json Viewer 下载地址:http://download.csdn.n ...

  2. 人口热力图java_数据可视化API之热力图实现

    前言 数据可视化API(Web),是基于腾讯位置服务JavaScript API GL实现的专业地理空间数据可视化渲染引擎. 通过这套API,可以实现轨迹数据.坐标点数据.热力.迁徙.航线等空间数据的 ...

  3. linux中IGV的运行,[Linux] IGV使用说明

    做完mapping之后,需要知道mapped reads在基因组上的分布,所以尝试在ubuntu系统中使用IGV软件,简单记录使用过程,如下: 1.下载IGV软件 去IGV官网: http://www ...

  4. UCSC寻找基因位置并用IGV画出基因转录示意图

    寻找基因位置并寻找出基因转录示意图 我们这里教大家采用两种方法进行基因转录示意图的查找,一种是直接在UCSC里面进行查找,不过不够方便,并且背景有条纹.另一种是下载基因注释文件到IGV基因可视化软件里 ...

  5. ATAC-seq【Harvard FAS Informatics】

    ATAC-seq数据质量评估注意 ENCODE的ATACseq数据标准. Uniform Processing Pipeline Restrictions The read length prior ...

  6. 转录组分析综述A survey of best practices for RNA-seq data analysis

    转录组分析综述 转录组 文献解读 Trinity cufflinks 转录组研究综述文章解读 今天介绍下小编最近阅读的关于RNA-seq分析的文章,文章发在Genome Biology 上的A sur ...

  7. Lecture 3——DNA-seq-1

    本文图片来自于学习视频--新一代测序技术数据分析第三讲 DNA-seq Review Alignment srategies Smith-Waterman(speed too slow to use) ...

  8. tcga数据下载_肿瘤研究不能不知道的TCGA数据库挖掘工具大全,TCGA再也不愁

    TCGA数据库的挖掘工具层出不穷,从数据下载到数据挖掘,这里小编给大家整理一份官网的数据挖掘工具大全: 1. http://www.cancerimagingarchive.net/ The Canc ...

  9. Google Docs Viewer在线打开阅读PDF、Word、PPT、Excel等

    Google Docs Viewer在线打开阅读PDF.Word.PPT.Excel等 免费网站 2009年10月19日 17:57 评论» 朋友发来了一个PDF.Word.PPT.Excel文件或者 ...

最新文章

  1. 解决Kali Linux没有声音
  2. 买股票的最佳时机(六种题解dp)
  3. numpy如何直接在数组上进行添加_NumPy 论文登上了 Nature!
  4. 操作篇 了解学习NAT实验(关于静态NAT、动态NAT、EasyIP、端口映射的配置方法))
  5. Ant Design 3.15.0 发布,企业级 UI 设计语言和 React 实现
  6. 股票的科创板,新三板,创业板到底哪个能让你赚钱
  7. [基础]全局组、域本地组、通用组的区别与联系
  8. Lua5.2 请求 luasocket 相关模块时的 multiple-lua-vms-detected
  9. 离散数学复习--第二章:一阶逻辑
  10. 6120:学长的沙包
  11. java销毁servlet_servlet一般什么时候销毁
  12. 给技术经理找了几款Docker开源镜像仓库,为什么经理选中了Sonatype Nexus(下)
  13. 转载 | 深度学习中的遥感影像数据集
  14. 自定义注解,实现业务处理
  15. java代码耗尽内存_windows server 2008 环境下,运行java程序,内存耗尽问题
  16. 自建低成本代码托管与文档管理平台经验分享
  17. SQL存储过程根据,循环取出数据
  18. python随机生成生日测试生日悖论
  19. liunx基础及系统移植
  20. 有哪些靠谱的服务器安全软件?

热门文章

  1. c# 错误 两个输出文件名解析为同一个输出路径
  2. JQ基础练习---图片划过变暗
  3. 巧用“傍术”选择陈列点
  4. 当某个快捷键不能用时很可能是热键冲突
  5. WPF学习笔记(三)
  6. java tomcat jms_JavaWeb之使用Tomcat、JNDI与ActiveMQ实现JMS消息通信服务
  7. 解决IntelliJ创建Maven项目一直显Loading archetype list
  8. Android kotlin实现底部导航栏
  9. jdk和maven配置
  10. FileZilla Server 搭建FTP服务器