This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to know the isoform length and exon number of particular gene from rna sample

I have 4 sample of rna seq data(control and disease). i did alignment with STAR and RSEM to get the count matrix. now i want to know the size of isoforms in both cases of particular gene? how can i get this?

rna-seq
• 2,259 views • link
written by priyadhawan11234 • priyadhawan11234

Please use ADD COMMENT and not the answer field to keep the thread logically organized.

• 0 views • link
written by ATpoint ATpoint

1 answer

You simply look it up in the the annoation you used tu build the index for RSEM. If you used a fasta file you can look up the isoform length (as it is the same length as the number of nucleotides for that entry). If you want the number of exons you will have to download the corresponding GTF file and extract the information from there.

• 0 views • link
written by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup

I used fasta file for RSEM. I got the transcript length from fasta (that is 6400 bp long). how can i get that sequence from bam file ?

• 0 views • link
written by priyadhawan11234 • priyadhawan11234

Sorry but I do not understand what you want to do. Could you elaborate - and also please comment on why the nucleotide sequence in the fasta file is not good enougth?

• 0 views • link
written by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup

i just want to know which isoforms of particular gene are expressed in both sample and what is length of those isoforms. so how can i get that information from bam file generated by STAR?

• 0 views • link
written by priyadhawan11234 • priyadhawan11234

Log in to answer this question.

More posts like this