Question: Get introns in bam format from alignment reads result
gravatar for xd_d
3.7 years ago by
xd_d110 wrote:

Hey all,

I have alignt rna-seq data with Star against the reference genome hg38. Now I have the results in bam format.

Does anyone know if there is a program that can that detect the introns from the alignment result. Finally I want the only introns in a bed file that are mapped against the reference genome. Before I start programing I want ask if somebody know about nice tools ?

rna-seq intron • 1.0k views
ADD COMMENTlink written 3.7 years ago by xd_d110

What is your actual end goal?

ADD REPLYlink written 3.7 years ago by Devon Ryan97k

I usually extract reads mapped to exons from BAM files with featureCounts using a gtf file from ENSEMBL. So my guess is that you can also extract the reads that map to the introns only with featureCounts.

ADD REPLYlink written 3.7 years ago by Benn8.0k

This question is not entirely clear if you are just looking get known exons as a BED file or are looking to find new exons that may be represented in your file file. If it is the former then these directions are still valid. Otherwise you are going to need to look for split-read mappings in your file.

ADD REPLYlink written 3.7 years ago by genomax92k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1020 users visited in the last hour