view read_quality.xml @ 58:1a052c827e88 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rseqc commit 57f71aa633a43ab02bbf05acd0c6d7f406e01f1e"
author iuc
date Thu, 28 Nov 2019 15:56:37 -0500
parents 34e4c586e3c0
children dbedfc5f5a3c
line wrap: on
line source

<tool id="rseqc_read_quality" name="Read Quality" version="@WRAPPER_VERSION@">
    <description>determines Phred quality score</description>

    <macros>
        <import>rseqc_macros.xml</import>
    </macros>

    <expand macro="requirements">
        <!--
            Required due to conda solver bug: https://github.com/conda/conda/issues/6269
            See: https://github.com/galaxyproject/tools-iuc/pull/1578 for more info
        -->
        <requirement type="package" version="3.4.1">r-base</requirement>
    </expand>

    <expand macro="stdio" />

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

    <command><![CDATA[
        read_quality.py
            --input-file '${input}'
            --out-prefix output
            -r ${reduce}
            --mapq ${mapq}
        ]]>
    </command>

    <inputs>
        <expand macro="bam_sam_param" />
        <param name="reduce" type="integer" label="Ignore Phred scores less than this amount (only applies to 'boxplot', default=1000)" value="1000" help="(--reduce)"/>
        <expand macro="mapq_param" />
        <expand macro="rscript_output_param" />
    </inputs>

    <outputs>
        <data format="pdf" name="outputheatpdf" from_work_dir="output.qual.heatmap.pdf" label="${tool.name} on ${on_string} (Heatmap pdf)" />
        <data format="pdf" name="outputboxpdf" from_work_dir="output.qual.boxplot.pdf" label="${tool.name} on ${on_string} (Boxplot pdf)" />
        <expand macro="rscript_output_data" filename="output.qual.r" />
    </outputs>

    <tests>
        <test>
            <param name="input" value="pairend_strandspecific_51mer_hg19_random.bam"/>
            <param name="rscript_output" value="true" />
            <output name="outputr" file="output.qual.r"/>
            <output name="outputheatpdf" file="output.qual.heatmap.pdf" compare="sim_size" />
            <output name="outputboxpdf" file="output.qual.boxplot.pdf" compare="sim_size" />
        </test>
    </tests>

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

According to SAM specification, if Q is the character to represent "base calling quality"
in SAM file, then Phred Quality Score = ord(Q) - 33. Here ord() is python function that
returns an integer representing the Unicode code point of the character when the argument
is a unicode object, for example, ord('a') returns 97. Phred quality score is widely used
to measure "reliability" of base-calling, for example, phred quality score of 20 means
there is 1/100 chance that the base-calling is wrong, phred quality score of 30 means there
is 1/1000 chance that the base-calling is wrong. In general: Phred quality score = -10xlog(10)P,
here P is probability that base-calling is wrong.

Inputs
++++++++++++++

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

Ignore phred scores less than this number (default=1000)
    To avoid making huge vector in R, nucleotide with certain phred score represented less than this number will be ignored. Increase this number save more memory while reduce precision. This option only applies to the 'boxplot'.

Output
++++++++++++++

1. output.qual.r
2. output.qual.boxplot.pdf
    .. image:: $PATH_TO_IMAGES/36mer.qual.plot.png
        :height: 600 px
        :width: 600 px
        :scale: 80 %
3. output.qual.heatmap.pdf
    .. image:: $PATH_TO_IMAGES/36mer.qual.heatmap.png
        :height: 600 px
        :width: 600 px
        :scale: 80 %

Heatmap: use different color to represent nucleotide density ("blue"=low density,"orange"=median density,"red"=high density")

@ABOUT@

]]>
    </help>

    <expand macro="citations" />

</tool>