报告导航与解读:从全局到细节
报告生成了,一大堆图表、表格、警告标志,怎么下手?
这是我最常被问到的问题。说白了很多同学跑完 MultiQC,打开 report 一看,满屏花花绿绿,直接懵了。
看报告,最忌讳的就是一头扎进细节里。我个人的习惯,是先瞄一眼全局统计表,心里有个谱。然后再按模块往下走,有问题的地方我会优先点开。
这一章,我带你把 MultiQC 报告的骨架摸清楚。怎么导航、怎么看、怎么快速定位问题。
先看全局统计表
打开 MultiQC 报告,第一眼看到的就是一个表格。这个表列着所有样本,每一列是关键的质控指标。
像是 FastQC 的原始读长、比对率、插段大小、Duplication Rate 等等。我在项目中遇到过,数据量大跑一整天,结果报告一打开,发现有个样本的 % Duplicates 飙到了 90% 以上。嗯,马上就标记出来,看都不用看下面那些图了。
个人建议:全局统计表是巡警。从上往下扫一遍,谁数值偏离平均值太多,谁颜色异常,一目了然。不要跳过这步。
有的指标列还有颜色标记,绿色=正常,橙/红=警告或失败。这是 MultiQC 的核心设计——让眼睛替大脑过滤信息。
其实你还可以点击表头排序。找最高值、最低值、离群值,很方便。
# 举个例子,在自定义配置里调整警告阈值
general_stats:
# 调整FastQC per_base_quality的警告阈值
fastqc_per_base_sequence_quality:
warn: 25
fail: 20
这个配置文件我经常用。数据质量基线不一样,我就把警告线改一改,让它更贴合我的项目。
左侧导航:模块化视图
往下滑,或者点左侧的导航栏,你就能看到每一个模块的图表了。
MultiQC 把每一种分析工具的结果放在一个独立区块里。FastQC 一个区块,Cutadapt 一个区块,STAR 一个区块……清清楚楚。
我刚开始用的时候总想一次把全部看完,结果脑子过载。后来学会一个技巧:只看有问题的模块。
技巧:在报告左侧菜单里,点开有橙色或红色标记的模块,优先处理。绿色模块可以先忽略。
每个模块内部,通常包含多个图表。比如 FastQC 区块下面,有 Per Base Sequence Quality、GC Content、Adapter Content 等等。每个图表都配有样本名和图例,鼠标悬浮就能看到具体数值。
我把常见模块的应用场景列了一下:
| 模块 | 主要用途 |
|---|---|
| FastQC | 基础测序质量:碱基质量、GC含量、接头污染 |
| Cutadapt / Trimmomatic | 看切了多少数据、去接头效果如何 |
| STAR / HISAT2 | 比对率,mapping 结果是否正常 |
| Picard / SAMtools | Duplication 情况、覆盖率分布 |
嗯,不同项目的关注点不一样。RNA-seq 你就死磕比对率,WGS 你就盯住覆盖度和 duplication。不要贪多,抓核心。
交互式图表:让数据动起来
MultiQC 的图最爽的地方,是它们都是交互式的。
什么意思?就是你可以缩放、平移、点击显示数值、隐藏某些样本。这在排查异常时特别有用。
举个例子,我在一个 RNA-seq 项目里看到一个样本的 GC 含量曲线歪了,跟其他样本差别很大。如果是静态图,我只能猜。但是交互式图里,我可以点那个样本的名字,一下子把它单拎出来看,其他样本全部隐藏掉。哦,原来是那个样本混入了 rRNA,污染了。
注意:交互式功能依赖 JavaScript。你是离线打开报告,浏览器可能需要启用 JS。如果报告完全空白,先检查你是不是用 file:// 协议打开的,有时需要 --interactive 选项。
图表常见的交互操作:
- 滚动缩放:在图表区域滚动鼠标,可以缩放 y 轴或 x 轴。
- 拖拽平移:按住鼠标左键拖动,可以移动视图区域。
- 双击重置:双击图表区域,回到默认视图。
- 图例点击:点击图例中的样本名,可以隐藏或显示该样本。
还有一个容易被忽略的功能:鼠标悬浮显示数值。你不需要猜图的纵坐标精度,鼠标移上去,会弹出一个小框,精确到小数点后两位。这在写报告或者做图表对比时,很实用。
我曾经在团队内审时,直接把鼠标悬停的数值截屏贴到报告中。省去了手动读图估数的麻烦。效率可高。
看完了这些,你应该对报告的结构心里有底了。其实说到底,MultiQC 的报告就是一个高度集成的数据仪表盘。它的核心价值,就是帮你从海量的质控报告中,快速找到少数真正需要关注的问题。
最后提一句:如果你在报告中遇到看不懂的图或指标,别慌。MultiQC 的每个图表右上角都有个小问号,点开会有帮助说明。这个小功能,我一开始也没注意,后来发现真是省了不少找文档的时间。