view diagram.xml @ 5:8f642a41cd94 draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/cnvkit commit fc1282ec68b346988203ead860e9b9d6a47e9efb
author iuc
date Sat, 01 Mar 2025 11:58:15 +0000
parents 27ecb3477d34
children
line wrap: on
line source

<tool id="cnvkit_diagram" name="CNVkit Diagram" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="21.05">
    <description>Draw copy number on chromosomes as a diagram</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro="xrefs"/>
    <expand macro="creators"/>
    <expand macro="requirements"/>
    <command detect_errors="exit_code"><![CDATA[  
        ln -s '$input_cnr_file' ./tumor.cnr &&
        #if $advanced_settings.segment
            ln -s '$advanced_settings.segment' ./segment.cns &&
        #end if
        cnvkit.py diagram
            ./tumor.cnr
            --output sample-diagram.pdf
            #if $advanced_settings.segment
                --segment ./segment.cns
            #end if
            #if str($advanced_settings.threshold)
                --threshold $advanced_settings.threshold
            #end if
            #if str($advanced_settings.min_probes)
                --min-probes $advanced_settings.min_probes
            #end if
            $advanced_settings.male_reference
            #if str($advanced_settings.Sample_sex.sex) == "yes":
                #if '$advanced_settings.Sample_sex.sample_sex' == "Male"
                    #set $sample_sex_val = "Male"
                    --sample-sex '$segment_method_val'
                #else
                    --sample-sex '$advanced_settings.Sample_sex.sample_sex'
                #end if
            #end if
            $advanced_settings.no_shift_xy
            #if $plot_aesthetics.title
                --title '$plot_aesthetics.title'
            #end if
            #if $advanced_settings.chromosome
                --chromosome '$advanced_settings.chromosome'
            #end if
            $plot_aesthetics.no_gene_labels
    ]]></command>
    <inputs>
        <param name="input_cnr_file" type="data" format="cnr" label="Bin-Level log2 Ratios/Coverages cnr file" help="" />
        <section name="advanced_settings" title="Advanced settings" expanded="false">
            <expand macro="diagram_optional" />
            <expand macro="sample_sex_condition" />
        </section>    
        <section name="plot_aesthetics" title="Plot aesthetics" expanded="false">
            <expand macro="diagram_plot" />
        </section>
    </inputs>
    <outputs>
        <data name="out_diagram_file" format="pdf" label="${tool.name} on ${on_string}: diagram pdf file" from_work_dir="sample-diagram.pdf" />
    </outputs>
    <tests>
        <test expect_num_outputs="1">
            <param name="input_cnr_file" ftype="cnr" value="tumor.cnr" />
            <section name="advanced_settings">
                <param name="segment" ftype="cns" value="sample.cns" />
                <param name="threshold" value="10" />
                <param name="male_reference" value="1" />
            </section>
            <output name="out_diagram_file" file="sample-diagram.pdf" ftype="pdf" compare="sim_size" > 
                <assert_contents><has_size value="6000" delta="2000" /></assert_contents>
            </output>
        </test>
        <test expect_num_outputs="1">
            <param name="input_cnr_file" ftype="cnr" value="tumor.cnr" />
            <section name="advanced_settings">
                <conditional name="Sample_sex">
                    <param name="sex" value="yes" />
                </conditional>
                <param name="segment" ftype="cns" value="sample.cns" />
                <param name="threshold" value="10" />
                <param name="male_reference" value="1" />
            </section>
            <output name="out_diagram_file" file="sample-diagram.pdf" ftype="pdf" compare="sim_size" > 
                <assert_contents><has_size value="6000" delta="2000" /></assert_contents>
            </output>
        </test>
    </tests>
    <help><![CDATA[
         Draw copy number (log2 coverages, segments) on chromosomes as a diagram. If both the raw probes and
         segments are given, show them side-by-side on each chromosome (segments on the left side, probes on the
         right side).

-----

**Bin-level log2 ratios (.cnr)**

Tabular file containing normalized log2 ratios for small genomic bins (divided regions of the genome). Used to detect raw copy number variations (CNVs) before segmentation.

.. csv-table::
   :header-rows: 0

    "chromosome","Genomic chromosome (e.g., chr1, chrX)"
    "start","Start position of the bin."
    "end","End position of the bin."
    "gene","Gene name(s) overlapping the bin (if applicable)."
    "log2","Normalized log2 ratio (sample coverage / reference coverage)."
    "depth","Average read depth in the bin."
    "weight","Reliability weight of the bin (higher = more reliable)."

-----

**Segmented log2 ratios (.cns)**

Tabular file with smoothed, merged segments of stable copy number, derived from the .cnr file. Represents final CNV calls.

.. csv-table::
   :header-rows: 0

    "chromosome","start, end: Genomic coordinates of the segment"
    "gene","Gene(s) overlapping the segment."
    "log2","Mean log2 ratio of the segment."
    "probes","Mean log2 ratio of the segment."
    "depth","Average read depth."
    "weight","Reliability weight."
    "p_value","Statistical confidence (lower = more significant)."

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