Question: How can i get gene name instead of gene ID from prepDE.py after STRINGtie?
0
gravatar for k.kathirvel93
10 weeks ago by
k.kathirvel93120
India
k.kathirvel93120 wrote:

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.

rna-seq next-gen gene • 182 views
ADD COMMENTlink written 10 weeks ago by k.kathirvel93120

Can you add the command you ran along with a few lines from the input GTF files ?

ADD REPLYlink written 10 weeks ago by Vijay Lakhujani3.1k

I am facing the same issue. Please help me in solving this issue.

ADD REPLYlink written 26 days ago by arshad129220

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

ADD REPLYlink modified 26 days ago • written 26 days ago by k.kathirvel93120
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1130 users visited in the last hour