featureCounts的MultiQC结果解读
2023-02-26 10:16:33 时间
本文共2456字,阅读约需7分钟。
在前面我们介绍了MultiQC这个软件,接下来将分期解读MultiQC对多个软件处理后的结果。
本期解读转录组上游分析定量软件featureCounts的结果。
在定量完成以后,我们可以通过如下命令拿到经MultiQC处理的结果。
multiqc result.txt.summary
其中result.txt.summary为多个样本定量的summary文件,如下
Status ../genome/W0d_H5S_1.sam ../genome/W0d_H5S_2.sam ../genome/W0d_H5S_3.sam ../genome/W12h_H5S_XM_1.sam ../genome/W12h_H5S_XM_2.sam ../genome/W12h_H5S_XM_3.sam ../genome/W12h_H5S_ZJ_1.sam ../genome/W12h_H5S_ZJ_2.sam ../genome/W12h_H5S_ZJ_3.sam ../genome/W4d_H5S_XM_1.sam ../genome/W4d_H5S_XM_2.sam ../genome/W4d_H5S_XM_3.sam ../genome/W4d_H5S_ZJ_1.sam ../genome/W4d_H5S_ZJ_2.sam ../genome/W4d_H5S_ZJ_3.sam
Assigned 19014653 16380425 16001137 16397324 17022447 17257563 20357095 17418628 16731705 18129190 17029418 18299244 17599984 20027615 19097183
Unassigned_Unmapped 3104042 3092936 3046480 2780825 2849854 2843823 3996902 2829139 2725304 2796492 2473837 2715289 2784885 3880517 2871043
Unassigned_Read_Type 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_Singleton 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_MappingQuality 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_Chimera 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_FragmentLength 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_Duplicate 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_MultiMapping 7002729 6685469 5779516 6012470 6271135 6299551 7183720 6072958 5792098 7572696 6686149 6845406 6580915 7348777 6398929
Unassigned_Secondary 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_NonSplit 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_NoFeatures 4357138 3544196 3595771 3579708 3704939 3686892 4470750 3849596 3708928 4190702 3926947 4122400 4093841 4545252 4264114
Unassigned_Overlapping_Length 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Unassigned_Ambiguity 19762222 17364090 16896847 17453455 17450339 17164163 21053322 17376653 16844437 19114973 16903418 18388310 18504273 20907632 18764383
处理后得到文件multiqc_report.html,在浏览器中打开即可。
featureCounts的结果一共有两个部分,General Stats和featureCounts。
General Stats
在此部分,我们能直观的看到每个样本有多少reads单一比对到基因注释上,所占的百分比为多少。
举个栗子,W0d_H5S_1有19.0 millions条reads单一比对到基因注释上,占所有reads的35.7%。
General Stats
通过小标题下方的小工具,我们可以对结果作图,看看样本的整体情况。
featureCounts
在此部分,assignments的情况被进一步图表化的展示出来,并且可以选择Number of Reads和Percentages两种展示方式。
Number of Reads
Percentages
举个例子,样品W0d_H5S_1中,有19 014 653条(比General Statistics更精确)reads单一比对到基因注释上,比例为35.7%,其他几项分别为:
- Unassigned:Unmapped:没有比对到基因组上的reads(可能为测序过程中污染的reads)。
- Unassigned:Multi Mapping:覆盖了多个feature的reads,官方解释为【A multi-overlap read or fragment is one that overlaps more than one feature,or more than one meta-feature when summarizing at the meta-feature level】。
- Unassigned:No Features:比对到基因组上尚未注释部分的reads。
- Unassigned:Ambiguity:无法明确定义位置(模糊不清)的reads。
参考资料:
- featureCounts: an efficient general purpose program forassigning sequence reads to genomic features.(2013)
相关文章
- Jgit的使用笔记
- 利用Github Action实现Tornadofx/JavaFx打包
- 叹息!GitHub Trending 即将成为历史!
- 微软软了?开源社区讨论炸锅,GitHub CEO 亲自来答
- GitHub Trending 列表频现重复项,前后端都没去重?
- Photoshop Elements 2021版本软件安装教程(mac+windows全版本都有)
- (ps全版本)Photoshop 2020的安装与破解教程(mac+windows全版本都有)
- (ps全版本)Photoshop cc2018的安装与破解教程(mac+windows全版本,包括2023
- 环境搭建:Oracle GoldenGate 大数据迁移到 Redshift/Flat file/Flume/Kafka测试流程
- 每个开发人员都要掌握的:最小 Linux 基础课
- 来撸羊毛了!Windows 环境下 Hexo 博客搭建,并部署到 GitHub Pages
- 超实用!手把手入门 MongoDB:这些坑点请一定远离
- 【GitHub日报】22-10-09 zustand、neovim、webtorrent、express 等4款App今日上新
- 【GitHub日报】22-10-10 brew、minio、vite、seaweedfs、dbeaver 等8款App今日上新
- 【GitHub日报】22-10-11 cobra、grafana、vue、ToolJet、redwood 等13款App今日上新
- Photoshop 2018 下载及安装教程(mac+windows全版本都有,包括最新的2023)
- Photoshop 2017 下载及安装教程(mac+windows全版本都有,包括最新的2023)
- Photoshop 2020 下载及安装教程(mac+windows全版本都有,包括最新的2023)
- Photoshop 2023 资源免费下载(mac+windows全版本都有,包括最新的2023)
- 最新版本Photoshop CC2018软件安装教程(mac+windows全版本都有,包括2023