<tool id="ribotaper_ribosome_profiling" name="ribotaper part 3: ribosome profiling" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile='20.01'>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro='bio_tools'/>
    <expand macro='requirements'/>
    <stdio>
        <exit_code range="1:" />
    </stdio>
    <command><![CDATA[
        #if $annotation_path.is_of_type("tgz")
            tar -xzvf '${annotation_path}' &&
        #else
             tar -xvf '${annotation_path}' &&
        #end if
        Ribotaper.sh
            '${ribo_bam}'
            '${rna_bam}'
            './annotation_path'
            '${read_lenghts_ribo}'
            '${cutoff}'
            \${GALAXY_SLOTS:-12}
    ]]></command>
    <inputs>
        <param name="annotation_path" type="data" format="tgz,tar" label="Annotation path" help="Please run 'ribotaper part 1' to generate the archive."/>
        <param name="ribo_bam" type="data" format="BAM" label="Ribo-seq alignment file" help="Ribo-seq alignment file in BAM format."/>
        <param name="rna_bam" type="data" format="BAM" label="RNA-seq alignment file" help="RNA-seq alignment file in BAM format."/>
        <param name="read_lenghts_ribo" type="text" value="26,28,29"  label="Read length" help="Read lengths, comma-separated values, which are used for P-site calculation. Example, 26,28,29.
            Please run 'ribotaper part 2' to deterimine appropriate values.">
            <sanitizer invalid_char="">
                <valid initial="string.digits">
                    <add value="," />
                </valid>
            </sanitizer>
            <validator type="regex">[0-9,]+</validator>
        </param>
        <param name="cutoff" type="text" value="9,12,12"  label="Cutoff" help="Cutoffs, comma-separated values, which are used for P-sites calculation. Example, 9,12,12.
            Please run 'ribotaper part 2' to deterimine appropriate values.">
            <sanitizer invalid_char="">
                <valid initial="string.digits">
                    <add value="," />
                </valid>
            </sanitizer>
            <validator type="regex">[0-9,]+</validator>
        </param>
    </inputs>
    <outputs>
        <data name="output1" format="pdf" from_work_dir="quality_check_plots.pdf" label="${tool.name} on ${on_string}: QC plots"/>
        <data name="output2" format="tabular" from_work_dir="ORFs_genes_found" label="${tool.name} on ${on_string}: Summary of translated ORFs"/>
        <data name="output3" format="tabular" from_work_dir="ORFs_max" label="${tool.name} on ${on_string}: Translated ORFs (max)"/>
        <data name="output4" format="tabular" from_work_dir="ORFs_max_filt" label="${tool.name} on ${on_string}: Translated ORFs (max_filt)"/>
        <data name="output5" format="bed" from_work_dir="translated_ORFs_sorted.bed" label="${tool.name} on ${on_string}: Translated ORFs (sorted)"/>
        <data name="output6" format="bed" from_work_dir="translated_ORFs_filtered_sorted.bed" label="${tool.name} on ${on_string}: Translated ORFs (filtered/sorted)"/>
        <data name="output7" format="fasta" from_work_dir="protein_db_max.fasta" label="${tool.name} on ${on_string}: Protein DB"/>
        <data name="output8" format="pdf" from_work_dir="Final_ORF_results.pdf" label="${tool.name} on ${on_string}: ORF categories (length/coverage)"/>
    </outputs>
    <tests>
        <test expect_num_outputs="8">
            <param name="annotation_path" value="annotation_path.tgz"/>
            <param name="ribo_bam" value="test_ribo.bam"/>
            <param name="rna_bam" value="test_rna.bam"/>
            <output name="output2" file="ORFs_genes_found"/>
        </test>
    </tests>
    <help><![CDATA[
RiboTaper is an analysis pipeline for Ribosome Profiling
(Ribo-seq) experiments,
which exploits the triplet periodicity of
ribosomal footprints to call translated regions.
See
https://ohlerlab.mdc-berlin.de/software/RiboTaper_126/ for details.


The Ribotaper Galaxy tool set consists of three tools:

  - ``ribotaper part 1``: creation of annotation files
  - ``ribotaper part 2``: metagene analysis for P-sites definition
  - ``ribotaper part 3``: ribosome profiling

The order of execution should follow:
``ribotaper part 1, part 2 and part 3``.

The current tool is ``ribotaper part 3``,
ribosome profiling.

Outputs
--------

**QC plots**:
This plot provides the user statistics about the Ribo-seq and RNA seq data used, together with the assessment of the P-sites calculations.
Important values are the pie chart showing the agreement between the frame (defined by the P-sites position) and the annotated frame. Reliable P-sites calculations produce an agreement above 90%.
Very important are also the length/coverage statistics for the Ribo-seq (bottom right):
This shows how the P-site calculations can be used to detect active translation in regions of different length and coverage, in a way the user can estimate the precision of the Ribo-seq data, and understand the level or resolution the data allows.

**Summary of translated ORFs**:
Tab-separated values for the number of ORFs found and their corresponding genes, for the different ORF categories.

**Translated ORFs (max, max_filt)**:
Tab-separated file containing information about detected ORFs.
Translated ORFs (max_filt) contains ORFs filtered for excessive multimapping and ORFs in non-coding genes overlapping known coding regions (recommended for further analysis).

**Protein DB**:
Fasta file of the detected ORFs peptide sequence, suitable as an alternative protein database (not filtered for multimapping)

**Translated ORFs  (sorted, filtered/sorted)**:
BED files with genomic coordinates for the detected ORFs. The total number of P-sites along the ORF is reported on the 5th column.

**ORF categories (length/coverage)**:
PDF file containing info about the number of ORFs found, together with their length and coverage per category/annotation.

]]></help>
    <expand macro="citations"/>
</tool>