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

登录

2: 关注

1151: 浏览

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

相关问题

安捷伦（Agilent）的数据处理软件Feature Extraction Software的默认的数据标准化的算法是什么？最后得到的值是否做了log2处理？是否是RMA或者MASS的一种？: 455 浏览 1 关注 0 回答 0 评论

【坐标北京】近期准备做一些RNA-seq，想请问一下大家哪些公司比较靠谱？-w-: 754 浏览 2 关注 1 回答 0 评论

请问芯片数据分析的一般流程及涉及的常用算法？: 1048 浏览 2 关注 3 回答 0 评论

请问一下，我的噬菌体基因组fasta文件还是config打头的，是不是需要进一步拼接成scaffold？还是挑选最大的config进行后续分析？: 552 浏览 2 关注 1 回答 0 评论

请问有谁知道这张图是如何绘制的吗，在NCBI的Conserved Domians模块没找到可以直接导出的地方: 726 浏览 2 关注 1 回答 0 评论

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

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

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

请问什么是基因组共线性分析？有什么意义？: 480 浏览 1 关注 0 回答 1 评论

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

推荐内容

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

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

问如何利用利用TPM或者FPKM完成DESeq2完成的工作?: 1154 浏览 2 关注 2 回答 0 评论

问麻烦问下，转录组差异分析，我需要筛掉比对率低的bam文件或者外显子比对率的样本吗: 899 浏览 2 关注 1 回答 0 评论

问如何根据转录组数据得到新转录本？如何验证一个基因的多个转录本？: 1020 浏览 2 关注 2 回答 0 评论

问无对照的转录组数据如何寻找差异表达基因: 1006 浏览 2 关注 2 回答 0 评论

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

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

问转录组定量结果为0: 1745 浏览 2 关注 3 回答 0 评论

问用log2(fold_change)数据做热图，遇到想要表达的目的基因log2(fold_change)值为inf或-inf时，怎么办？: 574 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024