SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 12-22-2011, 05:19 PM   #1
gene coder
Member
 
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
raonyguimaraes
Member
 
Location: Belo Horizonte - Brazil

Join Date: Jun 2010
Posts: 38
Default

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

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: 699
Default

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

Join Date: Jul 2011
Posts: 18
Default

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
Reply

Tags
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 04:55 AM.


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