SEQanswers

Go Back   SEQanswers > Applications Forums > Genomic Resequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA and sub-optimal alignments LauraR Bioinformatics 0 02-10-2012 07:46 AM
multiple alignments with BWA efoss Bioinformatics 0 08-07-2011 11:37 AM
BWA: suboptimal alignments ElMichael Bioinformatics 2 03-29-2011 01:31 PM
BWA alignments and time cggj Bioinformatics 5 05-13-2010 10:41 AM
Merging alignments in mosaik - problems 454andSolid Bioinformatics 0 05-13-2010 08:44 AM

Reply
 
Thread Tools
Old 06-24-2010, 12:08 PM   #1
arnkas
Junior Member
 
Location: Seattle WA

Join Date: Nov 2009
Posts: 5
Default merging bwa alignments

Hello,
I am using bwa to map solid data to the mouse genome, but get a segmentation fault while running bwa aln. If I separate the genome into 2 halves, it works all right, so I assume that the segfault was a memory problem. My question is: Is there a way to merge the result of mapping to 2 halves of the genome?

I have a similar question about dividing the reads into smaller subsets, while using the full database.

I appreciate any help on this.
arnkas is offline   Reply With Quote
Old 07-06-2010, 12:29 PM   #2
Lee Sam
Member
 
Location: Ann Arbor, MI

Join Date: Oct 2008
Posts: 57
Default

I've been doing BWA alignments and merging them using picard-tools "MergeSamFiles.jar"
Lee Sam is offline   Reply With Quote
Old 07-06-2010, 05:00 PM   #3
arnkas
Junior Member
 
Location: Seattle WA

Join Date: Nov 2009
Posts: 5
Default merging files which map reads to separate parts of the genome

Thank you for you response. I have also used samtools merge.

As I understand it, bwa eliminates reads that map to more than one place in the genome. Since my 2 files involves mapping reads to separate halves of the genome, reads that are mapped in both files should be eliminated. Does the merge process recognize this?
arnkas is offline   Reply With Quote
Old 07-06-2010, 06:23 PM   #4
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

No it does not.
nilshomer is offline   Reply With Quote
Old 07-07-2010, 05:20 AM   #5
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

One should never split the genome. Split the read set. It is true for essentially all NGS aligners.
lh3 is offline   Reply With Quote
Old 07-07-2010, 07:26 AM   #6
arnkas
Junior Member
 
Location: Seattle WA

Join Date: Nov 2009
Posts: 5
Default

Well, the reason I split the genome was that I ran into a segmentation fault running "bwa aln -c" against the entire mouse genome. I was able to isolate this to a single read. I tried it against various genomes, and it worked all right for smaller genomes, eg human mrna and half the mouse genome.

Here is the problem read:
@559_9_257_F3
GNANNNGGGGGNNNNNCNNNAGCACGAGAANNNNCGNGNGNTNCNNCNN
+
2!'!!!)*'&)!!!!!&!!!')&&.*)&+'!!!!(,!&!)!'!&!!&!!

I also assumed that there must be a way to postprocess a merged file from a split genome, although I haven't done it yet.
arnkas is offline   Reply With Quote
Reply

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 08:23 AM.


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