2
关注
3493
浏览

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

查看全部 2 个回答

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

hisat2 + stringtie + deseq2 标准流程

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

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

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

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

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

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

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

关于作者

问题动态

发布时间
2023-07-31 10:34
更新时间
2023-08-15 17:13
关注人数
2 人关注

相关问题

16s rRNA分析中的标准化
请问哪位老铁有ADMIXTURE的软件包吗?原网址UCLA的软件下载界面访问不了
请问一下,如何根据基因名批量下载CDS区的序列?
植物TWAS流程
由cutesv流程鉴定到的SV有特别多的缺失基因型
请问如何使用R语言绘制散点图
尿代谢组正负离子数据标准化是否可以均用正离子检测出来的肌酐峰
请问hemmer是否可以使用转录组搜索HMM库?
请问在不考虑测序质量的情况下,如何根据bam文件获得比对到指定position的碱基?
DiffBind 标准化数据

推荐内容

麻烦问下,转录组差异分析,我需要筛掉比对率低的bam文件或者外显子比对率的样本吗
samtools 有错误
二代转录组去除批次效应
基因表达
植物TWAS流程
转录组
HISAT2建立索引
转录组和代谢组联合分析的方法有那些?
普通转录组的跨物种分析
肿瘤亚型分析
All Rights Reserved Powered BY WeCenter V4.1.0 © 2025