Hi,

I was trying to get the gene details of refseq ids (these are output from blastn) with the following code:

=================

  library(biomaRt)
  mart <- useMart("ensembl", dataset="hsapiens_gene_ensembl")
  attr <- c("refseq_mrna",'hgnc_symbol', 'chromosome_name','start_position','end_position')
  ids <- c("NM_001270493", "NM_001185076","NM_001185075","NM_001185082")
  
  xx <- getBM(attributes=attr,
              filters = "refseq_mrna", values = ids, mart = mart)
  
  print(xx)
===============

The results are:
>   print(xx)
   refseq_mrna hgnc_symbol chromosome_name start_position end_position
1 NM_001185075        FMR1               X      146993469    147032645
2 NM_001185076        FMR1               X      146993469    147032645


How can I get details of the other two refseq ids that do not show up in the results? 


thanks!

Tim

	[[alternative HTML version deleted]]

