polished shred_seq

[biopieces.git] / bp_bin / shred_seq

diff --git a/bp_bin/shred_seq b/bp_bin/shred_seq
index 5b7f7ff7760b42bee078ee0c51abc86a770d65a4..aaff5b79dbcc52b56c57c8cd123f7f0e36fa25e9 100755 (executable)
--- a/bp_bin/shred_seq
+++ b/bp_bin/shred_seq
@@ -42,6 +42,9 @@ class Seq
   # seq.shred { |subseq| } -> Seq
   # seq.shred              -> []
   def shred(size, max_cov)
+    raise SeqError, "Size: #{size} < 1"            if size    < 1
+    raise SeqError, "Max coverage: #{max_cov} < 1" if max_cov < 1
+
     entries = []
     sum     = 0
     strand  = '+'