How to change stringtie gene_id names with ref_gene_id when performing featureCount
1
0
Entering edit mode
@iraiamaialen-19272
Last seen 6.0 years ago

Hi everybody!

I would like to know if there is a way to change stringtie "geneid" with "refgene_id" when performing featureCount to construct the count table for edgeR.

The thing is that I have use stringtie to create a merged gtf file, using all the files I want to compare in the DEG analysis. After that I have use featureCount to create the count matrix, and as featureCount uses gtf file "geneid" by default, in the result count table, I have a very big list of MSTRG------ named genes. I would like to know if using an option in featureCount gives the posibility to change from geneid to refgeneid that is the ENSMBLE ID.

Maybe -g < string > (GTF.attrType), option and ther use "refgeneid"? Would this change be correct? Or maybe I have realized now that there is another option named −−extraAttributes < string > (GTF.attrType.extra) and I can add there refgeneid?

Thanks in advance,

Iraia

featureCount rna-seq gene_id re_gene_id stringtie • 2.2k views
ADD COMMENT
2
Entering edit mode
Wei Shi ★ 3.6k
@wei-shi-2183
Last seen 7 days ago
Australia/Melbourne

If "refgeneid" is included as an attribute type in your GTF file, you can set GTF.attrType="refgeneid" to instruct featureCounts to count reads to this type of genes.

The parameter "GTF.attrType.extra" does not work for this purpose. It just extracts extra annotation data from your gtf file and adds them to featureCounts output.

ADD COMMENT

Login before adding your answer.

Traffic: 1066 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6