How To Tell How Many Reads Aligned To Intron In Tophat2

Hi All,

I installed tophat2 and cufflink2 and now I want to obtain the reads aligned in the exon regions and reads aligned in the intron regions. Is there a quick way to get these information? Thanks a lot.

Yours, YING

tophat2 rna-seq

• 3,205 views • link

updated 12.9 years ago by Mikael Huss 480 • | written 12.9 years ago by dingying85 1

1 answer

There are many ways, for example Picard CollectRNASeqMetrics or RSeQC. Or you could use BEDTools on a GTF annotation file, possibly after creating an intron annotation by going through the GTF file and choosing genomic intervals between the exons.

• 265 views • link

written 12.9 years ago by Mikael Huss 480

Thank you so much for the answer!

• 0 views • link

written 12.9 years ago by dingying85 1

More posts like this

How to obtain full alignment results?
written 2.8 years ago by nrqstudent 0

Hi, When I prepare to `align my short reads` to the `human pan-genome` graph genome, the result is only a path, 1. what I want …
How to extract intergenic regions
written 5.1 years ago by Ak 6

I'm trying to obtain intergenic regions from the fasta file containing the genome sequence. I previously identified the exon-intron boundaries of a gene family (using …
Genomic position to genomic feature
written 5.7 years ago by irfanwustl 9

Hi, We have some Methylation samples and obtained some differentially methylated regions (DMR) from those samples. Now we are asking are those DMR makes biological …
Checking whether the reads aligned to specific gene have the specific exon in BAM file
written 6.2 years ago by dh.han 0

Hi guys. So I did a cellranger count on my sample with my own new gene in the gtf annotation file. And from the bam …
genomic reads count
written 7.0 years ago by 13263617869xf 0

Hi everyone,i want to count exon、intron、intergenic、exon-exon、exon-intron reads,now i konw how to use the bedtools like tools to solve tie problem,but i think my solution will …
Count intron retention reads from BAM files
written 8.5 years ago by Chirag Nepal 240

Hi there, How to find/filter intron-retained reads from BAM files? Which flag in BAM file describes information about retained intron. I know there are some …
Filtering a Tophat BAM file
written 11.3 years ago by Mehmet 82

Dear All, I have a BAM file after Tophat2. When I checked it on Artemis after gene prediction using exon, intron and exon-intron combined, I …
Output unmapped reads to a file when using tophat2
written 11.6 years ago by biolab 143

Hi, everyone, I am using tophat2 for reads mapping. I know there is a` --un` option in bowtie2 to write unmapped reads to a file. …
Blast Output With Species Name
written 12.8 years ago by Federico Giorgi 74

Dear all, my question is very simple. I have reads from a human RNASeq and I would like to check for contamination on a specific …
Extracting Exon/Intergenic/Intron Region Of Human Genome
written 14.7 years ago by Dataminer 277

Hi! Is their a way (already available in form of tool) to obtain all sequences or regions in human (hg 18) or mouse genome belonging …