Entering edit mode
6.2 years ago
k.kathirvel93
▴
310
Hi EveryOne,
I am trying STRINGtie for RNA-seq mRNA quantification after mapping with hisat2. STRINGtie can't produce gene raw counts, so i have used prepDE.py script to obtain gene counts. But, here i am getting only ensembl ID. I want gene name instead of ID. Any solution ? Thanks in advance.
Can you add the command you ran along with a few lines from the input GTF files ?
I am facing the same issue. Please help me in solving this issue.
Hi arshad1292,
I found the way by myself, First go to normal stringtie gene abundance file and reorder the ensemblID column in alphabet order with extent selection and copy gene ID column. Second go to PrepDE.py output file gene_Count_matrix.csv, and reorder the ensemblID column in alphabet order with extent selection and replace the geneID column by the copied geneID column from previous. This you can do in R also very easily. all the best
Could you also, share how do we do that on R? I think I have gone way ahead of that PrepDE.py step!! Thank You!!
Its faster if you do manually.