<tool id="nanopore_qc" name="NanoporeQC" version="0.002">

    <description>Quality report for nanopore data</description>

    <requirements>
        <requirement type="package" version="3.3.2">r-base</requirement>
        <requirement type="package" version="2.2.1">r-ggplot2</requirement>
        <requirement type="package" version="1.8.4">r-plyr</requirement>
        <requirement type="package" version="1.4.2">r-reshape2</requirement>
        <requirement type="package" version="0.2.2">r-readr</requirement>
        <requirement type="package" version="2.1.13">r-yaml</requirement>
        <requirement type="package" version="0.4.1">r-scales</requirement>
        <requirement type="package" version="1.4.1">r-futile.logger</requirement>
        <requirement type="package" version="1.10.4">r-data.table</requirement>
        <requirement type="package" version="1.3.2">r-optparse</requirement>
        <requirement type="package" version="1.8_16">r-mgcv</requirement>
        <requirement type="package" version="0.66">perl-yaml-libyaml</requirement>
    </requirements>

    <stdio>
        <exit_code range="1:" level="fatal" description="Unknown error" />
        <exit_code range=":-1" level="fatal" description="Unknown error" />
    </stdio>

    <version_command>
    <![CDATA[
        echo '0.001'
    ]]>
    </version_command>

    <command>
    <![CDATA[

        mkdir -p '${html_file.files_path}'

        &&

        Rscript '${__tool_directory__}/nanopore_qc.R'
            -i '$input'
            -o '${html_file.files_path}'
            -q '$q_cutoff'
            $discard_failed

        &&

        perl '${__tool_directory__}/yaml_to_html.pl'

            '${html_file.files_path}/summary.yaml'
            '${html_file.files_path}'
            '$html_file'

    ]]>
    </command>

    <inputs>
        <param name="input" type="data" format="tabular" label="Basecall summary file from Albacore"/>
        <param name="q_cutoff" type="float" value="10" min="0" label="Quality cutoff for QC calculations" />
        <param name="discard_failed" type="boolean" checked="false" truevalue="--discard_failed T" falsevalue="" label="Ignore reads failing quality filter" />
    </inputs>
    <outputs>
        <data format="html" name="html_file" label="NanoporeQC on ${on_string}" />
    </outputs>
    <tests>
        <!-- binary image blobs can differ, which is why we use line matching
        instead of file matching -->

        <!-- quality cutoff 10 -->
        <test>
            <param name="input" value="sequencing_summary.txt.small" ftype="tabular" />
            <param name="q_cutoff" value="10" />
            <output name="html_file">
                <assert_contents>
                    <has_text_matching expression="&lt;td&gt;Total Reads&lt;/td&gt;\s*&lt;td&gt;9990&lt;/td&gt;" />
                    <has_text_matching expression="&lt;td&gt;Total Reads&lt;/td&gt;\s*&lt;td&gt;7952&lt;/td&gt;" />
                </assert_contents>
            </output>
        </test>

        <!-- quality cutoff 6 -->
        <test>
            <param name="input" value="sequencing_summary.txt.small" ftype="tabular" />
            <param name="q_cutoff" value="6" />
            <output name="html_file">
                <assert_contents>
                    <has_text_matching expression="&lt;td&gt;Total Reads&lt;/td&gt;\s*&lt;td&gt;9990&lt;/td&gt;" />
                    <has_text_matching expression="&lt;td&gt;Total Reads&lt;/td&gt;\s*&lt;td&gt;9254&lt;/td&gt;" />
                </assert_contents>
            </output>
        </test>

        <!-- quality cutoff 6, discard filtered -->
        <test>
            <param name="input" value="sequencing_summary.txt.small" ftype="tabular" />
            <param name="q_cutoff" value="6" />
            <param name="discard_failed" value="true" />
            <output name="html_file">
                <assert_contents>
                    <has_text_matching expression="&lt;td&gt;Total Reads&lt;/td&gt;\s*&lt;td&gt;9107&lt;/td&gt;" />
                    <has_text_matching expression="&lt;td&gt;Total Reads&lt;/td&gt;\s*&lt;td&gt;9107&lt;/td&gt;" />
                </assert_contents>
            </output>
        </test>

    </tests>
    <help>
<![CDATA[
    `NanoporeQC` produces a QC report for nanopore sequencing data. It utilizes
    the summary table produced by Albacore, which saves time by avoiding the
    need to re-parse the FAST5 read data. `NanoporeQC` was originally forked
    from the `minion_qc` codebase (https://github.com/roblanf/minion_qc).

    The only tuneable parameter is the quality value cutoff. Two summaries are
    produced for each table and plot -- one for all reads (with length >= 0)
    and one for reads with a quality value greater than the cutoff.
]]>
    </help>
    <citations>
    </citations>
</tool>