SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa samse -s parameter rahul Bioinformatics 6 03-01-2012 03:44 PM
BWA - samse giverny Bioinformatics 6 07-01-2010 07:09 AM
Why no multithreading for BWA sampe/samse? krobison Bioinformatics 5 02-20-2010 11:14 AM
BWA samse -n option kellyv Bioinformatics 0 01-19-2010 12:15 PM
BWA samse results combiochem Bioinformatics 1 11-16-2009 11:22 AM

Reply
 
Thread Tools
Old 08-28-2014, 11:13 PM   #1
scami
Member
 
Location: italy

Join Date: Sep 2010
Posts: 55
Default Bwa, sampe vs samse

Hi guys,

this is probably a silly question. If I run bwa sampe I will get an alignment file which report informations about the paired end reads being in the expected pairing mode on the reference genome. However also not properly mapped paired reads will be reported, such as reads that maps although the mate does not map. Does this mean that the alignment sam file will provide the same information reported by the sam file generated by bwa sampe?
If I run sampe with paired end reads I am treating all the reads as unpaired and I will map them on the reference. Since for samse I will map the reads even if they are unpaired I would expect the same result than samse. Am I correct?
thanks
scami is offline   Reply With Quote
Old 08-29-2014, 12:48 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

I'm guessing that you swapped "sampe" and "samse" a few times in your post.

In general, if a pair are aligned as singletons with sampe then they should be given the same alignment with samse. The only difference with by the FLAG field, since in one case they would have been paired in sequencing (bit 0x1) and have had an order (bits 0x40 and 0x80), though obviously bit 0x2 would be unset.
dpryan is offline   Reply With Quote
Old 08-29-2014, 01:16 AM   #3
scami
Member
 
Location: italy

Join Date: Sep 2010
Posts: 55
Default Thanks

Thanks a lot! Yes sorry I used sampe and samse all the way around. You answer clarifies everything.
Cheers
scami 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 06:53 PM.


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