Question: Genomic Bam To Transcriptome Bam
1
gravatar for dfernan
5.4 years ago by
dfernan640
United States
dfernan640 wrote:

Hi,

I was looking for a tool to convert a bam file in genomic coordinates to transcriptomic coordinates (using a GTF reference). Is there a fast way to do this for large bam files? Such tool would filter out reads that do not map to any transcript.

Let me know if you have any suggestions,

Thanks,

Daniel F.

bam samtools • 3.0k views
ADD COMMENTlink modified 3.4 years ago by ngsbioinformatics30 • written 5.4 years ago by dfernan640
2

You'll probably get the best answer by simply mapping your reads to the transcripts directly if having transcript-based BAM files is really your goal. If, instead, you want read counts for your genomic BAM files, that can be done with several tools. If the latter is what you want, then perhaps you can clarify your question.

ADD REPLYlink written 5.4 years ago by Sean Davis25k

Hi Sean, thanks. The problem is that if I map directly to the transcripts I loose the ability to build new transcripts, find splice junctions, etc. Also in many cases I am dealing with totalRNA experiments so mapping only to transcripts I loose many reads. Let us say I already mapped all my data and just want to transform coordinates, there is no way to do this? Maybe you are right and it does not make much sense that's why it does not exist...

ADD REPLYlink written 5.4 years ago by dfernan640
1
gravatar for Biomonika (Noolean)
5.4 years ago by
State College, PA, USA
Biomonika (Noolean)3.0k wrote:

I believe that you may use bedtools intersect:

http://bedtools.readthedocs.org/en/latest/content/tools/intersect.html

ADD COMMENTlink modified 5.4 years ago • written 5.4 years ago by Biomonika (Noolean)3.0k
0
gravatar for ngsbioinformatics
3.4 years ago by
United States
ngsbioinformatics30 wrote:

I believe you want this tool sam-xlate from UBU: https://github.com/mozack/ubu/wiki

It doesn't take a GTF, but rather a BED file in refFlat (I believe) format.

ADD COMMENTlink written 3.4 years ago by ngsbioinformatics30
Please log in to add an answer.

Help
Access

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