请问hisat2+stringtie+deseq2的标准流程 - 做专业的生信问答网站

登录

2: 关注

3237: 浏览

请问hisat2+stringtie+deseq2的标准流程

在网上搜流程，不知道什么时候要不要stringtie --merge，另外我做出来的count的ID都是MSTRG，不知道后续怎么分析

好问题 0 评论收藏举报

查看全部 2 个回答

chatGPT机器人机器人用户来自于: IANA
2023-07-31 10:36

hisat2 + stringtie + deseq2 标准流程

标准的流程包括以下几个步骤：

从NCBI或其他数据库下载参考基因组序列和注释文件
使用hisat2对测序数据进行比对，生成SAM/BAM文件
使用stringtie对比对结果进行转录本组装，生成转录本表达量矩阵
使用deseq2进行差异表达分析

对于你的第一个问题，是否需要使用stringtie --merge，要看你的实验目的和数据类型。

如果你有多个样本或多个测序数据集，想要将它们合并成一个转录本注释文件，那么你需要使用stringtie --merge。

如果你只有一个测序数据集，或者你只对单个样本进行分析，那么你可以不使用stringtie --merge。

对于你的第二个问题，关于count的ID都是MSTRG，你可以按照以下步骤进行后续分析：

将得到的转录本表达量矩阵导出为普通的表格文件，比如CSV格式
打开表格文件，查找MSTRG开头的列，这些列是stringtie生成的新的转录本ID
使用转录本注释文件，将新的转录本ID映射到已知的基因或转录本ID
根据映射结果，将表格中的MSTRG ID替换为对应的基因或转录本ID
继续使用deseq2等工具进行差异表达分析

赞同 1 0评论

关于作者

: 陈文豪7709 注册会员
这家伙很懒，还没有设置简介

1: 回答

0: 文章

5: 问题

问题动态

发布时间: 2023-07-31 10:34

更新时间: 2023-08-15 17:13

关注人数: 2 人关注

相关问题

请问在不考虑测序质量的情况下，如何根据bam文件获得比对到指定position的碱基？: 3841 浏览 2 关注 2 回答 0 评论

请问如何使用R语言绘制散点图: 2372 浏览 1 关注 2 回答 0 评论

聚类数据标准化: 2617 浏览 2 关注 2 回答 0 评论

请问cytoscape是怎么构建网络的？: 1750 浏览 3 关注 1 回答 0 评论

请问RNA-seq采用poly A(+)策略建库，处理数据时若不去除rRNA会对后续分析有何影响？: 1873 浏览 2 关注 1 回答 0 评论

请问UCSC上下载的参考基因组中这个文件是什么含义: 1670 浏览 2 关注 1 回答 0 评论

单细胞RNA-seq分析流程: 2419 浏览 2 关注 2 回答 0 评论

请问无参转录组的GO分析如何进行: 3622 浏览 7 关注 6 回答 0 评论

请问hemmer是否可以使用转录组搜索HMM库？: 2018 浏览 2 关注 1 回答 0 评论

耐药细菌的比较基因组分析流程和机制探索: 1725 浏览 2 关注 1 回答 0 评论

推荐内容

问进行转录组数据分析时，进行cuffdiff后的输出文件gene_exp.diff中，一个基因出现了两个不同的表达量数据，应该如何处理？: 1757 浏览 2 关注 1 回答 0 评论

问植物TWAS流程: 2872 浏览 2 关注 1 回答 0 评论

问 samtools 有错误: 2207 浏览 2 关注 1 回答 0 评论

问转录组测序问题: 2696 浏览 2 关注 1 回答 0 评论

问链特异性文库（mRNA/lncRNA/circRNA）如何将RNA类型分开？: 1989 浏览 2 关注 1 回答 0 评论

问肿瘤亚型分析: 1971 浏览 2 关注 1 回答 0 评论

问 log2(fold_change)可以用来做热图吗？遇到想要表达的目的基因log2(fold_change)值为0时，怎么办？: 2650 浏览 2 关注 1 回答 0 评论

问 ciriRNA表达量如何计算？: 2578 浏览 2 关注 1 回答 0 评论

问 synapser无法加载: 2261 浏览 2 关注 1 回答 0 评论

问不同比对软件出的结果能进行比较吗？: 1852 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025