IGV-基因组浏览器-改造记录(二)

写在前面

还有两个小时....似乎还是有点紧张。想想,那就写写微信推文。
IGV改造系列,一直会继续,直到或许是我们觉得似乎没有改造的必要。
在前述的改造记录中(https://www.jianshu.com/p/7e90e454d487),我们主要是给IGV增加了一个ReadLen Track,这样就可以实时地查看任何区域的sRNA读段的丰度和分布情况,在鉴定PHAS和数据观察上,有较大的用处 。
但是很多时候,更多人关注的不是phasiRNA,而是miRNA,或者有一定二级结构的siRNA前体。
这个时候,查看IGV的同时,观察一些二级结果,显得尤为重要。

调整RegionOfInterest弹出菜单,增加RNAfold功能

直接从IGV内部启动RNAfold,只能修改源代码。大体我使用了三步完成这个事情:

  1. 增加RNAfold调用源码

  2. 修改菜单,并增加响应调用事件

  3. 考虑正反链(一个是SenseFold,一个是AntisenseFold)

具体的使用比较简单,选定了感兴趣的区域之后,鼠标右击即可,

可以看到,已经增加了两个新的选项。
其中,SenseFold,是直接截取感兴趣的区域序列,进行RNAfold
而AntiSenseFold,则是截取后的区域序列,先做反向互补,随后再进行折叠

弹出的RNAfold窗口有较好的交互性

鼠标双击二级结构,即可调出菜单

直接在文本框中输入感兴趣的序列子串,
即可直接高亮出对应的碱基

对于选中的碱基,可以直接修改其一些特征,比如填充颜色

当然,还有其他功能....感兴趣的可以自行摸索

IGV与RNAfold窗口的动态交互

IGV可以显示测序数据在基因组上的覆盖情况和相关的详细信息,而RNAfold窗口,则展示对应区域的二级结构,结合两者分析,可以更好地发现一些东西

大体的使用是,fold的时候

当鼠标位于IGV的某个碱基位置时,RNAfold窗口会实时高亮出对应的碱基

如此,我们可以更好的观察到
每个碱基所在的二级结构位置和该位点的覆盖深度

写在后面

更多的观察数据,可以发现更多的生物信息。
数据的敏感度很重要,而工具也很重要。
顺手的工具,常常对应适合的方法。
OK....提前祝大家新年快乐.....

(0)

相关推荐

  • 【直播】我的基因组73:在IGV看看indel是啥样子

    前面我们特意用scalpel软件来找indel,期待它会有一些出彩的表现,当然我还没来得及比较它找到的INDEL跟GATK等工具区别在哪里,不过我们先在IGV里面看看找出来的是什么吧. [直播]我的基 ...

  • m6A图文复现05-比对结果以及可视化

    下面是MeRIP-seq 图表复现笔记 在上一期中我们得到hisat2的比对结果,见:m6A图文复现03-测序数据去除rRNA序列并且比对到参考基因组 现在一起来看下比对情况,并且进行一些常规的可视化 ...

  • IGV基因组注释文件说明

    IGV基因组注释文件说明

  • 【软件介绍】IGV软件的安装和基本介绍

    [软件介绍]IGV软件的安装和基本介绍 - 目录 1. IGV 下载与安装 2. Java 安装及环境变量设置 3. IGV 基本介绍 3.1 IGV 界面布局 3.2 IGV 结果界面 3.3 序列 ...

  • Molecular Plant | 20000个公共拟南芥RNA-Seq的综合在线数据库(超实用、免费、可访问)

    编译:阿菲,编辑:夏甘草.江舜尧. 原创微文,欢迎转发转载. 导读 在过去的十年里,由于测序成本的持续下降和最新的易于使用的文库构建工具的开发,RNA测序(RNA-seq)已经超过微阵列成为基因表达谱 ...

  • 修改IGV基因组浏览器源码,做一个自己需要的浏览器

    写在前面 从图片来看,即小RNA的reads展示在染色体上,并用长度做可视化.这个需求,其实拿起任何一门语言,手写一个.svg就可以搞定了.只是读取read alignment数据想想比较辛苦.有比较 ...

  • IGV-基因组浏览器-改造记录(一)

    我一直都相信, 只要有足够长的时间,没有什么事情是我们做不到的. 我也明白, 人的生命是有限的,所以追求的东西应该尽可能少. 我只是没想到, 很多时候,不是你没时间,而是你身边的人,等不了.而事实上, ...

  • IGV-基因组浏览器-改造记录(三)

    写在前面 前面已经做了两次IGV浏览器的改造.在进行第二次改造后,我们会发现一个明显的需求.即,二级结构上同时附加每个碱基的测序深度(类似表达量),那么将能更好的帮助判断miRNA或者特殊的sRNA二 ...

  • IGV-基因组浏览器-改造记录(四)

    写在前面 前述,我们已经对IGV进行了超过三次改造,同时我也写了三个推送.写IGV改造系列推送的主要原因,事实是作为课题组成员使用改造的IGV时的参考教程. 在上一次RNAfold的窗口上,我们增加了 ...

  • IGV-基因组浏览器-改造记录(五)

    写在前面 IGV是目前广泛使用的基因组浏览器.但在前述的推文中,我们已经提到过,其并没有较多我们做小RNA数据分析或可视化所需的特性.为此,我们已经介绍了四处我们课题组对IGV修改.我的意思是,我们做 ...

  • IGV-基因组浏览器-改造记录(六)

    写在前面 经过了前面的五次IGV改造,我以为这个系列可能就到此为止了. 然而在近期的数据分析过程中,我们仍然发现了一个新的需求: 在展示sRNA reads的时候,只看某个长度的reads,比如植物的 ...

  • IGV-基因组浏览器-改造记录(七)

    写在前面 活着的软件,不会停止更新.只是每一个软件更新(对外开放新版本)的频率不同.曾经TBtools平均每天一更,当然现在已经很长时间不更新了,或者确切地说,并不那么频繁地推出新版本了.这或许有两个 ...

  • 改造-基因组浏览器-IGV(九)- 快速[t]BlastN结果

    写在前面 高通量测序技术普及,大量物种基因组被测序.组装质量逐步提高,但基因结构注释的质量则一如既往的让人担心.错误的基因结构注释,可能会使得你得到错误的结果.比如: 不准确的基因结构注释,使得编码蛋 ...

  • IGV-sRNA,改造出一个适合小RNA分析的基因组浏览器

    写在前面 IGV是我认知范围内,得到最广泛使用的基因组浏览器.作为一个强大的软件,其支持 各种测序数据的可视化,但对于小RNA测序数据的支持,却仍然一般.于是,我对其进行了几次小的修改,完成一个相对适 ...