<tool id="abricate" name="ABRicate" version="0.3">
    <requirements>
        <requirement type="package" version="0.3">abricate</requirement>
    </requirements>

    <version_command>abricate --version</version_command>

    <command detect_errors="exit_code"><![CDATA[
        abricate '$fasta_input'
        #if $settings.advanced == "advanced"
            $settings.no_header
            #if str($settings.min_dna_id)
                --minid=$settings.min_dna_id
            #end if
            #if str($settings.cull_limit)
                --cull=$settings.cull_limit
            #end if
        #end if
        > '$report'
    ]]></command>

    <inputs>
        <param name="fasta_input" type="data" format="fasta" label="FASTA file" help="To screen for antibiotic resistant genes" />
        <conditional name="settings">
            <param name="advanced" type="select" label="Specify advanced parameters">
                <option value="simple" selected="true">No, use program defaults.</option>
                <option value="advanced">Yes, see full parameter list.</option>
            </param>
            <when value="simple">
            </when>
            <when value="advanced">
                <param name="no_header" argument="--noheader" type="boolean" truevalue="--noheader" falsevalue="" label="Suppress header" help="Suppress output file's column headings" />
                <param name="min_dna_id" argument="--minid" type="float" value="75" min="0" max="100" optional="true" label="Minimum DNA identity" />
                <param name="cull_limit" argument="--cull" type="integer" value="1" min="0" optional="true" label="Culling limit (encompassing hits)" />
            </when>
        </conditional>
    </inputs>

    <outputs>
        <data name="report" format="tabular" label="${tool.name} on ${on_string} report file" />
    </outputs>

    <tests>
        <!-- Basic test - will produce no results. -->
        <test>
            <param name="fasta_input" value="Acetobacter.fna"/>
            <param name="advanced" value="simple"/>
            <output name="report" ftype="tabular" file="output_noresults.txt" />
        </test>

        <!-- Basic test - will produce results. -->
        <test>
            <param name="fasta_input" value="MRSA0252.fna"/>
            <param name="advanced" value="simple"/>
            <output name="report" ftype="tabular" file="output_mrsa.txt" compare="contains"/>
        </test>

        <!-- Advanced test - No header. -->
        <test>
            <param name="fasta_input" value="MRSA0252.fna"/>
            <param name="advanced" value="advanced"/>
            <param name="no_header" value="True"/>
            <output name="report" ftype="tabular" file="output_noheader.txt" compare="contains"/>
        </test>

        <!-- Advanced test - Cull limit 100 -->
        <test>
            <param name="fasta_input" value="MRSA0252.fna"/>
            <param name="advanced" value="advanced"/>
            <param name="cull_limit" value="100"/>
            <output name="report" ftype="tabular" file="output_cull100.txt" compare="contains"/>
        </test>

        <!-- Advanced test - Min DNA ID 100 -->
        <test>
            <param name="fasta_input" value="MRSA0252.fna"/>
            <param name="advanced" value="advanced"/>
            <param name="min_dna_id" value="100"/>
            <output name="report" ftype="tabular" file="output_minid100.txt" compare="contains"/>
        </test>
    </tests>

    <help><![CDATA[
**What it does**

Given a FASTA contig file, ABRicate will perform a mass screening of contigs and identify the presence of antibiotic resistance genes. The database of these genes is built from ResFinder.

**Output**

ABRicate will produce a tab-seperated output file with the following outputs:

+---------------+------------------------------------------------+
|  Column       | Description                                    |
+===============+================================================+
|  FILE         | The filename this hit came from                |
+---------------+------------------------------------------------+
|  SEQUENCE     | The sequence in the filename                   |
+---------------+------------------------------------------------+
|  START        | Start coordinate in the sequence               |
+---------------+------------------------------------------------+
|  END          | End coordinate                                 |
+---------------+------------------------------------------------+
|  GENE         | ABR gene                                       |
+---------------+------------------------------------------------+
|  COVERAGE     | What proportion of the gene is in our sequence |
+---------------+------------------------------------------------+
|  COVERAGE_MAP | A visual represenation                         |
+---------------+------------------------------------------------+
|  GAPS         | Was there any gaps in the                      |
|               | alignment - possible pseudogene?               |
+---------------+------------------------------------------------+
|  %COVERAGE    | Proportion of gene covered                     |
+---------------+------------------------------------------------+
|  %IDENTITY    | Proportion of exact nucleotide matches         |
+---------------+------------------------------------------------+


**Example Output**

::

	#FILE   SEQUENCE    START   END GENE    COVERAGE    COVERAGE_MAP    GAPS    %COVERAGE   %IDENTITY
	Processing: 6008.fna
	Found 12 ABR genes in 6008.fna
	Klebsiella.fna  NC_021232.1 872545  872964  fosA    1-420/420   =============== 0   100.00  100.00
	Klebsiella.fna  NC_021232.1 1381252 1382427 oqxA    1-1176/1176 =============== 0   100.00  99.32
	Klebsiella.fna  NC_021232.1 2584899 2585759 blaSHV1 1-861/861   =============== 0   100.00  99.88
    ]]></help>

    <citations>
        <citation type="bibtex">
@UNPUBLISHED{Seemann2016,
    author = {Seemann, Torsten},
    title = {ABRicate: mass screening of contigs for antiobiotic resistance genes},
    year = {2016},
    url = {https://github.com/tseemann/abricate},
}
        </citation>
    </citations>
</tool>