Cannot map some GENEID to SEQNAME from EnsDb.Hsapiens.v86 database
0
0
Entering edit mode
2.0 years ago
tien ▴ 10

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.

annotation cellranger • 363 views
ADD COMMENT

Login before adding your answer.

Traffic: 2056 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6