Entering edit mode
                    3.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.