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 c7bfe1bc7007dbca8cd24c98d8c434b554585882..611ca56609c63b236722ded3769725f94675e8e2 100755 (executable)
--- a/bp_bin/assemble_pairs2
+++ b/bp_bin/assemble_pairs2
@@ -51,14 +51,15 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
         entry1 = Seq.new_bp(record)
       elsif entry2.nil?
         entry2 = Seq.new_bp(record)
-        entry2.type = :dna
-        entry2.reverse!.complement!
       end
 
       if entry1 and
          entry2 and
          entry1.length >= options[:overlap_min] and
          entry2.length >= options[:overlap_min]
+        entry2.type = :dna
+        entry2.reverse!.complement!
+
         merged = Assemble.pair(
           entry1,
           entry2,