refactoring of assemble_pairs

[biopieces.git] / bp_bin / find_gaps

diff --git a/bp_bin/find_gaps b/bp_bin/find_gaps
index 300ae0952b9cb9c0c194eab4d3bf7087d490aff2..ade831ad9c9e1d28060050bec35a850a8cb4572a 100755 (executable)
--- a/bp_bin/find_gaps
+++ b/bp_bin/find_gaps
@@ -83,11 +83,13 @@ sub find_gaps
 
     # Returns a list.
 
-    my ( @gaps, $block, $beg, $end );
+    my ( $gap_char, @gaps, $block, $beg, $end );
+
+    $gap_char = 'N';
 
     $seq = uc $seq;
 
-    $block = 'N' x $min;
+    $block = $gap_char x $min;
 
     $beg = 0;
 
@@ -99,7 +101,7 @@ sub find_gaps
 
         $end = $beg;
 
-        while ( substr( $seq, $end, 1 ) eq 'N' ) {
+        while ( substr( $seq, $end, 1 ) eq $gap_char ) {
             $end++;
         }