Reverse engineering BAM files: BAM -> FASTQ

gene coder · 12-22-2011, 05:19 PM

How can I possibly extract the reads from a BAM file and put them into a FASTQ file for simulation (maq simutrain, then maq simulate)?

Should I just extract col. 1, 10 and 11 from a BAM file and put them in a text file along with a '+'? That is the output by read simulators.

But does it not happen that DNA sequences and base quality sequences are reversed and/or transformed depending on the direction and strand of reads relative to reference genomes? I would have to fix that too in that case.

raonyguimaraes · 12-22-2011, 05:57 PM

http://biostar.stackexchange.com/que...g-bam-to-fastq

Richard Finney · 12-22-2011, 06:39 PM

Check out http://seqanswers.com/forums/showthread.php?t=16395 for bampe2fq.c and bamse2fq.c for fast implementations of bam to fastq programs. It handles your concerns about reverse complimenting the sequence and reversing the quality string.

gene coder · 01-03-2012, 03:42 PM

Thanks. I found that SamToFastq in Picard did the job on a chromosome of NA12878 from the 1000 Genomes Project.

1. I separated SE reads and PE reads by library into separate BAM files using SamTools.
2. For PE reads, I also had to get rid of read-pairs that were unmapped (bit 3) or whose mate was unmapped (bit 4) or that were not properly aligned (bit 2).
3. I further separated the BAM files by read length.
4. Each BAM file now contained SE or PE reads only of the same read length and the same library.
5. Now I could run SamToFastq to convert a SAM file to a DAT file for use with MAQ simutrain.

Those were the steps to do to get what I wanted.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Convert BAM file to FASTQ	malachig	Bioinformatics	28	02-02-2016 04:04 AM
Convert merged BAM back to per lane BAM or FASTQ file	danielsbrewer	Bioinformatics	6	10-03-2013 08:29 AM
NEw to Chip-seq and have .bam/.sam/.bam.bai files... then what?	NGS newbie	Bioinformatics	11	05-25-2011 08:48 AM
Consensus part from sequence read(fastq) and align(BAM) files	culmen	Bioinformatics	5	12-21-2010 04:57 AM
Stand Alone Bam to FASTQ	dcfargo	Bioinformatics	11	08-02-2010 06:56 AM

12-22-2011, 05:19 PM	#1
gene coder Member Location: Ireland Join Date: Jul 2011 Posts: 18	Reverse engineering BAM files: BAM -> FASTQ How can I possibly extract the reads from a BAM file and put them into a FASTQ file for simulation (maq simutrain, then maq simulate)? Should I just extract col. 1, 10 and 11 from a BAM file and put them in a text file along with a '+'? That is the output by read simulators. But does it not happen that DNA sequences and base quality sequences are reversed and/or transformed depending on the direction and strand of reads relative to reference genomes? I would have to fix that too in that case.

12-22-2011, 05:57 PM	#2
raonyguimaraes Member Location: Belo Horizonte - Brazil Join Date: Jun 2010 Posts: 38	http://biostar.stackexchange.com/que...g-bam-to-fastq

12-22-2011, 06:39 PM	#3
Richard Finney Senior Member Location: bethesda Join Date: Feb 2009 Posts: 700	Check out http://seqanswers.com/forums/showthread.php?t=16395 for bampe2fq.c and bamse2fq.c for fast implementations of bam to fastq programs. It handles your concerns about reverse complimenting the sequence and reversing the quality string.

01-03-2012, 03:42 PM	#4
gene coder Member Location: Ireland Join Date: Jul 2011 Posts: 18	Thanks. I found that SamToFastq in Picard did the job on a chromosome of NA12878 from the 1000 Genomes Project. 1. I separated SE reads and PE reads by library into separate BAM files using SamTools. 2. For PE reads, I also had to get rid of read-pairs that were unmapped (bit 3) or whose mate was unmapped (bit 4) or that were not properly aligned (bit 2). 3. I further separated the BAM files by read length. 4. Each BAM file now contained SE or PE reads only of the same read length and the same library. 5. Now I could run SamToFastq to convert a SAM file to a DAT file for use with MAQ simutrain. Those were the steps to do to get what I wanted.