SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa sampe segmentation fault papori Bioinformatics 5 09-22-2013 10:05 PM
BWA Alignment Segmentation Fault adkostic Bioinformatics 24 09-09-2013 12:52 AM
segmentation fault in BWA sampe papori Illumina/Solexa 0 07-28-2011 08:12 AM
bwa aln Segmentation fault DNAjunk Bioinformatics 4 03-02-2011 06:28 AM
BWA Segmentation Fault (aln) raela Bioinformatics 0 05-18-2010 06:41 AM

Reply
 
Thread Tools
Old 02-15-2012, 02:55 AM   #1
attilav
Junior Member
 
Location: Hungary

Join Date: Mar 2011
Posts: 7
Default BWA 0.6.1 fail to insert infer isize: weird pairing + Segmentation Fault

Hi all!

Using the command: bwa sampe refseq.fasta seq_1.aln.sai seq_2.aln.sai seq_1.fastq.gz seq_2.fastq.gz > output.sam , sometimes I get the error that you can see in the title. It looks like this:

[bwa_sai2sam_pe_core] 5505024 sequences have been processed.
[bwa_sai2sam_pe_core] convert to sequence coordinate...
[infer_isize] (25, 50, 75) percentile: (410, 437, 469)
[infer_isize] fail to infer insert size: weird pairing
[bwa_sai2sam_pe_core] time elapses: 0.82 sec
[bwa_sai2sam_pe_core] changing coordinates of 643 alignments.
[bwa_sai2sam_pe_core] align unmapped mate...
Segmentation fault

I did several very similar alignment tasks with BWA, most of them goes through fine, but with a few fastq files this happens. My data is from the 1000genomes, so that shouldn't be the problem I guess.
Please let me know, if you guys have an idea what's causing this!

Yours,
Attila Vikor
attilav is offline   Reply With Quote
Old 02-15-2012, 08:47 AM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

The simplest thing is to redo the alignment. If you mistyped a file name somewhere, you might get that result.

Also, try making the single end .sams separately, and spot check a few pairs of reads, see if thei rcoordiantes and relative orientation make sense. Weird pairing might mean that the reads don't run into each other like they should.
swbarnes2 is offline   Reply With Quote
Old 02-15-2012, 09:01 AM   #3
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

Is your genome index made with bwa 0.6.1?
How much memory do you have ?
Is the error consistent (reproduceable)?
Richard Finney is offline   Reply With Quote
Old 02-15-2012, 01:46 PM   #4
attilav
Junior Member
 
Location: Hungary

Join Date: Mar 2011
Posts: 7
Default

I made the genome index with bwa 0.6.1.
I have 2GB RAM on this computer, which is not optimal, but the segfaults only happen, if there is a weird pairing before.
Its consistent.

I will align the files separately, and see if something comes up.
attilav is offline   Reply With Quote
Old 02-15-2012, 02:56 PM   #5
Richard Finney
Senior Member
 
Location: bethesda

Join Date: Feb 2009
Posts: 700
Default

I think you're going to need more that 2GB.
Try one with 8GB on the problem samples.
Richard Finney is offline   Reply With Quote
Old 02-20-2012, 02:54 PM   #6
attilav
Junior Member
 
Location: Hungary

Join Date: Mar 2011
Posts: 7
Default

Apparently the problem was some of the files getting corrupted during downloading. Sorry to have wasted your time with this... Guess I have learned a lesson to always check the MD5 sums when downloading huge files from the internet.
attilav 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 10:18 AM.


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