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

写在前面

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

展示二级结构的同时展示位点测序深度

一开始,这是一个相对复杂的事情,但事实上也是简单的事情。
复杂在于:序列是基因组的,而参考的基因组,只有一个;而测序深度是针对每个数据集不同的。
这个就意味着,同时展示结构域信息和测序深度信息,是一个直接不可调和的事情。
解决办法,前述不可调和,主要原因在于我的思维被RegionOfInterest的RNAfold实现限制。事实上,我们只需要对每个数据集合开放菜单接口即可。

使用

  1. 调整到进行Fold的区间(独立于RegionOfInterest)

  2. 在ReadLenTrack上调出菜单
    鼠标右击ReadLenTrack,即可看到菜单

  1. 点击则会弹出Fold窗口
    分别点击,可看到推文最初的图片

图片说明

从图片中科看到,

  1. 左侧是IGV窗口,按照正常IGV浏览器操作即可

  2. 右侧上端是当前链折叠的结果;下端是反向互补链折叠的结果

  3. 鼠标在IGV窗口移动的同时,竖线会指向某个碱基位置,同时,对应的折叠窗口会高亮对应的碱基。(图中,碱基A,对应正向折叠窗口红圈高亮碱基A,对应反向折叠窗口红圈高亮碱基U)

  4. 凡是有reads覆盖的碱基位置,外圈为橙色,否则为黑色

  5. 覆盖率高低在折叠窗口上以红色的深浅程度表示

写在后面

改造和增强是没有终点的,只是这个改造或者增强是否确实存在意义。
前述的两个改造推文均有朋友问到需要这个版本,
那么请联系园艺植物小分子RNA与基因组研究-夏瑞课题组

课题组主页:http://xialab.scau.edu.cn/

(0)

相关推荐

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

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

  • 直播我的基因组(第一阶段)完整目录

    最近的全国巡讲不少人问到我两年前的直播基因组系列教程的完整目录,这里先放出直播我的基因组(第一阶段)完整目录.(悄悄告诉你,后台回复直播可以拿到精排版EXCEL表格!)(然后,点击阅读原文也可以拿到可 ...

  • IGV基因组注释文件说明

    IGV基因组注释文件说明

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

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

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

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

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

    写在前面 还有两个小时....似乎还是有点紧张.想想,那就写写微信推文. IGV改造系列,一直会继续,直到或许是我们觉得似乎没有改造的必要. 在前述的改造记录中(https://www.jianshu ...

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

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

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

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

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

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

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

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

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

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

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

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