可视化工具之 IGV 使用方法
整合基因组浏览器(IGV)是一种高性能的可视化工具,用来交互式地探索大型综合基因组数据。它支持各种数据类型,包括array-based的和下一代测序的数据和基因注释。IGV这个工具很牛,发了NB:
使用IGV的基本思路
安装和配置
要用linux上的IGV(首先需要去官网下载一个Linux下的IGV),然后需要一个远程控制工具:
这里使用 VNC,网上下载破解版,安装(输入注册码即可破解,文章底部有注册码)。
运行IGV
开启 linux 的 VNCserver,开启了之后,linux上会显示节点和端口,用于Windows端的连接。
vncserver -geometry 1366x768 #调节分辨率
开启服务后,Linux会提示可以登录的节点,登录VNC需要你的ip、端口及VNC密码。
运行 VNC 客户端(VNC viewer)里的 IGV.sh 开启可视化程序(打开Windows的VNC会自动跳转到Linux下的IGV安装目录,双击IGV.sh即可打开)。
1.载入参考基因组数据
可以点击选择一个,也可以通过菜单栏Genomes载入文件。
2.载入排序后的bam文件
通过菜单栏的File,load from file。。。
IGV可视化数据准备
不是什么数据都可以拿IGV看的,参考基因组必须为FASTA格式;
IGV只是负责将比对结果可视化,并没有比对过程,所以不能直接载入reads;
需要将待比对的reads与前面指定的参考基因组用bwa进行比对;
比对后的sam文件也不能直接载入(麻烦),要转bam;bam排序;bam建索引(可以一步完成);
命令流程如下:
BWA出来的sam文件无法直接被 IGV 利用的,必须经过 samtools 处理后才能被 IGV 显示出来。
samtools view -@ -bS -F ./contigs_sequence_align_to_public_genome.sam > ./contigs_sequence_align_to_public_genome.bam
samtools sort -@ ./contigs_sequence_align_to_public_genome.bam contigs_sequence_align_to_public_genome.sorted
samtools index contigs_sequence_align_to_public_genome.sorted.bam contigs_sequence_align_to_public_genome.sorted.bai
samtools depth contigs_sequence_align_to_public_genome.sorted.bam >depth_reads.txt
wc -l depth_reads.txt > Coverage-aln_reads.txt
更多的使用细节参照 samtools 使用方法。
IGV 需要设置内存,在看大型基因组时很吃内存.
示例:
下面显示的是组装异常区,这是手动连接的地方,Pacbio全部不能顺利跨过去,说明组装有问题。
以下是单细胞转录组比对可视化:
IGV基本操作
如何去掉右边自动弹出的框框?
蓝色框起来的reads是啥?
…
IGV用户指南
目录:
- 用户界面
- 开始使用IGV
- 导航
- 加载一个基因组
- 外部控制IGV
- 查看参考基因组
- 加载数据和属性
- 查看数据
- 查看比对(核心)
- 查看Variants
- 基因列表视图
- 感兴趣的区域
- 样本的属性
- 排序、分组和过滤
- 保存和恢复会话
- 服务器配置
- Motif查找器
- igvtools
- BLAT 搜索
1.用户界面
1.1 Main Window(主窗口)
1.tool bar(工具栏),menu bar(菜单栏),pop-up menus(弹出式菜单)。
2.染色体上的红色盒子表示显示这部分染色体,显示完整染色体是红框会消失。
3.尺度显示了染色体的可见部分,刻度线显示了染色体的位置,跨度列表显示了当前显示的碱基的数量。
4.IGV在水平行显示的数据称为tracks。通常,每个tracks代表一个样本或实验。这个例子展示了甲基化、基因表达、拷贝数,LOH和突变数据。
5.IGV也显示某些特性,比如在tracks中的基因。默认情况下,IGV在一个面板显示数据,在另一个面板显示数据特性。拖放一个track名称,将一个track从一个面板移动到另一个地方。
6.Track名称列在最左边面板。名字的易读性取决于 tracks的高度,例如,track越小,它的名字的可读性越小。
7.属性名称被列在顶部的属性面板。彩色块代表属性值,每个独特的值被都有一个独特的颜色。鼠标放在一个颜色块的附近来查看其属性值。
参考资料
Viewing Alignments(查看比对,核心文档)
IGV软件使用指南(博客)
软件安装:
RealVNC是一款免费的远程控制程序,它是VNC (Virtual Network Computing)众多操作平台版本中的一员,这里小编提供的RealVNC适用于windows操作平台,软件采用2048 位 RSA 服务器验证和 128 位 AES 会话加密术,安装此软件的电脑可在全球各地远程操控自己的电脑,所操作平台没有任何限制。
功能特色
1.Windows 验证 —— 允许访问 VNC 服务器以便控制基于标准的 Windows 用户证书, 无论其是本地的、NT 域或者基于 Active Directory。
2.单一登录 —— 允许在一个域内部访问 VNC 服务器而不需要重新输入用户的名称和口令。
3.简单但有效的部署 —— 用 VNC Tool,即 VNC Deployment Tool for Windows 安装和 管理 VNC Enterprise Edition 服务器。
4.桌面缩放 —— 通过精确的比例缩放,或动态缩放为任何大小。
5.Windows 防火墙集成 —— 使 VNC Server 可以更直接的部署。
6.文件传输 —— 与 Windows 剪贴板结合,允许通过 VNC 连接在服务器和查看器计算机之间 复制文件而不需要额外的配置。
最好的就是,与 VNC Enterprise 系列的其他产品一样,VNC Enterprise Edition for Windows是由原始 VNC team 设计、开发和维护的。
RealVNC注册码
3YHED-MNEHC-RMJT5-4UAAK-6A5HA
FBV9V-7Z3V9-MED3U-47SEU-85T3A
HA7MG-J8R3J-R3528-HC6P5-HTQ6A
B7SLM-7MAX5-B4M74-UTDBE-K5WFA
7736G-H4YBP-6ZCEU-AM3H4-QBTWA
NUQQ6-UP89V-9XGPU-GUYU9-U6VTA
可视化工具之 IGV 使用方法的更多相关文章
- java桌面应用开发可视化工具windowbuilder详细使用方法
http://blog.csdn.net/qq_28859405/article/details/52562131
- 99%的人都搞错了的java方法区存储内容,通过可视化工具HSDB和代码示例一次就弄明白了
https://zhuanlan.zhihu.com/p/269134063 番茄番茄我是西瓜 那是我日夜思念深深爱着的人啊~ 已关注 6 人赞同了该文章 前言 本篇是java内存区域管理系列教 ...
- 漫谈可视化Prefuse(五)---一款属于我自己的可视化工具
伴随着前期的基础积累,翻过API,读过一些Demo,总觉得自己已经摸透了Prefuse,小打小闹似乎已经无法满足内心膨胀的自己.还记得儿时看的<武状元苏乞儿>中降龙十八掌最后一张居然是空白 ...
- 可视化工具gephi源码探秘(二)---导入netbeans
在上篇<可视化工具gephi源码探秘(一)>中主要介绍了如何将gephi的源码导入myeclipse中遇到的一些问题,此篇接着上篇而来,主要讲解当下通过myeclipse导入gephi源码 ...
- 【MongoDB】2.可视化工具的安装和使用
首先:关于 能支持MongoDB新版本的可视化工具,争议不断,各人都有各人的支持. 因此之前选择安装的时候就安装了MongoDB 3.0.14版本的. 最终,确定使用Robomongo作为我第一个 ...
- mongodb下载、安装、配置服务启动、及可视化工具下载、使用
MongoDB: 1.下载地址:http://www.mongodb.org/downloads(32位还是64位自行选择).我下载的是:mongodb-win32-x86_64-3.2.4-sign ...
- MYSQL系列1_MySQL的安装,可视化工具的使用,以及建库建表等
大家都知道MYSQL是开源的数据库,现在MYSQL在企业中的使用也越来越多,本人之前用过SQL SERVER数据库,因业务需要和自己的兴趣想要学习MYSQL,对于MYSQL,本人还是新手,请大家多多指 ...
- mysq数据库管理工具navicat基本使用方法
navicat是mysql数据库的客户端查询管理工具,本文详细的介绍了该软件的基本使用方法 本文转自 http://hejiawangjava.iteye.com/blog/2245758 sql是操 ...
- 高效使用 Python 可视化工具 Matplotlib
Matplotlib是Python中最常用的可视化工具之一,可以非常方便地创建海量类型的2D图表和一些基本的3D图表.本文主要介绍了在学习Matplotlib时面临的一些挑战,为什么要使用Matplo ...
随机推荐
- 湘潭 A simple problem
A simple problem Accepted : 30 Submit : 303 Time Limit : 15000 MS Memory Limit : 655360 KB Probl ...
- Android HandlerThread 完全解析
HandlerThread本质上是一个线程类,它继承了Thread: HandlerThread有自己的内部Looper对象,可以进行looper循环: 通过获取HandlerThread的loope ...
- Mysql-学习笔记(==》函数的建立与使用 十)
函数的建立与使用 USE db;SELECT sname,sscore,CASE WHEN sscore>=90 THEN '优秀'WHEN sscore>=70 THEN '良好'WHE ...
- java写的web服务器
经常用Tomcat,不知道的以为Tomcat很牛,其实Tomcat就是用java写的,Tomcat对jsp的支持做的很好,那么今天我们用java来写一个web服务器 //首先得到一个server, S ...
- HDU4549 M斐波那契数列 矩阵快速幂+欧拉函数+欧拉定理
M斐波那契数列 Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)Total Sub ...
- BQ27510 电量计的校准 的 C语言实现
点击打开链接 根据TI官方MSP430平台移植修改过来的,在Omap37xx(wince)平台测试,理论上和平台无关,伸手党,赶紧复制粘贴代码吧.如果这篇文章帮助了你,给了好评也无妨. [cpp] v ...
- 在JSP页面中输出字符" * "组成的金字塔
<%@ page language="java" contentType="text/html; charset=UTF-8" pageEncoding= ...
- Cheatsheet: 2014 01.01 ~ 01.14
.NET 15 reasons why I can't work without JetBrains ReSharper Web Web scraping with Node.js Koa.js : ...
- FreeSWITCH的传真发送
详细的学习请参考:https://wiki.freeswitch.org/wiki/Mod_spandsp 我只说说WIKI上没有直接写出来却又很实用的东西. 一.传真收发 环境大致如下: FreeS ...
- 机器学习——SVM详解(标准形式,对偶形式,Kernel及Soft Margin)
(写在前面:机器学习入行快2年了,多多少少用过一些算法,但由于敲公式太过浪费时间,所以一直搁置了开一个机器学习系列的博客.但是现在毕竟是电子化的时代,也不可能每时每刻都带着自己的记事本.如果可以掏出手 ...