<tool id="rseqc_deletion_profile" name="Deletion Profile" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@GALAXY_VERSION@">
    <description>
     calculates the distributions of deleted nucleotides across reads
    </description>
    <macros>
        <import>rseqc_macros.xml</import>
    </macros>

    <expand macro="bio_tools"/>

    <expand macro="requirements" />

    <expand macro="stdio" />

    <version_command><![CDATA[deletion_profile.py --version]]></version_command>

    <command><![CDATA[
        @BAM_SAM_INPUTS@
        deletion_profile.py -i 'input.${extension}' -o output -l ${read_align_length} -n ${read_num} -q ${mapq}
        ]]>
    </command>

    <inputs>
        <expand macro="bam_param" />
        <expand macro="readlength_param" />
        <expand macro="readnum_param" />
        <expand macro="mapq_param" />
        <expand macro="rscript_output_param" />
    </inputs>

    <outputs>
        <expand macro="pdf_output_data" filename="output.deletion_profile.pdf" />
        <expand macro="xls_output_data" filename="output.deletion_profile.txt" />
        <expand macro="rscript_output_data" filename="output.deletion_profile.r" />
    </outputs>

    <tests>
        <test>
            <param name="input" value="pairend_strandspecific_51mer_hg19_chr1_1-100000.bam" />
            <param name="read_align_length" value="101" />
            <param name="rscript_output" value="true" />
            <output name="outputpdf" file="output.deletion_profile.pdf" compare="sim_size" />
            <output name="outputxls" file="output.deletion_profile.txt" ftype="tabular"/>
            <output name="outputr" file="output.deletion_profile_r" />
        </test>
    </tests>

    <help><![CDATA[
deletion_profile.py
+++++++++++++++++++

Calculate the distributions of deleted nucleotides across reads.

Inputs
++++++

Input BAM/SAM file
    Alignment file in BAM/SAM format.

Alignment length of read
    It is usually set to the orignial read length. For example, all these cigar
    strings ("101M", "68M140N33M", "53M1D48M") suggest the read alignment
    length is 101. [required]

Number of aligned reads used
    Number of aligned reads with deletions used to calculate the deletion
    profile. default=1000000

Minimum mapping quality
    Minimum mapping quality for an alignment to be considered as "uniquely
    mapped". default=30

Sample Output
++++++++++++++

.. image:: $PATH_TO_IMAGES/out.deletion_profile.png
   :height: 600 px
   :width: 600 px
   :scale: 80 %

@ABOUT@

]]>

    </help>

    <expand macro="citations" />

</tool>