【直播】我的基因组 34:Qualimap对比对好的bam文件进行统计及可视化

上一讲我们说道Qualimap这个软件,那我们现在仔细了解下这个软件吧!

比如下面这个覆盖度与测序深度的图来说:

如果我们自己要对55G的bam文件进行统计,可以sort好之后用samtools的mpileup功能来自己写脚本来区分染色体,分别统计覆盖度和测序深度,再绘图。但其实有很多优秀的工具也可以做到,而且是一步到位,还附送非常多的其它统计指标,比如Qualimap,有点类似于fastqc,也非常好用!

软件安装如下:

## Download and install Qualimap ## http://qualimap.bioinfo.cipf.es/ cd ~/biosoft mkdir Qualimap && cd Qualimap wget https://bitbucket.org/kokonech/qualimap/downloads/qualimap_v2.2.1.zip ## readme http://qualimap.bioinfo.cipf.es/doc_html/index.html ## example results :http://kokonech.github.io/qualimap/HG00096.chr20_bamqc/qualimapReport.html ~/biosoft/bamtools/bamtools/bin/bamtools ~/biosoft/Qualimap/qualimap_v2.2.1/qualimap --help

软件安装之后使用非常简单,而且我们只使用它的bamqc的功能命令如下:

~/biosoft/Qualimap/qualimap_v2.2.1/qualimap --java-mem-size=35G bamqc -bam ~/data/project/myGenome/bamFiles/P_jmzeng.final.bam

而且它不需要再对每条染色体分别处理了,因为它本身就自带这个统计机制。至于这个软件做了什么呢,我已经把它的报告上传到网络了,大家可以点击查看原文来了解或者复制链接到浏览器 http://www.biotrainee.com/jmzeng/blogMyGenome/P_jmzeng_Qualimap_results/qualimapReport.html,其实也就是我们前面对bam文件的一些解析而已,没什么意思,讲这两个软件,只是告诉大家,对于生信工作者来说,学习一个新的软件的用法应该是非常简单的事情。

文:Jimmy、吃瓜群众

图文编辑:吃瓜群众

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

原文发表时间:2016-12-31

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

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏FreeBuf

2014年漏洞最多的操作系统:苹果Mac OS X、iOS和Linux

微信号:freebuf 无论是Mac OS X还是iPhone iOS,苹果系统一直被认为是最安全的操作系统。然而在2014年,苹果系统却获得了漏洞最多的“殊荣...

1866

使用Bluemix,NoSQL DB和Watson创建云应用程序

大家好,自从我上次写文章已经很久了。事实上,这几年,我总是忙于工作。我现在是IBM的Bluemix平台的云架构师。我在Tomcat服务器上用Web应用程序编写了...

1836
来自专栏IT笔记

三分钟读懂TT猫分布式、微服务和集群之路

针对新手入门的普及,有过大型网站技术架构牛人路过,别耽误浪费了时间,阅读之前,请确保有一定的网络基础,熟练使用Linux,浏览大概需要3-5分钟的时间,结尾有彩...

28910
来自专栏北京马哥教育

上班第一天,一个合格的运维应该做什么?

运维行业正在变革,推荐阅读:30万年薪Linux运维工程师成长魔法 作为一名运维工程师,如果你在春节放假期间没有被报警电话和邮件吵醒过,那说明你在放假前的准备...

3598
来自专栏编程一生

架构师之路--从业务角度谈缓存的选型

1595
来自专栏FreeBuf

什么是“移动端应用协作”(MAC)攻击?

Intel安全团队最近表示,他们检测到恶意代码被发布在数以千计的安卓包里。黑客对这些代码进行组合后,对没有防备的手机用户发起攻击。 这种攻击被称为“移动端应用协...

1719
来自专栏FreeBuf

窃听电话的Hacking Team RCSAndroid木马

安卓设备小心:4.0-4.3版本都可以被RCSAndroid 搞定。 安卓平台上的远程控制木马RCSAndroid是目前曝光的安卓中最专业、最复杂的恶意程序之一...

2148
来自专栏安恒信息

携程曝重大安全漏洞 客户信用卡信息或遭泄露

3月22日,乌云平台连续披露了两个携程网安全漏洞,漏洞发现者称由于携程开启了用户支付服务借口的调试功能,导致携程安全支付日志可被任意还可读取,日志...

2668
来自专栏黑白安全

漏洞评估指南 {Vulnerability Assessment}

漏洞评估的四步指南 这是一个建议的四步法,使用任何自动或手动工具启动有效的漏洞评估流程。

1212
来自专栏FreeBuf

每16台Android手机中,就有一台受BadKernel漏洞的影响

根据国外媒体的最新报道,安全研究专家在Google的V8 JavaScript引擎中发现了一个安全漏洞(BadKernel),该漏洞将会间接影响到Android...

1885

扫码关注云+社区