From 94c5176532db69ab360414e082111f1d8e47a925 Mon Sep 17 00:00:00 2001 From: martinahansen Date: Wed, 13 Mar 2013 09:59:23 +0000 Subject: [PATCH] added check for minimum seq length git-svn-id: http://biopieces.googlecode.com/svn/trunk@2140 74ccb610-7750-0410-82ae-013aeee3265d --- bp_bin/assemble_pairs2 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/bp_bin/assemble_pairs2 b/bp_bin/assemble_pairs2 index c7bfe1b..611ca56 100755 --- 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, -- 2.39.2