fixed bug in remove_primers

[biopieces.git] / bp_bin / remove_primers

diff --git a/bp_bin/remove_primers b/bp_bin/remove_primers
index a68fb89ce732b4bbf0e53b3ab8637d6aca6d18f8..4e4e2e7fe48abdb101f4de2ae0a7e59e1097a1d3 100755 (executable)
--- a/bp_bin/remove_primers
+++ b/bp_bin/remove_primers
@@ -56,10 +56,10 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
       end
 
       if reverse = seq.patscan(options[:reverse].to_s, pos, options[:mismatches], options[:insertions], options[:deletions])
-        record[:REVERSE_POS] = reverse.last.pos
-        record[:REVERSE_LEN] = reverse.last.length
-        pos = reverse.last.pos
-        seq.subseq!(pos)
+        record[:REVERSE_POS] = reverse.first.pos
+        record[:REVERSE_LEN] = reverse.first.length
+        pos = reverse.first.pos
+        seq.subseq!(0, pos)
       end
 
       if pos > 0