Google searching for NM_002084 gives the following result:
This, I assume, is the latest version v4, hence the .4
suffix.
Searching for previous versions I get the following results, along with notes saying it was updated or removed.
This sequence has been updated. See current version.
Record removed. This record was replaced or removed.
Using biomaRt I can get the latest(?) version as follows:
library("biomaRt")
# define db
ensembl <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
# get refseqs
getBM(attributes = c('refseq_mrna',
'chromosome_name',
'transcript_start',
'transcript_end',
'strand'),
filters = c('refseq_mrna'),
values = list(refseq_mrna = "NM_002084"),
mart = ensembl)
# refseq_mrna chromosome_name transcript_start transcript_end strand
#1 NM_002084 5 151020438 151028992 1
But querying for specific versions gives nothing:
getBM(attributes = c('refseq_mrna',
'chromosome_name',
'transcript_start',
'transcript_end',
'strand'),
filters = c('refseq_mrna'),
values = list(refseq_mrna = c("NM_002084.1", "NM_002084.2", "NM_002084.3", "NM_002084.4")),
mart = ensembl)
# [1] refseq_mrna chromosome_name transcript_start transcript_end strand
# <0 rows> (or 0-length row.names)
Question: How can I get all versions (preferably using R)?