refactoring of revcomp in seq.rb

[biopieces.git] / bp_bin / find_adaptor

diff --git a/bp_bin/find_adaptor b/bp_bin/find_adaptor
index 4fd559dd281b45b06b372134e6215cc600ef55a2..ac3681a090c8e80c6139a4d1640ce4dd97664430 100755 (executable)
--- a/bp_bin/find_adaptor
+++ b/bp_bin/find_adaptor
@@ -53,11 +53,11 @@ casts << {:long=>'deletions',   :short=>'d', :type=>'uint',   :mandatory=>false,
 options = Biopieces.options_parse(ARGV, casts)
 
 if options[:forward_rc]
-  options[:forward] = Seq.new("test", options[:forward_rc], 'dna').revcomp.seq
+  options[:forward] = Seq.new("test", options[:forward_rc], 'dna').reverse.complement.seq
 end
 
 if options[:reverse_rc]
-  options[:reverse] = Seq.new("test", options[:reverse_rc], 'dna').revcomp.seq
+  options[:reverse] = Seq.new("test", options[:reverse_rc], 'dna').reverse.complement.seq
 end
 
 raise ArgumentError, "no adaptor specified" unless options[:forward] or options[:reverse]