how to know the isoform length and exon number of particular gene from rna sample
1
0
Entering edit mode
19 months ago

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 • 452 views
0
Entering edit mode

0
Entering edit mode

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

1
Entering edit mode
19 months ago

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

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

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

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?