added check for minimum seq length

[biopieces.git] / bp_bin / assemble_pairs2

diff --git a/bp_bin/assemble_pairs2 b/bp_bin/assemble_pairs2
index 3741222a0f2a03e9983f1d390017d4870bd092e0..c7bfe1bc7007dbca8cd24c98d8c434b554585882 100755 (executable)
--- a/bp_bin/assemble_pairs2
+++ b/bp_bin/assemble_pairs2
@@ -55,7 +55,10 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
         entry2.reverse!.complement!
       end
 
-      if entry1 and entry2
+      if entry1 and
+         entry2 and
+         entry1.length >= options[:overlap_min] and
+         entry2.length >= options[:overlap_min]
         merged = Assemble.pair(
           entry1,
           entry2,