From: martinahansen Date: Thu, 17 Nov 2011 14:50:32 +0000 (+0000) Subject: fixed bug in remove_primers X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;h=bc41bf1f0f46fbf4c6468a0bc55bdb093c496084;p=biopieces.git fixed bug in remove_primers git-svn-id: http://biopieces.googlecode.com/svn/trunk@1659 74ccb610-7750-0410-82ae-013aeee3265d --- diff --git a/bp_bin/remove_primers b/bp_bin/remove_primers index a68fb89..4e4e2e7 100755 --- 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