added seed_size to soap_seq

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

diff --git a/code_perl/Maasha/Biopieces.pm b/code_perl/Maasha/Biopieces.pm
index ee38e282d4392a56f181713480cfe66b59d40576..80d1c18ae1dd2c216de777f62976211897fe45d1 100644 (file)
--- a/code_perl/Maasha/Biopieces.pm
+++ b/code_perl/Maasha/Biopieces.pm
@@ -574,6 +574,7 @@ sub get_options
         @options = qw(
             in_file|i=s
             genome|g=s
+            seed_size|s=s
             mismatches|m=s
             gap_size|G=s
             cpus|c=s
@@ -3582,6 +3583,7 @@ sub script_soap_seq
 
     my ( $genome, $tmp_in, $tmp_out, $fh_in, $fh_out, $record, $line, @fields, $entry );
 
+    $options->{ "seed_size" }  ||= 10;
     $options->{ "mismatches" } ||= 2;
     $options->{ "gap_size" }   ||= 0;
     $options->{ "cpus" }       ||= 1;
@@ -3604,7 +3606,7 @@ sub script_soap_seq
 
     close $fh_out;
 
-    Maasha::Common::run( "soap", "-r 2 -a $tmp_in -v $options->{ 'mismatches' } -g $options->{ 'gap_size' } -p $options->{ 'cpus' } -d $options->{ 'in_file' } -o $tmp_out > /dev/null 2>&1", 1 );
+    Maasha::Common::run( "soap", "-s $options->{ 'seed_size' } -r 2 -a $tmp_in -v $options->{ 'mismatches' } -g $options->{ 'gap_size' } -p $options->{ 'cpus' } -d $options->{ 'in_file' } -o $tmp_out > /dev/null 2>&1", 1 );
 
     unlink $tmp_in;