Hello,
I'm trying to simply convert an Ensembl ID to a uniprot gene symbol. I've been trying to find something for this but it's like finding a needle in a haystack.
I think I should use $gene->external_name() from the BioMart API, but I'm not sure how to correctly incorporate this in my script (which modules to use, what elements are needed etc).
I need this because I have hundreds of files with RNAseq data and for each file I remove the unexpressed genes. In the resulting list of expressed genes, there are a lot of pseudogenes, which I do not need. When I convert the Ensembl ID to gene name, the actual genes (with a gene symbol) can be easily extracted. I can do this conversion with biomart or some other tool but that's too much work for hunderds of files. That's why I need a simple function in my script to convert the ID in each line to gene symbol, and insert it in the next column.
Does anyone have an idea how to do this?
Thanks in advance.