view make.contigs.xml @ 0:d9268e6b2ea0 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mothur commit 180a403421967d36f995941b1a4561349d75cfc5
author iuc
date Fri, 24 Jun 2016 16:37:19 -0400
parents
children 2de7963eb3ba
line wrap: on
line source

<tool profile="16.07" id="mothur_make_contigs" name="Make.contigs" version="@WRAPPER_VERSION@.0">
    <description>Aligns paired forward and reverse fastq files to contigs as fasta and quality</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro="requirements"/>
    <expand macro="version_command"/>
    <command detect_errors="aggressive"><![CDATA[
        ## create symlinks to input datasets
        ln -s "$ffastq" ffastq.dat &&
        ln -s "$rfastq" rfastq.dat &&
        #if $oligo.add == "yes":
            ln -s "$oligo.oligos" oligo.oligos.dat &&
            ln -s "$oligo.findex" oligo.findex.dat &&
            ln -s "$oligo.rindex" oligo.rindex.dat &&
        #end if

        echo 'make.contigs(
            ffastq=ffastq.dat,
            rfastq=rfastq.dat,
            align=$align,
            #if $oligo.add == "yes":
                oligos=oligo.oligos.dat,
                bdiffs=$oligo.bdiffs,
                pdiffs=$oligo.pdiffs,
                tdiffs=$oligo.tdiffs,
                #if $oligo.findex:
                    findex=oligo.findex.dat,
                #end if
                #if $oligo.rindex:
                    rindex=oligo.rindex.dat,
                #end if
            #end if
            match=$match,
            mismatch=$mismatch,
            gapopen=$gapopen,
            gapextend=$gapextend,
            processors='\${GALAXY_SLOTS:-8}'
        )'
        | sed 's/ //g'  ## mothur trips over whitespace
        | mothur
    ]]></command>
    <inputs>
        <param name="ffastq" type="data" format="fastq" label="ffastq - Forward Fastq Sequence file"/>
        <param name="rfastq" type="data" format="fastq" label="rfastq - Reverse Fastq Sequence file"/>
        <param name="align" type="select" label="align - Select a pairwise alignment method" help="">
            <option value="needleman" selected="true">needleman (default) </option>
            <option value="gotoh">gotoh</option>
            <option value="kmer">kmer</option>
        </param>
        <conditional name="oligo">
            <param name="add" type="select" label="Trim with an oligos file?" help="">
                <option value="no">no</option>
                <option value="yes">yes</option>
            </param>
            <when value="no"/>
            <when value="yes">
                <param name="oligos" type="data" format="mothur.oligos" optional="true" label="oligos - barcodes and primers" help="a file that can contain the sequences of the forward and reverse primers and barcodes and their sample identifier. Each line of the oligos file can start with the key words &quot;forward&quot;, &quot;reverse&quot; and &quot;barcode&quot; or it can start with a &quot;#&quot; to tell mothur to ignore that line of the oligos file"/>
                <param name="bdiffs" type="integer" value="0" min="0" label="bdiffs - number of differences to allow in the barcode (default 0)"/>
                <param name="pdiffs" type="integer" value="0" min="0" label="pdiffs - number of differences to allow in the primer (default 0)"/>
                <param name="tdiffs" type="integer" value="0" min="0" label="tdiffs - total number of differences to allow in primer and barcode (default 0)"/>
                <param name="findex" type="data" format="fastq" optional="true" label="forward index file (fastq)"/>
                <param name="rindex" type="data" format="fastq" optional="true" label="revese index file (fastq)"/>
            </when>
        </conditional>
        <param name="match" type="integer" value="1" label="match - Pairwise alignment reward for a match"/>
        <param name="mismatch" type="integer" value="-1" label="mismatch - Pairwise alignment penalty for a mismatch"/>
        <param name="gapopen" type="integer" value="-2" label="gapopen - Pairwise alignment penalty for opening a gap"/>
        <param name="gapextend" type="integer" value="-1" label="gapextend - Pairwise alignment penalty for extending a gap"/>
    </inputs>
    <outputs>
        <expand macro="logfile-output"/>
        <data name="fasta" format="fasta" from_work_dir="ffastq*.trim.*.fasta" label="${tool.name} on ${on_string}: trim.contigs.fasta"/>
        <data name="qual" format="qual" from_work_dir="ffastq*.trim.*.qual" label="${tool.name} on ${on_string}: trim.contigs.qual"/>
        <data name="scrapfasta" format="fasta" from_work_dir="ffastq*.scrap.*.fasta" label="${tool.name} on ${on_string}: scrap.contigs.fasta"/>
        <data name="scrapqual" format="qual" from_work_dir="ffastq*.scrap.*.qual" label="${tool.name} on ${on_string}: scrap.contigs.qual"/>
        <data name="report" format="txt" from_work_dir="ffastq*.contigs.report" label="${tool.name} on ${on_string}: report"/>
    </outputs>
    <tests>
        <test>
            <param name="ffastq" value="Mock_S280_L001_R1_001_small.fastq"/>
            <param name="rfastq" value="Mock_S280_L001_R2_001_small.fastq"/>
            <output name="fasta" file="Mock_S280_L001_R1_001_small.trim.contigs.fasta" ftype="fasta"/>
            <output name="qual" file="Mock_S280_L001_R1_001_small.trim.contigs.qual" ftype="qual"/>
            <output name="report" file="Mock_S280_L001_R1_001_small.contigs.report" ftype="txt"/>
            <expand macro="logfile-test"/>
        </test>
        <test>
            <param name="ffastq" value="Mock_S280_L001_R1_001_small.fastq"/>
            <param name="rfastq" value="Mock_S280_L001_R2_001_small.fastq"/>
            <param name="align" value="gotoh"/>
            <param name="match" value="2"/>
            <param name="mismatch" value="-2"/>
            <param name="gapopen" value="-3"/>
            <param name="gapextend" value="-2"/>
            <output name="fasta" md5="48e32c65bd9f064c5c0b4ea7695cabe9" ftype="fasta"/>
            <output name="qual" md5="1e7778cee0d86bfa2759a07bb4356165" ftype="qual"/>
            <output name="report" md5="5274725ef45890fd6da4650d5d536173" ftype="txt"/>
            <expand macro="logfile-test"/>
        </test>
    </tests>
    <help>
<![CDATA[

@MOTHUR_OVERVIEW@

**Command Documenation**

The make.contigs_ command reads a forward fastq file and a reverse fastq file and outputs new fasta and quality files.

.. _make.contigs: http://www.mothur.org/wiki/Make.contigs

v.1.27.0: Updated to use Mothur 1.33. Added findex and rindex parmaeters, optionally used with the oligos file.
]]>
    </help>
    <expand macro="citations"/>
</tool>