<tool id="chopper" name="Chopper" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="23.2">
    <description>Filtering and trimming of long reads.</description>
    <macros>
        <token name="@TOOL_VERSION@">0.12.0</token>
        <token name="@VERSION_SUFFIX@">1</token>
    </macros>
    <requirements>
        <requirement type="package" version="@TOOL_VERSION@">chopper</requirement>
    </requirements>
    <version_command>chopper --version</version_command>
    <command detect_errors="exit_code"><![CDATA[
 chopper

 --input $input

 #if $contam
    --contam $contam    
 #end if

 #if $option_params.quality
    --quality $option_params.quality.value
 #end if

 #if $option_params.maxqual
    --maxqual $option_params.maxqual.value
 #end if

 #if $option_params.minlength
    --minlength $option_params.minlength.value
 #end if

 #if $option_params.maxlength
    --maxlength $option_params.maxlength.value
 #end if

 #if $option_params.headcrop
    --headcrop $option_params.headcrop.value
 #end if

 #if $option_params.tailcrop
    --tailcrop $option_params.tailcrop.value
 #end if

 #if $option_params.mingc
    --mingc $option_params.mingc.value
 #end if

 #if $option_params.maxgc
    --maxgc $option_params.maxgc.value
 #end if

 ##output capture
 #if $output_params.inverse == "yes"
    $output_params.inverse
 #end if 


 #if str($input.ext).endswith("gz")
    | gzip > $fq_filt
 #else
    > $fq_filt
 #end if

 ]]></command>
    <inputs>
        <param argument="--input" type="data" label="FASTQ file to check" format="fastqsanger.gz,fastqsanger" />
        <param argument="--contam" type="data" format="fasta" optional="True" label="Reference FASTA" help="FASTA file with reference to check potential contaminants against."/>    

        <section name="option_params" title="Optional Parameters" expanded="True">
            <param argument="--quality" type="integer" label="Minimal quality score" value="0" min="0" max="60" help="Sets a minimum Phred average quality score."/>
            <param argument="--maxqual" type="integer" label="Maximal quality score" value="60" min="0" max="60" help="Sets a maximum Phred average quality score."/>
            <param argument="--minlength" type="integer" label="Sets a minimum read length" value="1" min="1"  help="Minimal length of read to keep."  optional="True"/>
            <param argument="--maxlength" type="integer" label="Sets a maximum read length" help="Maximal length of read to keep" optional="True"/>
            <param argument="--headcrop" type="integer" optional="True" label="Headcrop" value="0" min="0" help="Trim N nucleotides from the start of a read."/>
            <param argument="--tailcrop" type="integer" optional="True" label="Tailcrop" value="0" min="0" help="Trim N nucleotides from the end of a read."/>
            <param argument="--mingc" type="float" optional="True" label="Minimum GC content" value="0.0" min="0.0" max="1.0" help="Sets a minimum GC content for reads to keep."/>
            <param argument="--maxgc" type="float" optional="True" label="Maximum GC content" value="1.0" min="0.0" max="1.0" help="Sets a maximum GC content for reads to keep."/>
            <param argument="--trim" type="integer" label="Q-score cutoff to trim read ends" value="0" min="0" max="60" help="Takes a quality score and will trim the ends of the reads if they are below the specified cut-off (window-size = 1)."/>
        </section>
        <section name="output_params" title="Output Parameters" expanded="False">
            <param argument="--inverse" type="boolean" checked="false" truevalue="--inverse" falsevalue="" label="Output the opposite of the normal results" help="Reverse the output results (aka, the 'failed reads')"/>
        </section>
    </inputs>
    <outputs>
	    <data name="fq_filt" format_source="input" label="${tool.name} on ${on_string}" />
    </outputs>
    <tests>
        <!-- 1) -->
        <test expect_num_outputs="1">
            <param name="input" ftype="fastqsanger" value="other-test.fastq"/>
            <output name="fq_filt" ftype="fastqsanger">
                <assert_contents>
                    <has_text text="@35febf09-dcbc-424c-987e-9f3f80fe73a5"/>
                    <has_text text="@3fda06e9-62ef-4448-9993-b90124a793d5"/>
                    <has_text text="@19d9337f-4fb6-46e5-b484-14d05f562506"/>
                </assert_contents>
            </output>
        </test>
        <!-- 2) -->
        <test expect_num_outputs="1">
            <param name="input" ftype="fastqsanger" value="other-test.fastq"/>
            <param name="contam" value="random_contam.fa"/>
            <output name="fq_filt" ftype="fastqsanger">
                <assert_contents>
                    <has_text text="@35febf09-dcbc-424c-987e-9f3f80fe73a5"/>
                    <has_text text="@3fda06e9-62ef-4448-9993-b90124a793d5"/>
                    <has_text text="@19d9337f-4fb6-46e5-b484-14d05f562506"/>
                </assert_contents>
            </output>
        </test>
        <!-- 3) -->
        <test expect_num_outputs="1">
            <param name="input" ftype="fastqsanger" value="testGC.fastq"/>
            <section name="option_params">
                <param name="mingc" value="0.3"/>
                <param name="maxgc" value="0.8"/>
            </section>
            <output name="fq_filt" ftype="fastqsanger"> 
                <assert_contents>
                    <not_has_text text='@GC20'/>
                    <not_has_text text='@GC0'/>
                    <not_has_text text='@GC100'/>
                    <has_text text='@GC50'/>
                    <has_text text='@GC80'/>
                </assert_contents>
            </output>
        </test>
        <!-- 4) fastqsanger -->
        <test expect_num_outputs="1">
            <param name="input" ftype="fastqsanger" value="wrapping_as_sanger.fastqsanger"/>
            <section name="option_params">
                <param name="maxlength" value="136"/>
            </section>
            <output name="fq_filt" ftype="fastqsanger"> 
                <assert_contents>
                    <not_has_text text='@SRR014849.203935'/>
                    <has_text text='@SRR014849.50939'/>
                    <has_text text='@SRR014849.110027'/>
                </assert_contents>
            </output>
        </test>
        <!-- 5) fastqsanger.gz -->
        <test expect_num_outputs="1">
            <param name="input" ftype="fastqsanger.gz" value="wrapping_as_sanger.fastqsanger.gz"/>
            <section name="option_params">
                <param name="maxlength" value="136"/>
            </section>
            <output name="fq_filt" ftype="fastqsanger.gz" decompress="true">
                <assert_contents>
                    <not_has_text text='@SRR014849.203935'/>
                    <has_text text='@SRR014849.50939'/>
                    <has_text text='@SRR014849.110027'/>
                </assert_contents>
            </output>
            <assert_stderr>
                <has_text text="Kept 2 reads out of 3 reads"/>
            </assert_stderr>
        </test>
    </tests>
    <help><![CDATA[
**Chopper**

Rust implementation of NanoFilt+NanoLyse, both originally written in Python. This tool, intended for long read sequencing such as PacBio or ONT, filters and trims a fastq file.
Filtering is done on average read quality and minimal or maximal read length, and applying a headcrop (start of read) and tailcrop (end of read) while printing the reads passing the filter.

Compared to the Python implementation the scope is to deliver the same results, almost the same functionality, at much faster execution times. At the moment this tool does not support filtering using a sequencing_summary file. 

**More Information** 

- **Official Repository**: `Chopper on GitHub`_

.. _Chopper on GitHub: https://github.com/wdecoster/chopper

    ]]></help>
    <citations>
        <citation type="doi">10.1093/bioinformatics/btad311</citation>
    </citations>
</tool>