SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SAM flag forward/reverse first/second nupurgupta Bioinformatics 3 01-20-2016 03:47 AM
rmdup can not move duplicates in forward and reverse strand for single-end reads ct586 Bioinformatics 4 03-11-2012 05:01 PM
how to define a forward or reverse read file poorphd Illumina/Solexa 3 11-22-2011 12:34 PM
How to separate coverage of forward and reverse reads on same axis? Kennels Bioinformatics 6 05-04-2011 12:12 AM
forward and reverse sequance asankaf General 5 05-27-2009 07:48 AM

Reply
 
Thread Tools
Old 02-20-2013, 05:39 AM   #1
jwhite
Member
 
Location: Boston

Join Date: Jun 2012
Posts: 33
Default forward and reverse reads with BWA

Hi,

For those of you who use BWA:

Does BWA treat forward reads different from reverse reads?

The reason I ask is to determine whether the forward and reverse reads can be concatenated into a single file to supply to BWA.

Joe White
jwhite is offline   Reply With Quote
Old 02-20-2013, 05:50 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

By forward and reverse reads, are you referring to paired-end reads or reads which will map to the forward and reverse strand of the reference sequence? Presuming the former, don't concatenate them together. See, instead, the sampe command.
dpryan is offline   Reply With Quote
Old 02-20-2013, 05:57 AM   #3
jwhite
Member
 
Location: Boston

Join Date: Jun 2012
Posts: 33
Default

Quote:
Originally Posted by dpryan View Post
By forward and reverse reads, are you referring to paired-end reads or reads which will map to the forward and reverse strand of the reference sequence? Presuming the former, don't concatenate them together. See, instead, the sampe command.
I was indeed referring to the former. A colleague is doing this, and I thought it unusual, so wanted to ask.

Thanks.
jw
jwhite is offline   Reply With Quote
Old 02-20-2013, 08:41 AM   #4
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Concatenating them together is almost certainly the wrong thing to do.

For one thing, there is some unknown amount of sequence between the two ends of a paired end read. You might have some empirical data about what the likely size distribution of your library is, but that won't tell you how far apart any particular pair of reads is.
swbarnes2 is offline   Reply With Quote
Reply

Tags
sequence alignment

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 01:13 PM.


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