Mercurial > repos > iuc > meryl_trio_mode
view trio-mode.xml @ 0:c85aec7804a4 draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meryl commit a95f5b04b1219489a327a622184633a561fe5ac0
author | iuc |
---|---|
date | Thu, 25 Apr 2024 21:01:24 +0000 |
parents | |
children |
line wrap: on
line source
<tool id="meryl_trio_mode" name="Meryl" version="@TOOL_VERSION@+@GALAXY_TOOL_VERSION@@SUFFIX_VERSION@" profile="@PROFILE@"> <description>build hap-mers databases for trios</description> <macros> <import>macros.xml</import> </macros> <expand macro="edam_ontology"/> <xrefs> <xref type="bio.tools">meryl</xref> </xrefs> <expand macro="requirements"/> <version_command>meryl --version</version_command> <command detect_errors="exit_code"><![CDATA[ export GALAXY_MEMORY_GB=\$((\${GALAXY_MEMORY_MB:-8192}/1024)) && export MERQURY=\$(dirname \$(command -v merqury.sh))/../share/merqury/ && #if $options_kmer_size.kmer_size == 'estimate' #from math import log #set size=int(log(int($options_kmer_size.genome_size)*(1-float($options_kmer_size.collision_rate))/float($options_kmer_size.collision_rate))/log(4)) #elif $options_kmer_size.kmer_size == 'provide' #set size=$options_kmer_size.input_kmer_size #end if #for $i, $read in enumerate($paternal_reads): mkdir 'paternal${i}.meryl' && ln -s '$read' ./input_paternal_${i}.${read.ext} && meryl count k=$size memory=\$GALAXY_MEMORY_GB threads=\${GALAXY_SLOTS:-1} ./input_paternal_${i}.${read.ext} 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' && ln -s '$read' ./input_maternal_${i}.${read.ext} && meryl count k=$size memory=\$GALAXY_MEMORY_GB threads=\${GALAXY_SLOTS:-1} ./input_maternal_${i}.${read.ext} 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' && ln -s '$read' ./input_child_${i}.${read.ext} && meryl count k=$size memory=\$GALAXY_MEMORY_GB threads=\${GALAXY_SLOTS:-1} ./input_child_${i}.${read.ext} output 'child${i}.meryl' && #end for meryl union-sum child*.meryl output child.meryl && meryl histogram child.meryl > read-db.hist && tar -czf 'read-db.meryldb' child.meryl && ## mat specific kmers meryl difference mat.meryl pat.meryl output mat.only.meryl && meryl histogram mat.only.meryl > mat.only.hist && java -jar -Xmx1g \$MERQURY/eval/kmerHistToPloidyDepth.jar mat.only.hist > mat.only.ploidy && VAR=`sed -n 2p mat.only.ploidy | awk '{print \$NF}'` && meryl greater-than \$VAR output mat.only.filt.meryl mat.only.meryl && ## pat specific kmers meryl difference pat.meryl mat.meryl output pat.only.meryl && meryl histogram pat.only.meryl > pat.only.hist && java -jar -Xmx1g \$MERQURY/eval/kmerHistToPloidyDepth.jar pat.only.hist > pat.only.ploidy && VAR=`sed -n 2p pat.only.ploidy | awk '{print \$NF}'` && meryl greater-than \$VAR output pat.only.filt.meryl pat.only.meryl && ## shared kmers meryl intersect output shared.meryl mat.meryl pat.meryl && ## mat hapmers meryl intersect output mat.inherited.meryl child.meryl mat.only.filt.meryl && meryl histogram mat.inherited.meryl > mat.inherited.hist && java -jar -Xmx1g \$MERQURY/eval/kmerHistToPloidyDepth.jar mat.inherited.hist > mat.inherited.ploidy && VAR=`sed -n 2p mat.inherited.ploidy | awk '{print \$NF}'` && meryl greater-than \$VAR output mat.hapmer.meryl mat.inherited.meryl && tar -czf 'mat.meryldb' mat.hapmer.meryl && ## pat hapmers meryl intersect output pat.inherited.meryl child.meryl pat.only.filt.meryl && meryl histogram pat.inherited.meryl > pat.inherited.hist && java -jar -Xmx1g \$MERQURY/eval/kmerHistToPloidyDepth.jar pat.inherited.hist > pat.inherited.ploidy && VAR=`sed -n 2p pat.inherited.ploidy | awk '{print \$NF}'` && meryl greater-than \$VAR output pat.hapmer.meryl pat.inherited.meryl && tar -czf 'pat.meryldb' pat.hapmer.meryl && echo 'K-mer size: ${size}' ]]></command> <inputs> <param name="child_reads" type="data" format="fastq,fasta,fastq.gz,fasta.gz" multiple="true" label="F1 reads" help="Select F1 reads used for generating the assembly"/> <param name="paternal_reads" type="data" format="fastq,fasta,fastq.gz,fasta.gz" multiple="true" optional="true" label="Paternal reads" help="Select the paternal reads used for generating the assembly"/> <param name="maternal_reads" type="data" format="fastq,fasta,fastq.gz,fasta.gz" multiple="true" optional="true" label="Maternal reads" help="Select the maternal reads used for generating the assembly"/> <conditional name="options_kmer_size"> <param name="kmer_size" type="select" label="K-mer size selector"> <option value="provide">Set a k-mer size</option> <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="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."/> <param name="collision_rate" type="float" min="0.0001" max="0.01" value="0.001" label="Tolerable collision rate" help="Tolerable collision rate. By default is 0.001."/> </when> </conditional> </inputs> <outputs> <data name="read_db" format="meryldb" from_work_dir="read-db.meryldb" label="${tool.name} on ${on_string}: read-db.meryldb"/> <data name="read_db_hist" format="tabular" from_work_dir="read-db.hist" label="${tool.name} on ${on_string}: read-db histogram"/> <data name="pat_db" format="meryldb" from_work_dir="pat.meryldb" label="${tool.name} on ${on_string}: pat.meryldb"/> <data name="pat_db_hist" format="tabular" from_work_dir="pat.inherited.hist" label="${tool.name} on ${on_string}: paternal inherited histogram"/> <data name="mat_db" format="meryldb" from_work_dir="mat.meryldb" label="${tool.name} on ${on_string}: mat.meryldb"/> <data name="mat_db_hist" format="tabular" from_work_dir="mat.inherited.hist" label="${tool.name} on ${on_string}: maternal inherited histogram"/> </outputs> <tests> <test expect_num_outputs="6"> <conditional name="options_kmer_size"> <param name="kmer_size" value="provide"/> <param name="input_kmer_size" value="21"/> </conditional> <param name="maternal_reads" value="maternal.fasta"/> <param name="paternal_reads" value="paternal.fasta"/> <param name="child_reads" value="child.fasta"/> <output name="read_db" ftype="meryldb"> <assert_contents> <has_size value="3362942" delta="2000"/> <expand macro="meryldb_archive_assumptions"/> </assert_contents> </output> <output name="read_db_hist" file="output_23.read-db.hist"/> <output name="pat_db" ftype="meryldb"> <assert_contents> <has_size value="120610" delta="400"/> <expand macro="meryldb_archive_assumptions"/> </assert_contents> </output> <output name="pat_db_hist" file="output_23.pat.hist"/> <output name="mat_db" ftype="meryldb"> <assert_contents> <has_size value="67883" delta="300"/> <expand macro="meryldb_archive_assumptions"/> </assert_contents> </output> <output name="mat_db_hist" file="output_23.mat.hist"/> </test> </tests> <help> .. class:: infomark **Purpose** Meryl is the k-mer counter. This tool builds hap-mer databases for trios, in accordance with `merqury's recommended guidelines. <https://github.com/marbl/merqury/wiki/1.-Prepare-meryl-dbs#3-build-hap-mer-dbs-for-trios>`_ </help> <expand macro="citations"/> </tool>