Question: Paired End Mapping : What Is Bwa Solid Paired End Default Direction (Bwa Sampe) ?
5
gravatar for Hx
7.9 years ago by
Hx60
Hx60 wrote:

Hi,

SOLiDTM Library construct

For Paired End

F3 F5

---------><----------

Can anyone tell me the default directions on how BWA will map paired ends?

Thanks in advance.

paired solid bwa • 4.6k views
ADD COMMENTlink written 7.9 years ago by Hx60
4
gravatar for brentp
7.9 years ago by
brentp23k
Salt Lake City, UT
brentp23k wrote:

BWA expects mate-pair F3/R3 (as input to solid2fastq.pl).

You could reverse-complement your F5 reads (and reverse the quality) and send them to solid2fastq.pl and it should work more/less as expected. Then you just have to consider how having the low-quality bases at the start of the read affects BWA (lh3?) -- you may want to quality trim your reads before-hand.

EDIT: if this works, it would be just reverse, not reverse-complement of the F5 reads.

ADD COMMENTlink modified 7.3 years ago • written 7.9 years ago by brentp23k

Thanks brentp, solve my question mark on weird flagstats. However, regarding "reverse-complement your F5 reads (and reverse the quality)" as you mentioned, can I know how to reverse the reads?

ADD REPLYlink written 7.9 years ago by Hx60

just reverse complement the reads. You can do this with bioperl or biopython. In perl, it's like my $rc = reverse(seq); =~ tr/ACGTacgt/TGCAtgca/. So just loop over the csfasta and do that to every other line. And for the qual file. Just reverse it.

ADD REPLYlink written 7.9 years ago by brentp23k

SOLiD reads are colourspace reads where I cant perform reverse complement on them, are there any other ways?

ADD REPLYlink written 7.9 years ago by Hx60

Oh, yes of course. I think in colorspace, you can just reverse the colorspace sequence, without complementing to get (what would be in base-space) the reverse complement. I havent mapped F3/F5 Solid reads but maybe check the -c parameter to bwa (http://bio-bwa.sourceforge.net/bwa.shtml)

ADD REPLYlink modified 7.1 years ago by Istvan Albert ♦♦ 80k • written 7.9 years ago by brentp23k

I'm not sure whether this works, will give a try, thanks for the suggestions brentp, but I'm still hoping for a more concrete solution, perhaps other software may help

ADD REPLYlink written 7.9 years ago by Hx60
0
gravatar for Yaira04
7.7 years ago by
Yaira040
Yaira040 wrote:

Hi HX,

did you find any better solution eventually? I am confronting the same issue these days...

Thanks.

ADD COMMENTlink written 7.7 years ago by Yaira040

Yaira, it seems there's no concrete solution for bwa paired end mapping yet since bwa assumes paired end mapping has the same orientation as mate pair, maybe you can consider trying out Bioscope or try using bfast http://seqanswers.com/forums/showthread.php?t=7100 or bowtie http://bowtie-bio.sourceforge.net/manual.shtml#paired-end-colorspace-alignment

ADD REPLYlink modified 6.9 years ago by Istvan Albert ♦♦ 80k • written 7.3 years ago by Hx60
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1461 users visited in the last hour