From 9cd6824c09e22e968d8b2cac68a79f0b95a274d3 Mon Sep 17 00:00:00 2001
From: martinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Date: Thu, 24 Jul 2008 07:49:26 +0000
Subject: [PATCH] added seed_size to soap_seq

git-svn-id: http://biopieces.googlecode.com/svn/trunk@192 74ccb610-7750-0410-82ae-013aeee3265d
---
 code_perl/Maasha/Biopieces.pm | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/code_perl/Maasha/Biopieces.pm b/code_perl/Maasha/Biopieces.pm
index ee38e28..80d1c18 100644
--- 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;
 
-- 
2.39.5