Mercurial > repos > iuc > stacks_sstacks
view stacks_sstacks.xml @ 1:c4ab391b5274 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/stacks commit 64069b03d33644632275cb79a017f5425323c75e
author | iuc |
---|---|
date | Wed, 15 Jun 2016 05:39:55 -0400 |
parents | 17469af4e858 |
children | 344cfdb7dd48 |
line wrap: on
line source
<tool id="stacks_sstacks" name="Stacks: sstacks" version="@WRAPPER_VERSION@.0"> <description>match stacks to a catalog</description> <macros> <import>macros.xml</import> </macros> <expand macro="requirements"/> <expand macro="stdio"/> <command><![CDATA[ #import re mkdir stacks_inputs stacks_outputs && #set $catalog = "" #for $input_file in $input_cat: #set $filename = str($input_file.name) #if not filename.endswith('.tsv'): #set $filename = ".tsv" #end if #if re.search('catalog\.[a-z]+(\.tsv)?$', $filename): ln -s "${input_file}" "stacks_inputs/$filename" && #if $filename.endswith('.tags.tsv'): #set catalog += " -c \"stacks_inputs/"+$filename[:-17] + "\"" #end if #end if #end for #set $samples = "" #for $input_file in $input_tags: #set $filename = str($input_file.name) #if not filename.endswith('.tsv'): #set $filename = ".tsv" #end if #if not re.search('catalog\.[a-z]+(\.tsv)?$', $filename): ln -s "${input_file}" "stacks_inputs/$filename" && #if $filename.endswith('.tags.tsv'): #set samples += " -s \"stacks_inputs/"+$filename[:-9] + "\"" #end if #end if #end for sstacks -p \${GALAXY_SLOTS:-1} $catalog $samples ## Batch description -b 1 $g $check_haplo $gapped -o stacks_outputs ]]></command> <inputs> <param name="input_cat" format="tabular,txt" type="data_collection" collection_type="list" label="Catalog files" help="output from a previous Stacks pipeline steps e.g. denovo_map, refmap or cstacks" /> <param name="input_tags" format="tabular,txt" type="data_collection" collection_type="list" label="Samples stacks" help="output from previous Stacks pipeline steps e.g. denovo_map, refmap or ustacks/pstacks" /> <param name="g" argument="-g" type="boolean" checked="false" truevalue="-g" falsevalue="" label="Base catalog matching on genomic location, not sequence identity" /> <param name="check_haplo" argument="-x" type="boolean" checked="false" truevalue="-x" falsevalue="" label="Don't verify haplotype of matching locus" /> <param name="gapped" argument="--gapped" type="boolean" checked="false" truevalue="--gapped" falsevalue="" label="Perform gapped alignments between stacks" /> </inputs> <outputs> <collection name="matches" type="list" label="Matches to the catalog on ${on_string}"> <discover_datasets pattern="(?P<name>.+\.matches)\.tsv" ext="tabular" directory="stacks_outputs" /> </collection> </outputs> <tests> <test> <param name="input_cat"> <collection type="list"> <element name="batch_1.catalog.alleles.tsv" ftype="tabular" value="genotypes/batch_1.catalog.alleles.tsv" /> <element name="batch_1.catalog.snps.tsv" ftype="tabular" value="genotypes/batch_1.catalog.snps.tsv" /> <element name="batch_1.catalog.tags.tsv" ftype="tabular" value="genotypes/batch_1.catalog.tags.tsv" /> </collection> </param> <param name="input_tags"> <collection type="list"> <element name="PopA_01.alleles.tsv" ftype="tabular" value="genotypes/PopA_01.alleles.tsv" /> <element name="PopA_01.matches.tsv" ftype="tabular" value="genotypes/PopA_01.matches.tsv" /> <element name="PopA_01.snps.tsv" ftype="tabular" value="genotypes/PopA_01.snps.tsv" /> <element name="PopA_01.tags.tsv" ftype="tabular" value="genotypes/PopA_01.tags.tsv" /> <element name="PopA_02.alleles.tsv" ftype="tabular" value="genotypes/PopA_02.alleles.tsv" /> <element name="PopA_02.matches.tsv" ftype="tabular" value="genotypes/PopA_02.matches.tsv" /> <element name="PopA_02.snps.tsv" ftype="tabular" value="genotypes/PopA_02.snps.tsv" /> <element name="PopA_02.tags.tsv" ftype="tabular" value="genotypes/PopA_02.tags.tsv" /> </collection> </param> <output_collection name="matches"> <element name="PopA_01.matches"> <assert_contents> <has_text text="generated on " /> </assert_contents> </element> </output_collection> </test> <test> <param name="input_cat"> <collection type="list"> <element name="batch_1.catalog.alleles.tsv" ftype="tabular" value="genotypes/batch_1.catalog.alleles.tsv" /> <element name="batch_1.catalog.snps.tsv" ftype="tabular" value="genotypes/batch_1.catalog.snps.tsv" /> <element name="batch_1.catalog.tags.tsv" ftype="tabular" value="genotypes/batch_1.catalog.tags.tsv" /> </collection> </param> <param name="input_tags"> <collection type="list"> <element name="PopA_01.alleles.tsv" ftype="tabular" value="genotypes/PopA_01.alleles.tsv" /> <element name="PopA_01.matches.tsv" ftype="tabular" value="genotypes/PopA_01.matches.tsv" /> <element name="PopA_01.snps.tsv" ftype="tabular" value="genotypes/PopA_01.snps.tsv" /> <element name="PopA_01.tags.tsv" ftype="tabular" value="genotypes/PopA_01.tags.tsv" /> <element name="PopA_02.alleles.tsv" ftype="tabular" value="genotypes/PopA_02.alleles.tsv" /> <element name="PopA_02.matches.tsv" ftype="tabular" value="genotypes/PopA_02.matches.tsv" /> <element name="PopA_02.snps.tsv" ftype="tabular" value="genotypes/PopA_02.snps.tsv" /> <element name="PopA_02.tags.tsv" ftype="tabular" value="genotypes/PopA_02.tags.tsv" /> </collection> </param> <param name="gapped" value="true" /> <param name="check_haplo" value="true" /> <output_collection name="matches"> <element name="PopA_01.matches"> <assert_contents> <has_text text="generated on " /> </assert_contents> </element> </output_collection> </test> </tests> <help> <![CDATA[ .. class:: infomark **What it does** Sets of stacks constructed by the ustacks or pstacks programs can be searched against a catalog produced by cstacks. In the case of a genetic map, stacks from the progeny would be matched against the catalog to determine which progeny contain which parental alleles. -------- **Input files** Output from denovo_map, refmap or cstacks/ustacks/pstacks **Output files** - XXX.tags.tsv file: See `Stacks output description <http://catchenlab.life.illinois.edu/stacks/manual/#files>`_ Notes: For the tags file, each stack will start in the file with a consensus sequence for the entire stack followed by the flags for that stack. Then, each individual read that was merged into that stack will follow. The next stack will start with another consensus sequence. - XXX.snps.tsv file: See `Stacks output description <http://catchenlab.life.illinois.edu/stacks/manual/#files>`_ Notes: If a stack has two SNPs called within it, then there will be two lines in this file listing each one. - XXX.alleles.tsv file: See `Stacks output description <http://catchenlab.life.illinois.edu/stacks/manual/#files>`_ - XXX.matches.tsv file: See `Stacks output description <http://catchenlab.life.illinois.edu/stacks/manual/#files>`_ Notes: Each line in this file records a match between a catalog locus and a locus in an individual, for a particular haplotype. The Batch ID plus the Catalog ID together represent a unique locus in the entire population, while the Sample ID and the Stack ID together represent a unique locus in an individual sample. @STACKS_INFOS@ ]]> </help> <expand macro="citation" /> </tool>