.bam files to fastq including the unmapped reads · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

.bam files to fastq including the unmapped reads

Hello,

I have tophat output (.bam) generated by tophat and I want to convert them back to fastq files. Since it was mapped using tophat, I thought that the best tool to do that will be: bam2fastx. Though, I want to refer both the mapped reads (tophat output: accepted_hits.bam) and the unmapped ones (tophat output: unmapped.bam). How can it be done?

Thanks!!

rna-seq rna-seq alignment

• 3,241 views • link

updated 6.5 years ago by Biostar

• written 10.2 years ago by galozs •

2 answers

You can convert accepted_hits.bam and unmapped.bam separately and then just merge the two FASTQs. It's just a text file, so the order does not matter.

• 1 views • link

written 10.2 years ago by igor

If I only use the accepted reads (and afterwards turn it again to bam during the analysis) I can loose some information right? I wonder why there is no feature for doing this step backwards without loosing information...

• 0 views • link

written 10.2 years ago by galozs •

Try bedtools bamtofastq. I have recently used this with single end reads from RNA-Seq, so not sure how well it works if you have paired-end reads.

• 1 views • link

written 10.2 years ago by mastal511

Log in to answer this question.

More posts like this

Forum: chimeric vs unaligned reads
written 4.8 years ago by Galileo •

Hello, I'm trying to obtain the chemic alignments from a BAM file that originally was generated by STAR and that only has a list of …
from unmapped bam file picard samtofastq produced no output
written 6.4 years ago by Rashedul Islam

I generated bam of unmapped reads using: >$ samtools view -b -f 4 input.bam >unmapped.bam > >$ samtools flagstat unmapped.bam > > 2313037 + 4694809 …
Getting unmapped paired-end reads as fastq from unmapped.bam obtained with tophat2
written 7.5 years ago by cris.kgs •

Hi, I used tophat2 to align paired-end reads to Arabidopsis genome and I want to use the unmapped.bam file to obtain the unmapped reads so …
Realign Aligned.sortedByCoord.out.bam to Aligned.toTranscriptome.out.bam for RSEM input
written 7.7 years ago by mkvalluru •

Hi, I was using STAR-RSEM to quantify isoforms (primary input fastq files). Recently I have received Aligned.sortedByCoord.out.bam from the collaborators. However, I can't use these …
extract one-side unmapped reads from sam
written 8.2 years ago by liorglic

Hello, I've ran [TopHat][1] and aligned paired end (PE) reads to a reference genome. TopHat outputs two bam files: accepted_hits.bam and unmapped.bam. I'd like to …
How to change the quality score of reads ? (quick help is appreciated)
written 10.3 years ago by mbparsa •

So I have my FASTQ ran through Tophat and produced unmapped and mapped bam files, I would like to increase the quality score for all …
[Solved] Problems merging TopHat mapped and unmapped BAMs for submission to ENA
written 11.0 years ago by Chris Cole

<p>As part of their submission requirements ENA/ArrayExpress require BAMs to include all reads (e.g. also unmapped ones). TopHat in its wisdom splits unmapped reads from …
Tophat more multiple alignments than mapped reads?
written 11.7 years ago by bmpbowen •

I used bam2fastx to convert the unmapped.bam output from TopHat into new fastq files. Then I re-ran Tophat with those fastq files against a bacterial …
Why # Of Reads From Accepted_Hits.Bam + Unmapped.Bam > # Of Reads In Fastq File?
written 13.0 years ago by newDNASeqer

<p>a quick question:</p> <p>After running Tophat with a fastQ files, I found the # of reads from (accepted_hits.bam and unmapped.bam) is greater than the # …
What Kind Of Inputs Does Cuffdiff Take?
written 13.5 years ago by narges •

<p>Hi, </p> <p>I have run tophat over the RNA-seq data and I got the output files : accepted_hits.bam, insertions.bed, junctions.bed, deletions.bed, unmapped.bam.</p> <p>Now to find …