<?xml version='1.0' encoding='utf-8'?>
<tool profile="16.04" id="kissplice" name="KisSplice" version="2.4.0">
    <description>is a local transcriptome assembler for SNPs, indels and AS events</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <requirements>
        <requirement type="package" version="2.4.0p1">kissplice</requirement>
    </requirements>
    <version_command>kissplice --version</version_command>
    <command><![CDATA[

        ## prepare symlinks
        #set $samples = ""

        #for x in $list_reads
            #set $filename = os.path.basename(str($x)) + "." + $x.ext
            ln -sf '${x}' '${filename}' &&
            #set $samples += " -r '" + $filename + "'"
        #end for

        kissplice
        $samples
        -k ${k}
        -c ${c}
        -C ${C}
        ${output_context}
        --counts ${counts}

    ]]></command>
    <inputs>

        <param name='list_reads' argument="-r" multiple="true" format="fasta,fastq" type="data" label="Reads file" />
        <param argument="-k" type="integer" value="41" label="K-mer length" />
        <param argument="-c" type="integer" value="2" label="K-mers present strictly less than this number of times in the dataset will be discarded (integer value required)" />
        <param argument="-C" type="float" value="0.05" label="Edges with relative coverage below this number are removed (value in [0,1) required)" />
        <param name="output_context" argument="--output-context" type="boolean" checked="false" truevalue="--output-context" falsevalue="" label="Output the maximum non-ambiguous context of each event" />
        <param argument="--counts" type="select" label="Selecting how the counts will be reported">
            <option value="0">Total counts</option>
            <option value="1">Counts on junctions</option>
            <option value="2">All counts</option>
        </param>
        <param name="output_option" type="select" label="Output type">
            <option value="alternative_file">Generate only the Alternative Splicing file</option>
            <option value="all_files">Generate all files</option>
        </param>
    </inputs>
    <outputs>
        <data name="type0a" format="fasta" label="Single SNPs or sequencing errors" from_work_dir="results/*_coherents_type_0a.fa">
            <filter>output_option == 'all_files'</filter>
        </data>
        <data name="type0b" format="fasta" label="Multiple SNPs or sequencing errors" from_work_dir="results/*_coherents_type_0b.fa">
            <filter>output_option == 'all_files'</filter>
        </data>
        <data name="type1" format="fasta" label="Alternative splicing events" from_work_dir="results/*_coherents_type_1.fa"/>
        <data name="type2" format="fasta" label="Inexact tandem repeats" from_work_dir="results/*_coherents_type_2.fa">
            <filter>output_option == 'all_files'</filter>
        </data>
        <data name="type3" format="fasta" label="Short indels" from_work_dir="results/*_coherents_type_3.fa">
            <filter>output_option == 'all_files'</filter>
        </data>
        <data name="type4" format="fasta" label="Other events" from_work_dir="results/*_coherents_type_4.fa">
            <filter>output_option == 'all_files'</filter>
        </data>
    </outputs>

    <tests>
        <test>
            <param name="list_reads" value="kissplice/reads1,kissplice/reads2" ftype="fasta" />
            <param name="k" value="25"/>
            <param name="output_option" value="alternative_file" />
            <output name="type1" file="kissplice/alternative_splicing.fa"/>
        </test>
    </tests>

    <help><![CDATA[

**Description**

KisSplice is a piece of software that enables the **analysis of RNA-seq data with or without a reference genome**. It is an exact local transcriptome assembler that allows one to identify **SNPs, indels and alternative splicing events**. It can deal with an arbitrary number of biological conditions, and will quantify each variant in each condition. It has been tested on Illumina datasets of up to 1G reads. Its memory consumption is around 5Gb for 100M reads.

KisSplice is not a full-length transcriptome assembler. This means that it will output the variable regions of the transcripts, not reconstruct them entirely. However, KisSplice results can be further aligned to a reference transcriptome (if available), or to the output of a full-length transcriptome assembler like Trinity or Oases.

------

**Website**

References, documentation and FAQ here_.

.. _here: http://kissplice.prabi.fr/


    ]]></help>
    <expand macro="citations">
        <citation type="bibtex">
         @article{PMID:22537044,
          author       = {Sacomoto, G. A. and Kielbassa, J. and Chikhi, R. and Uricaru, R. and Antoniou, P. and Sagot, M. F. and Peterlongo, P. and Lacroix, V.},
          title        = {KISSPLICE: de-novo calling alternative splicing events from RNA-seq data.},
          journal      = {BMC Bioinformatics},
          year         = {2012},
          volume       = {13 Suppl 6},
          pages        = {S5},
          url          = {http://www.ncbi.nlm.nih.gov/pubmed/22537044},
        }
        </citation>
    </expand>

</tool>