refactoring of revcomp in seq.rb

[biopieces.git] / bp_bin / shred_seq

diff --git a/bp_bin/shred_seq b/bp_bin/shred_seq
index 5b2a886af36e0708ebf4066795e0396067eae86b..465f08855bec36eb769bf319e34b5c5ff25d0d61 100755 (executable)
--- a/bp_bin/shred_seq
+++ b/bp_bin/shred_seq
@@ -51,7 +51,7 @@ class Seq
 
     while coverage(sum) < max_cov
       entry = self.subseq_rand(size)
-      entry.revcomp if strand == '-'
+      entry.reverse!.complement! if strand == '-'
 
       if block_given?
         yield entry