Hi all,
I have RNA-Seq data (ctrl vs kd) and performed a bunch of analyses. One of the analyses I performed was via rMATS tool to find possible retained introns. In the rMATS output, I do get coordinates for individual retained introns but I think depending on the rMATS output, the exact positions listed might be including the last/first exonic nucleotide, or excluding it.
I want the true length of each retained intron to perform possible frameshift analyses (intron length divided by 3), thus even one nucleotide offset can make my analyses go wrong. With that said, how can I extract intron sequence so then I can calculate the length of it?
Thank you.
1 answer
I would use STAR to align the raw reads against the dataset, and check this manually in detail.
Get the BAMs, and use IGV to inspect. I think this is the only method if you need that accuracy.
Log in to answer this question.
-
usually retained introns
written by adi.gershon1 •
-
How to calculate the percentage of hydrophobic amino acids?
written by Apex92 •
-
How the "IJC" is calculated in rMATS tool?
written by Apex92 •
-
Differential gene expression analysis with no replicates using edgeR
written by Apex92 •
-
How robust is rMATS when swapping samples?
written by Apex92 •
-
Why do I see 3-nucleotide periodicity in RNA-Seq?
written by Apex92 •
-
ANNOVAR - RefSeq Annotation
written by Maguelonne •
-
How many of the reads end up being counted by HTSeq count
written by Apex92 •
-
Alternative splicing event annotations method
written by Rashedul Islam
-
Searching Introns On Ucsc Genome/Table Browser
written by jyryc •
