Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA sampe: wierd pairing zlu Bioinformatics 18 11-10-2015 02:05 AM
bwa sampe error Mali Salmon Bioinformatics 14 10-27-2014 11:25 AM
BWA 0.6.1 fail to insert infer isize: weird pairing + Segmentation Fault attilav Bioinformatics 5 02-20-2012 02:54 PM
bwa sampe problem sszelinger Bioinformatics 3 06-15-2011 05:34 AM
bwa sampe max insert size zlu Bioinformatics 0 10-27-2009 07:35 AM

Thread Tools
Old 06-14-2011, 07:16 AM   #1
Location: california

Join Date: Jan 2011
Posts: 26
Question BWA - sampe max occurrences of a read for pairing


Can anyone explain what the parameter -o means for sampe.According to BWA manual it says "A read with more occurrences will be treated as a single-end read.Reducing this parameter helps faster pairing".

Does it mean uneven number of reads in the two input fastq files?How do you decide what parameter to use?

Any help will be greatly appreciated.

madsaan is offline   Reply With Quote
Old 06-14-2011, 01:48 PM   #2
Location: Brock University

Join Date: Dec 2010
Posts: 66

BWA determines the suffix array intervals for each read during the 'aln' stage. If you're allowing errors (bwa aln -n #oferrors), a read can have many suffix array intervals of varying lengths.

The more (and larger) suffix array intervals a read has, the more work BWA has to go through during 'sampe' to pair it with its mate as there is a greater combinations of positions to consider for each mate... if the number of potential alignments exceeds the -o parameter (defaulted at 100 000), BWA just skips pairing altogether and leaves each mate aligned as if they were not mates at all in order to maintain an efficient alignment process.
dp05yk is offline   Reply With Quote

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 03:28 AM.

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