changeset 2:35ac9de60bec draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meryl commit 9c098f0bdeca38de2f282a85a20e9e6fb3fd4929"
author iuc
date Fri, 09 Apr 2021 13:21:51 +0000
parents e5b18909f73d
children 9f7a227436d6
files macros.xml meryl.xml
diffstat 2 files changed, 21 insertions(+), 11 deletions(-) [+]
line wrap: on
line diff
--- a/macros.xml	Thu Apr 08 19:22:17 2021 +0000
+++ b/macros.xml	Fri Apr 09 13:21:51 2021 +0000
@@ -1,6 +1,6 @@
 <macros>
     <token name="@TOOL_VERSION@">1.2</token>
-    <token name="@GALAXY_TOOL_VERSION@">galaxy0</token>
+    <token name="@GALAXY_TOOL_VERSION@">galaxy1</token>
     <xml name="edam_ontology">
         <edam_topics>                                                                                  
             <edam_topic>topic_0196</edam_topic>
--- a/meryl.xml	Thu Apr 08 19:22:17 2021 +0000
+++ b/meryl.xml	Fri Apr 09 13:21:51 2021 +0000
@@ -8,16 +8,20 @@
     <version_command>meryl --version</version_command>
     <command detect_errors='exit_code'><![CDATA[
         #if $operation_type.command_type == 'count-kmers'
-            meryl
             #if $operation_type.options_kmer_size.kmer_size == 'estimate'
                 #from math import log
-                $operation_type.options_kmer_size.input_kmer_size = int(log(int($operation_type.options_kmer_size.genome_size)*(1-float($operation_type.options_kmer_size.collision_rate))/float($operation_type.options_kmer_size.collision_rate))/log(4))
+                #set size=int(log(int($operation_type.options_kmer_size.genome_size)*(1-float($operation_type.options_kmer_size.collision_rate))/float($operation_type.options_kmer_size.collision_rate))/log(4))
+            #elif $operation_type.options_kmer_size.kmer_size == 'provide'
+                #set size=$operation_type.options_kmer_size.input_kmer_size
             #end if
+            meryl
             $operation_type.count_operations
-            k=$operation_type.options_kmer_size.input_kmer_size
+            k=$size
+            memory=\${GALAXY_MEMORY_GB:-1}
+            threads=\${GALAXY_SLOTS:-1}
             $operation_type.input_reads
             output read-db.meryl &&
-            echo 'K-mer size: ${$operation_type.options_kmer_size.input_kmer_size}' &&
+            echo 'K-mer size: ${size}' &&
             tar -zcf read-db.meryldb read-db.meryl
         #elif $operation_type.command_type == 'filter-kmers'
             mkdir  -p ./temp_db/ &&
@@ -58,19 +62,25 @@
             meryl histogram tmp.meryl > read-db.hist
         #elif $operation_type.command_type == 'trio-mode'
             export MERQURY=\$(dirname \$(command -v merqury.sh))/../share/merqury/ &&
+            #if $operation_type.options_kmer_size.kmer_size == 'estimate'
+                #from math import log
+                #set size=int(log(int($operation_type.options_kmer_size.genome_size)*(1-float($operation_type.options_kmer_size.collision_rate))/float($operation_type.options_kmer_size.collision_rate))/log(4))
+            #elif $operation_type.options_kmer_size.kmer_size == 'provide'
+                #set size=$operation_type.options_kmer_size.input_kmer_size
+            #end if
             #for $i, $read in enumerate($paternal_reads):
                 mkdir 'paternal{$i}.meryl' &&
-                meryl count k=$operation_type.options_kmer_size.input_kmer_size '${read}' output 'paternal{$i}.meryl' &&
+                meryl count k=$size memory=\${GALAXY_MEMORY_GB:-1} threads=\${GALAXY_SLOTS:-1} '${read}' output 'paternal{$i}.meryl' &&
             #end for
             meryl union-sum paternal*.meryl output pat.meryl &&
             #for $i, $read in enumerate($maternal_reads):
                 mkdir 'maternal{$i}.meryl' &&
-                meryl count k=$operation_type.options_kmer_size.input_kmer_size '${read}' output 'maternal{$i}.meryl' &&
+                meryl count k=$size memory=\${GALAXY_MEMORY_GB:-1} threads=\${GALAXY_SLOTS:-1} '${read}' output 'maternal{$i}.meryl' &&
             #end for
             meryl union-sum maternal*.meryl output mat.meryl &&
             #for $i, $read in enumerate($child_reads):
                 mkdir 'child{$i}.meryl' &&
-                meryl count k=$operation_type.options_kmer_size.input_kmer_size '${read}' output 'child{$i}.meryl' &&
+                meryl count k=$size memory=\${GALAXY_MEMORY_GB:-1} threads=\${GALAXY_SLOTS:-1} '${read}' output 'child{$i}.meryl' &&
             #end for
             meryl union-sum child*.meryl output child.meryl &&
             ## mat specific kmers
@@ -117,7 +127,7 @@
             meryl union-sum output child.inherited.meryl mat.inherited.meryl pat.inherited.meryl shared.inherited.meryl &&
             meryl difference output read.only.meryl child.meryl child.inherited.meryl &&
             tar -czf 'read-db.meryldb' read.only.meryl &&
-            echo 'K-mer size: ${$operation_type.options_kmer_size.input_kmer_size}'
+            echo 'K-mer size: ${size}'
         #end if
     ]]>    </command>
     <inputs>
@@ -143,7 +153,7 @@
                         <option value="estimate">Estimate the best k-mer size</option>
                     </param>
                     <when value="provide">
-                        <param name="input_kmer_size" type="integer" min="1" max="50" value="" label="K-mer size" help="Haploid genome size or diploid genome size, depending on what we evaluate. In bp. Only required if the k-mer size is not provided." />
+                        <param name="input_kmer_size" type="integer" min="1" max="50" value="" label="K-mer size" help="For a human genome, the best k-mer size is k=21 for both haploid (3.1G) or diploid (6.2G).." />
                     </when>
                     <when value="estimate">
                         <param name="genome_size" type="integer" min="1000" max="70000000000" value="1000" label="Genome size" help="Haploid genome size or diploid genome size, depending on what we evaluate. In bp. Only required if the k-mer size is not provided." />
@@ -203,7 +213,7 @@
                         <option value="estimate">Estimate the best k-mer size</option>
                     </param>
                     <when value="provide">
-                        <param name="input_kmer_size" type="integer" min="1" max="50" value="" label="K-mer size" help="Haploid genome size or diploid genome size, depending on what we evaluate. In bp. Only required if the k-mer size is not provided." />
+                        <param name="input_kmer_size" type="integer" min="1" max="50" value="" label="K-mer size" help="For a human genome, the best k-mer size is k=21 for both haploid (3.1G) or diploid (6.2G)." />
                     </when>
                     <when value="estimate">
                         <param name="genome_size" type="integer" min="1000" max="70000000000" value="1000" label="Genome size" help="Haploid genome size or diploid genome size, depending on what we evaluate. In bp. Only required if the k-mer size is not provided." />