finishing slice_align biopiece

git-svn-id: http://biopieces.googlecode.com/svn/trunk@2269 74ccb610-7750-0410-82ae-013aeee3265d

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]
       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! "-.~"
 
         compact = template ? template : Seq.new(seq: entry.seq.dup)
         compact.seq.delete! "-.~"