【直播】我的基因组 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

全球最大同性社交软件Grindr存在漏洞,泄露用户信息及位置

美国NBC的一份报道称,一款名为Grindr的交友应用程序存在两个安全问题,它可以暴露超过300万用户的信息,包括那些选择不共享这些信息的人的位置数据。此次漏洞...

1442
来自专栏FreeBuf

别让你的手机成为恶意软件的温床

手机早已成为人们生活中不可或缺的一部分,在娱乐,移动支付,社交等方面都起到举足轻重的作用,可以说在当今社会手机就是人们重要的“生存”工具之一。但是在享受时代带给...

1350
来自专栏小樱的经验随笔

网站渗透测试原理及详细过程

网站渗透测试原理及详细过程 零、前言 渗透测试在未得到被测试方授权之前依据某些地区法律规定是违法行为。 这里我们提供的所有渗透测试方法均为(假设为)合法的评估...

5359
来自专栏腾讯云安全的专栏

黑客教主 TK 说:咱们谁都没错,就是不该在一起

1946
来自专栏FreeBuf

担心投票系统被黑,详解荷兰政府为何决定在大选中使用手工计票

在3月的荷兰大选即将到来之际,安全专家Sijmen Ruwhof通过Youtube上有关荷兰大选使用的计票软件介绍视频,对荷兰大选中使用的投票计票软件OSV p...

2017
来自专栏安恒信息

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

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

2788
来自专栏FreeBuf

窃听电话的Hacking Team RCSAndroid木马

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

2288

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

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

1876
来自专栏漏斗社区

看!我手里有个Email收集神器

众所周知,在工作中,大家用到最多的便是Email来传输工作内容,因此,像一些VPN信息、服务器账号密码、公司人员清单等敏感数据,均会通过Email进行传输,并且...

1593
来自专栏FreeBuf

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

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

1779

扫码关注云+社区