Entering edit mode
2.5 years ago
tien
▴
40
Hello, I'm trying to map geneID
(from output of count matrix from cellranger
with latest reference version of human genome
) to seqname
from EnsDb.Hsapiens.v86
database. However, there exists 3k out of 33k geneID that cannot be mapped. Is there any other updated version of EnsDb.Hsapiens database that matches reference file used in cell ranger? Or can I use reference file from cell ranger with mapIds
instead of EnsDb.Hsapiens database?
gn_location <- mapIds(EnsDb.Hsapiens.v86
, keys = rownames(seurat_object_from_cellranger)
, column = "SEQNAME", keytype = "GENEID")
Thanks for your help.