finishing slice_align biopiece

[biopieces.git] / bp_bin / slice_align

diff --git a/bp_bin/slice_align b/bp_bin/slice_align
index 9c8ef6d034a112048e5351924e5ce7bb7427abf1..ed561981869c7549acfd4e5cea8e0d1b0fd5d6c5 100755 (executable)
--- a/bp_bin/slice_align
+++ b/bp_bin/slice_align
@@ -110,6 +110,7 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
       entry = Seq.new(seq: record[:SEQ])
 
       unless options[:beg]
+        raise "template length != alignment length" if template and template.length != entry.length
         compact = template ? template : Seq.new(seq: entry.seq.dup)
         compact.seq.delete! "-.~"