写在前面
前述,我们已经对IGV进行了超过三次改造,同时我也写了三个推送。写IGV改造系列推送的主要原因,事实是作为课题组成员使用改造的IGV时的参考教程。
在上一次RNAfold的窗口上,我们增加了二级结构上每个碱基覆盖深度的深浅着色。依据惯例,我会发一条朋友圈,同时附上评论信息:似乎已经没有特别值得改造的地方了。天真如我,我收到了这么一个评论:
增加实时计算Track的繁琐之处
在IGV中,伴随窗口滑动,实时更新的Track主要有三个:
- Alignment,即reads的比对位置(这个是默认的Track)
- Coverage,覆盖率,只要是输入测序数据bam/sam文件,那么可以调出这个Track
- Junction,read的断开位置的Track,可以很好的判断可变剪接或SV
实现这类Track的繁琐有二:1. 数据的获取和整理;2. Track的初始化和菜单的增加和调整
在前述,我们已经增加过一个ReadLenTrack,实现的过程是,获取当前区域所有的alignment信息,重新整理计算,随后更新。
而菜单的增加和调整,需要参考CoverageTrack....
在ReadLenTrack经验的基础上,我们可以简单一些地基于上次修改的位置,添加PhasScoreTrack.
PhasScoreTrack的开启方式
时间有限,我暂时并没有像ReadLenTrack一样,增加不同Track之间的互相开启关闭摁钮。在这次修改中,我仅仅在Alignment Track的右键菜单上增加一个开关。如下
从图中可见,绿色箭头对应的是前述增加的ReadLenTrack。从ReadLenTrack来看,这是一个非常明显的21nt的Phas位点;红色箭头对应的是本次增加的PhasScoreTrack。开启这个Track时,我们可以看到21nt的Phasing Score显示。伴随窗口的移动,会进行实时计算并展示。
相关使用方式
植物的phasiRNAs主要有21nt 和 24 nt的两类,但我们并不能保证,不存在22nt 更或者是 23 nt的... 于是在PhasScoreTrack的右键菜单,我们增加一个调整功能,用户可调整PhasLen,以从不同的phas长度观测数据
写在后面
......我想,IGV改造应该暂时告一段落了。事实上,有一个增加了的可以极大地加速sRNA数据可视化的隐藏功能,我暂时没有推出。留作后续调整。
自家广告
在家里呆了几天,我发现我增加了两个业务:
- 卖毛衣(毛衣生产 三十年 提供定制 <比如来样加工,你懂的,A货>)
- 卖手工橄榄菜(其实还卖酸菜),目前是¥30一罐(~500g),老少咸宜,优质美味下饭菜,微信CJ-Chen-0410。(这个算祖传的手艺,手工的,物有所值)
修改版IGV的获取方式
近日有多个朋友联系过来,想要使用这个改造后的IGV。嗯...
我个人的想法是:
- 付费,如资助XiaLab课题组出游一次,大体价格是3K,那么将获得本年度(如果还有更细的话)的我经手的IGV更新功能。
- 直接联系课题组PI即RX获取,课题组主页为 http://xialab.scau.edu.cn/