Are there any aligners that enable me to align sense strands only to a sense strand reference sequence
1
0
Entering edit mode
7.0 years ago
b10hazard • 0

Is there an aligner that allows me to ignore reverse complement alignments. For example, my reference sequence I'm aligning to will always be single stranded sense DNA (the antisense strand is selectivily degraded). Also, the fragments I'm aligning are always going to be single stranded sense DNA strands. I was using BWA but there doesn't seem to be an option to ignore reverse complement alignments. Does any aligner out there do this? Thanks!

next-gen alignment bwa sequencing • 2.4k views
ADD COMMENT
2
Entering edit mode
7.0 years ago
valdirbarth ▴ 20

Bowtie seems to be able to do it with this parameter:

--nofw/--norc If --nofw is specified, bowtie will not attempt to align against the forward reference strand. If --norc is specified, bowtie will not attempt to align against the reverse-complement reference strand. For paired-end reads using --fr or --rf modes, --nofw and --norc apply to the forward and reverse-complement pair orientations. I.e. specifying --nofw and --fr will only find reads in the R/F orientation where mate 2 occurs upstream of mate 1 with respect to the forward reference strand.

http://bowtie-bio.sourceforge.net/manual.shtml

ADD COMMENT
0
Entering edit mode

Just what I needed. Thanks!

ADD REPLY

Login before adding your answer.

Traffic: 2133 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6