transcript to gene mapping after sleuth run
0
1
Entering edit mode
7.0 years ago
bioinfo8 ▴ 230

Hi,

I have ran kallisto and sleuth for my RNA-seq data and used the following for transcript to gene mapping before running sleuth:

mart_11 <- biomaRt::useMart(biomart = "ENSEMBL_MART_ENSEMBL", dataset = "organism_gene_ensembl", host = 'mar2015.archive.ensembl.org')

t2g <- biomaRt::getBM(attributes = c("ensembl_transcript_id", "ensembl_gene_id", "external_gene_name"), mart = mart_11)

t2g <- dplyr::rename(t2g, target_id = ensembl_transcript_id, ens_gene = ensembl_gene_id, ext_gene = external_gene_name)

Now I have sleuth output and realized that some of the transcripts has NA in ens_gene and ext_gene. I found out that in the recent assembly there is information regarding this. Would anyone can guide me an easy way to update this information in my sleuth output and also adding gene names instead of transcript IDs. Otherwise, I have to run sleuth again by creating new t2g.

Thanks!

sleuth kallisto RNA-Seq R biomaRt • 1.9k views
ADD COMMENT

Login before adding your answer.

Traffic: 3327 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