Entering edit mode
4.8 years ago
Rituriya
▴
40
Hi all,
I have been using the script provided by StringTie named 'prepDE.py' for getting gene and transcript abundances. It was working fine with another set of mouse samples but failing now with these set of human samples.
/Softwares/prepDE.py -i grouping-info.txt
0 Sample1
1 Sample2
Traceback (most recent call last):
File "/Softwares/prepDE.py", line 255, in <module>
geneDict.setdefault(geneIDs[i],{}) #gene_id
KeyError: 'Sample2.101064.1'
It seems there is some issue with the Gene ID - 101064.1. I checked for this name in the gtf files, but only 101064 is found. I tried removing a few samples as well but still it fails after the first sample run.
Please help if anyone has seen such an issue. I am using python version 3.5.
Thank you, Rituriya.