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

登录

2: 关注

3492: 浏览

请问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 人关注

相关问题

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

16s rRNA分析中的标准化: 1702 浏览 1 关注 0 回答 0 评论

请问哪位老铁有ADMIXTURE的软件包吗？原网址UCLA的软件下载界面访问不了: 2892 浏览 3 关注 2 回答 0 评论

请问一下大家生物信息学的生物学具体应用: 2067 浏览 2 关注 1 回答 0 评论

请问下这种格式的R语言内容如何选择最小值: 2291 浏览 2 关注 1 回答 0 评论

blastp后，我想提取score的序列，请问使用python应该怎么写这个脚本: 2059 浏览 2 关注 1 回答 0 评论

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

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

由cutesv流程鉴定到的SV有特别多的缺失基因型: 2063 浏览 2 关注 1 回答 0 评论

DiffBind 标准化数据: 2159 浏览 2 关注 1 回答 0 评论

推荐内容

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

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

问转录组和代谢组联合分析的方法有那些？: 2936 浏览 3 关注 2 回答 0 评论

问 HISAT2建立索引: 3476 浏览 1 关注 2 回答 0 评论

问绵羊转录组测序后差异表达基因太多了: 3295 浏览 2 关注 2 回答 0 评论

问我测了一批寄生蜂的转录组样品，FPKM值感觉好奇怪，合理吗？: 3037 浏览 2 关注 3 回答 0 评论

问菌侵染植物的转录组分析中，若菌的一个基因在当前时间没有引起植物基因的变化，过了三个小时后才引起植物基因的变化，如何分析这种现象或如何将这一类基因关联起来，有没有类似的文献推荐？: 2021 浏览 2 关注 1 回答 0 评论

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

问请教一个问题，一个转录因子TF只会调控位于同一条链的基因，还是有可能调控反义链的基因？: 2595 浏览 2 关注 2 回答 0 评论

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

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025