Error: could not locate transcript ENST00000624361解决stringtie定量报错

生信喵

背景
       近期做的RNA-Seq数据上游，hisat2+stringtie流程，使用的stringtie官方的prepDE.py脚本定量。

1. prepDE.py -i sample_lst.txt

复制代码

      命令报错：

1. Error: could not locate transcript ENST00000624361 entry for sample Sample10
   
2. Traceback (most recent call last):
   
3.   File "/share/home/xiehs/tools/stringtie-2.1.4.Linux_x86_64/prepDE.py", line 281, in <module>
   
4.     geneDict.setdefault(geneIDs<i>,{}) #gene_id
   
5. KeyError: 'ENST00000624361'</i>

复制代码

解决方案
       经过一系列的google探索，得出了这是因为前面使用的stringtie版本不对，下面贴出官方的回答：

1. Thank you so much @Gin-Wang I took a look at the data you uploaded and I can confirm now that the problem is real and correctly reported by users in this thread (and in #234), stringtie 2.0 with the -e option outputs a few STRG entries which should NOT be there, as the output should only contain the MSTRG entries that were given in the -G file.
   
2. We will fix this -e problem soon and report back here, until then please use stringtie v1.3.6 for quantification via prepDE (for short reads there shouldn't be a significant difference).

复制代码

      意思就是选用1.3.6的版本，他们会尽快修复。
       那么如何快速的切换版本呢？自然是使用conda管理。

1. conda activate test
   
2. conda install stringtie=1.3.6-0

复制代码

      使用我们就用绝对路径写命令即可。

1. /share/Software/miniconda3/envs/test/bin/stringtie

复制代码

      跑完gtf文件，用脚本定量，果然没有报错了。
       特此分享给大家。