使用IGV看序列比对情况

本文从以下五个方面介绍了可视化序列比对数据和相关的tracks:

  • 文件格式:推荐的是BAM/SAM,其他格式,并且需要进行sorting&indexing
  • Read 覆盖率:整体视图,默认的覆盖率视图,和扩展覆盖率视图
  • 序列比对track:颜色、透明度、插入、缺失和排序
  • PE序列比对:将reads以pairs形式和颜色来区分,同时可以分为几个屏幕看。

文件格式

IGV推荐使用格式是:BAM以及SAM格式。

除了BAM,GOBY、VCF、PSL、BED、TDF等格式IGV也支持。

Sort和Index

BAM文件在载入IGV前,需要进行sort 和index。Index会生成一个以“.fai”结尾的辅助文件,这个文件会根据文件名自动关联序列比对数据(.bam),导入IGV中。

  • 序列比对数据需要以.bam扩展名结尾;
  • 在进行index时,文件名必须一致,index文件也必须在同一个文件夹下。 比如:test_xyz.bam文件的index文件名应该是test_xyz.bam.bai,或者是test_xyz.bai。
  • 这两步骤可以用samtools或Picard软件进行。

Tracks

载入bam文件后会产生3个相关的tracks:

  • Alignment track :显示每个的reads的比对情况
  • Coverage track:显示覆盖率和测序深度
  • Splice Junction Track:提供一个可选的横跨剪切位点(spanning splice junctions)的reads视图。

一般情况下,前两个tracks会自动出现。这些设置可以通过右键进行修改。

Coverage Track

默认情况下,IGV能动态计算和显示比对文件的覆盖率和测序深度。当IGV窗口放大到reads 可视化阈值大小(默认为30KB)时,这个track会以灰色条形图显示每个位点的测序深度。如果某核苷酸与参考序列不同(超过20%reads)时,IGV会标出不同的颜色。

即:A→绿色;C→蓝色;G→橙色;T→红色。

  • 将鼠标悬停在你需要查看的位点处可以看到详细的信息,右键可以复制。

覆盖率数据(TDF)

可用igvtools将BAM文件转化为TDF格式,这个文件是专门显示覆盖率

TDF文件是BAD文件的精简版,当只需要看覆盖率数据时,可用igvtools工具进行转换;方便快速查看。

Alignment Track

当IGV窗口放大到30KB(默认)时,就会显示出各个reads的情况。

  • 当窗口放大到一定程度时,IGV的窗口中心会出现一条线,再继续放大,中心线会变为两条,刚好可以框住一个碱基。

结构变异(Structural Variants)

IGV使用颜色和其他的标记来显示变异:SNP、SV、异倍体(aneuploidy)。

IGV使用reads透明度来表示其质量。

  • 灰色:和参考基因能比对上的reads;
  • 紫色 I:插入;(鼠标悬空在此处能显示插入的碱基信息)
  • 黑色横线(—):缺失;
  • 注意:那些透明或者白色的reads有着亮灰色边框的其比对质量(MQ)为0.

Paired-end 比对

对于PE比对,igv可通过右键选择将reads以成对的方式展示。(如下图)

用户可以将需要标记的reads用紫色标记(1),而不同于预期的将会同(3)一样进行标记。

  • Ctrl+鼠标单击(Mac:command+click)将paired reads用相同的颜色标记其轮廓,每对颜色都是不一样的。注:黑色外框的reads则该reads没有mate。
  • 而Ctrl+鼠标单击(Mac:command+click)任意的一个read都可以去掉这个轮廓。
  • 单击右键,选择Go to mate region可以找到PE中的另一条read。
  • 如果这对reads的插入片段大小较大,则标记一条read时,其另一条read将不会被标记。
  • 停留在一个read上或者单击它,可以看到这条read的信息,包括它paired的另一条read。

在PE比对中,还可以通过右击一个read,选择View mate region in split screen,分屏查看paired mate。(注:若这个read没有mapped mate,则这个选项为灰色)。双击窗口顶端的位点栏,即可恢复正常的单屏视图。

Interpreting Color by Insert size

IGV用颜色来标记异常的插入片段大小的reads,这个只用DNA比对(不适用于RNA-Seq)。

默认的颜色标记规则:如图

  • 红色:插入片段大小大于预期值(这里可能有缺失)
  • 蓝色:插入片段大小小于预期值(这里可能存在有插入)
  • 不同颜色显示PE reads中其mate可在其他染色体上找到(此处可能为染色体间重组)

Interpreting Color by Pair Orientation

  • inversions
  • duplications
  • Translocations

本文简单介绍了使用igv查看序列比对的情况,还有其他模式如Bisulfite Mode,详细可看原文网页(http://software.broadinstitute.org/software/igv/book/export/html/37)。期待大家的交流学习哇。

原文发布于微信公众号 - 生信技能树(biotrainee)

原文发表时间:2017-07-24

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏web编程技术分享

flex布局方法详解之flex-direction

1344
来自专栏小李刀刀的专栏

[译]CSS边框实现“无图化”设计

原文:http://www.yuiblog.com/blog/2011/03/08/css-border-tricks-with-collapsed-boxes...

3478
来自专栏编程

一行JS代码实现翻翻卡效果

本文章是由源码时代H5前端讲师分享:一行JS代码实现翻翻卡效果,希望对大家有一定帮助。 实现步骤 准备2张图片,一个在上一个在下 上方的图片不动,下方的图片先翻...

1668
来自专栏hightopo

实践 HTML5 的 CSS3 Media Queries

1082
来自专栏老司机的简书

CoreText实现图文混排

也好久没来写博客了,主要是最近也工作了,手头的事有点多,一时间也就断了,闲下来了我就来补博客了,刚好最近也做了很多东西,放在这里也算给自己做个笔记吧。

612
来自专栏HT

基于 CSS3 Media Queries 的 HTML5 应用

先来介绍下 media,确切的说应该是 CSS media queries(CSS 媒体查询),媒体查询包含了一个媒体类型和至少一个使用如宽度、高度和颜色等媒体...

1655
来自专栏Material Design组件

Material Design — App bars: topApp bars: top

1546
来自专栏python爬虫实战之路

极验验证码破解之selenium

这一篇写完很久了,因为识别率一直很低,没办法拿出来见大家,所以一直隐藏着,今天终于可以拿出来见见阳光了。

754
来自专栏IMWeb前端团队

移动端重构实战系列3——各种等分

本文作者:IMWeb 结一 原文出处:IMWeb社区 未经同意,禁止转载 ”本系列教程为实战教程,是本人移动端重构经验及思想的一次总结,也是对sand...

1.1K7
来自专栏偏前端工程师的驿站

CSS魔法堂:你真的理解z-index吗?

一、前言                                 假如只是开发简单的弹窗效果,懂得通过z-index来调整元素间的层叠关系就够了。但要将...

1695

扫描关注云+社区