I have aligned raw RNA-seq reads against a reference genome with STAR. I would like to extract the intronic coordinates from obtained bam files. Any suggestion?
This isn't information that one extracts, since it's not held within the BAM files. Instead your steps will be:
The second step may require some custom code, depending on whether you want introns that overlap exons (e.g., due to having multiple transcripts for a gene) or not.
Hi Ryan, The GTF file, that I obtained from the stringTie, is missing strand information. I need strand information to check the splice sites sequence (both at 5' and 3')
Do you have idea, how can I deal with that ?
Did you use a strand-specific library prep?
We didn't use strand-specific library prep.
Then you can't know which strand the transcripts are from.
Today, I got to know that mRNA library is stranded. In that case, how can I obtained the strand information or why could it be missing ?
You need to tell your aligner and stringTie that.
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy