How can I count one specific exon in RNA-seq data
I would like to quantify the number of exones of one particular exone in a specific gene. I have rnaseq data.
• 4,304 views
•
link
updated
by
Charles Warden
•
written
by
carl.h •
•
written
by
carl.h •
2 answers
you can use DEXSeq on R see this.
• 0 views
•
link
written
by
firatuyulur
DEXSeq is a good suggestion.
If you wanted to use JunctionSeq for splicing analysis, QoRTs can provide exon and splice junction counts.
• 0 views
•
link
written
by
Charles Warden
Log in to answer this question.
More posts like this
-
quantifying splicing transcripts for nanopore RNAseq
written by jshen •
-
Transcript quantification
written by firefox91 •
-
Accuracy of Methods for Gene Variant Detection and Quantification: Kallisto Pseudoalignment vs. STA…
written by sohail.mbio •
-
Gene enrichment or rescue dataset
written by Dirk
-
How to quantify mutations by type (e.g. A>T, T>A, G>C, ...) within reads mapping to genes
written by gdeniz •
-
Suggestions on open databases
written by carl.h •
-
How does regularized log transformed data from DESeq2 relate to read counts?
written by amymccarthy91 •
-
Specific sequence abundance using RNA-Seq
written by user230613
-
Assign sample identities based on known gene signatures in RNA-seq data
written by igor
-
How to compare correlations between tumor samples
written by carl.h •
I assume you mean counts for a particular exon? Align to the genome/transcriptome as usual. Use a GTF containing just the data you are interested in for counting with
featureCountsat exon level.Yes that is right. Thank you!
For a follow up question, you don't happen to know how to count all exones in one gene? I have the gtf annotations for all the exones in one gene and I would like to get the number of each exon for the gene.
As long as you have the annotations in your GTF file you can count the exons separately (e.g.
-t exon, which is the default for featureCounts. Take a look at the manual here.)Thank you it worked perfectly, I had to add -f to count feature.