fixed seq qual length check

[biopieces.git] / bp_bin / remove_indels

diff --git a/bp_bin/remove_indels b/bp_bin/remove_indels
index 59bbd4f92b932073a7fdaecd211d410712e18f7e..9a1ad033a0fbe508f4f12891ba06885179040814 100755 (executable)
--- a/bp_bin/remove_indels
+++ b/bp_bin/remove_indels
@@ -26,6 +26,7 @@
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
 
+use warnings;
 use strict;
 use Maasha::Biopieces;
 
@@ -42,7 +43,8 @@ $out = Maasha::Biopieces::write_stream( $options->{ "stream_out" } );
 
 while ( $record = Maasha::Biopieces::get_record( $in ) ) 
 {
-    $record->{ 'SEQ' } =~ tr/-~.//d if $record->{ "SEQ" };
+    $record->{ 'SEQ' }     =~ tr/-~.//d if $record->{ "SEQ" };
+    $record->{ 'SEQ_LEN' } = length $record->{ "SEQ" };
 
     Maasha::Biopieces::put_record( $record, $out );
 }