X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Fassemble_pairs2;h=f6d2f28594056689190bbb5e34d66ae183132836;hb=dd65a9e6c523ee61c85ff0bfe00eaf1cc80be4c6;hp=bd3c7fac3e18792c55762200804f1bc84ce3bdd4;hpb=d83e1ce27879cdb71a8cad530e138d96214619ae;p=biopieces.git diff --git a/bp_bin/assemble_pairs2 b/bp_bin/assemble_pairs2 index bd3c7fa..f6d2f28 100755 --- a/bp_bin/assemble_pairs2 +++ b/bp_bin/assemble_pairs2 @@ -64,11 +64,22 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output| overlap_max:options[:overlap_max] ) - output.puts merged.to_bp if merged + if merged + new_record = merged.to_bp + + if merged.seq_name =~ /overlap=(\d+):hamming=(\d+)$/ + new_record[:OVERLAP] = $1 + new_record[:HAMMING_DIST] = $2 + end + + output.puts new_record + end entry1 = nil entry2 = nil end + else + output.puts record end end end