目前主流差异表达分析软件(Deseq2)均是基于counts值,而Hisat2结合cufflinks只能得到FPKM值,那得想办法利用Hisat2软件mapping后得到的bam文件,获得counts值。有两种主流软件:一种是htseq软件(htseq2),但是运行时间长,这里使用一个性能相对较好软件:featurecounts(subread软件-featurecounts)。我使用了自己甘蔗转录组数据bam文件,测试结果表明,使用该软件,一个3.3G的bam文件,所需时间约2-3min就能运行获得。
2020-12-26-生信入坑笔记(三)-转录组分析-Hisat2与featurecounts结合获得转录组counts值
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
推荐阅读更多精彩内容
- 使用工具 fastp(质控), hisat2(比对), samtools(sam文件转bam文件), featur...
- 要求 实现这个功能的软件也很多,还是烦请大家先自己搜索几个教程,入门请统一用htseq-count,对每个样本都会...
- 5月以来,哪怕对市场风向再不敏感的人,也感觉到阵阵凉意。二级市场连续下挫,一级市场融资环境恶化,不论企业融资数量还...