I am trying to use the getGeneLengthAndGCContent() function from the EDASeq R package to obtain gene lengths and GC content for each gene in a honeybee (apis mellifera) dataset. An example subset of these genes are as follows:
"GB40001" "GB40002" "GB40003" "GB40004" "GB40005" "GB40006"
Upon running the command on an example gene as follows and using the code for apis mellifera ("ame"):
I receive the error:
Error in getGeneLengthAndGCContent("GB40001", org = "ame", mode = "org.db") : Only ENTREZ or ENSEMBL gene IDs are supported.
How can I convert my list of genes (in format of GB40001) into a list of ENTREZ or ENSEMBL gene IDs?
Thank you for any support on accomplishing this step.