SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 06-14-2011, 08:16 AM   #1
madsaan
Member
 
Location: california

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

Hello,

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.

Thanks,
Joji
madsaan is offline   Reply With Quote
Old 06-14-2011, 02:48 PM   #2
dp05yk
Member
 
Location: Brock University

Join Date: Dec 2010
Posts: 66
Default

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
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 11:52 AM.


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