Question: Transcript features and annotations
0
gravatar for Sergio Martínez Cuesta
18 months ago by
Cambridge, UK

I am planning a trancriptome alignment of iCLIP sequencing data. How can I link UCSC hg19 transcript ids with transcript features e.g. coding, non-coding, lncRNA, antisense, pseudogene ...?

I downloaded iGenome's UCSC hg19 reference genome and used the genes.gtf file available with the download to prepare reference sequences for RSEM:

rsem-prepare-reference --gtf genes.gtf --bowtie2 genome.fa ../RSEM_bowtie2/genome

The above command genererates a list of files, one of them being genome.transcripts.fa, a fasta file containing 51398 sequences, one for each transcript (NM_130786, NR_015380, NM_001198818 ...) as defined in the genes.gtf file.

Once I perform a transcriptome alignment using rsem-calculate-expression, how can I then link each transcript id with transcript features such as the ones mentioned above?

Any ideas would be helpful.

transcript alignment • 546 views
ADD COMMENTlink modified 18 months ago by Biostar ♦♦ 20 • written 18 months ago by Sergio Martínez Cuesta60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 730 users visited in the last hour