[Bioc-devel] UniProt.ws: Error in .mergeList when mapping to UNIPROTKB

[Thomas Dybdal Pedersen]

Hi

When trying to map from a vector of refseq accessions to uniprot id's I get an error:

refseq <- c("YP_139402", "YP_141320", "YP_820357", "YP_006002448", "YP_006040776", 
"YP_006340074", "YP_006340075", "YP_138830", "YP_819836", "YP_006001858"
)

select(UniProt.ws, keys=refseq, keytype="REFSEQ_PROTEIN", cols=c('UNIPROTKB'))

I have traced the error back to .mergeList. The problem seems to be that when mapping to UNIPROTKB a data frame with two columns named ACC+ID is produced making the by='ACC+ID' argument in merge() ambiguous.

The easy fix seems to be to only use 'ACC+ID' for internally keeping track of the mapping and use another name for the UNIPROTKB id?


best

Thomas Pedersen, PhD student at Technical University of Denmark