Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Stand Alone Bam to FASTQ dcfargo Bioinformatics 14 08-20-2021 01:47 AM
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

Thread Tools
Old 12-22-2011, 05:19 PM   #1
gene coder
Location: Ireland

Join Date: Jul 2011
Posts: 18
Default 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.
gene coder is offline   Reply With Quote
Old 12-22-2011, 05:57 PM   #2
Location: Belo Horizonte - Brazil

Join Date: Jun 2010
Posts: 38

raonyguimaraes is offline   Reply With Quote
Old 12-22-2011, 06:39 PM   #3
Richard Finney
Senior Member
Location: bethesda

Join Date: Feb 2009
Posts: 700

Check out 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.
Richard Finney is offline   Reply With Quote
Old 01-03-2012, 03:42 PM   #4
gene coder
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.
gene coder is offline   Reply With Quote

maq, simulation

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 01:06 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO