comparison lr.xml @ 1:d5124d5c8131 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/delly commit 3da28c7772346e1872b6d768b904305be0c61db7"

0:629a0066003d

<?xml version="1.0"?>
<tool id="delly_lr" name="Delly long-read (lr)" version="@TOOL_VERSION@+galaxy0" profile="18.01">
    <description>optimized calling and genotyping of structural variants</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro="requirements"/>
    <expand macro="version_command"/>
    <command detect_errors="exit_code"><![CDATA[
## initialize
#for $i, $current in enumerate($samples)
    ln -s '${current}' 'sample_${i}.bam' &&
    ln -s '${current.metadata.bam_index}' 'sample_${i}.bam.bai' &&
#end for
        
## run
delly lr
## generic options
--svtype $generic.svtype

#end if
--outfile 'result.bcf'
## discovery options
--mapqual $discovery.mapqual
--minclip $discovery.minclip
--minrefsep $discovery.minrefsep
--maxreadsep $discovery.maxreadsep
## genotyping options
#if $genotyping.vcffile
    --vcffile '$genotyping.vcffile'

#for $i, $current in enumerate($samples)
    'sample_${i}.bam'
#end for

## postprocessing
#if 'log' in $oo.out
    |& tee 'log.txt'
#end if
#if 'vcf' in $oo.out
    && test -f 'result.bcf' && bcftools view 'result.bcf' > 'result.vcf' || echo 'No results.'
#end if
#if 'dump' in $oo.out
    && test -f 'dump.tsv.gz' && bgzip -d 'dump.tsv.gz' || echo 'No dump file.'
#end if
    ]]></command>
    <inputs>
        <expand macro="samples"/>
        <section name="generic" title="Generic options" expanded="true">
            <expand macro="genome"/>

            </param>
        </section>
        <section name="discovery" title="Discovery options" expanded="true">
            <param argument="--mapqual" type="integer" value="1" label="Set minimum mapping quality"/>
            <expand macro="minclip"/>
            <expand macro="minrefsep"/>
            <expand macro="maxreadsep"/>
        </section>
        <section name="genotyping" title="Genotyping options" expanded="true">
            <expand macro="vcffile"/>
            <expand macro="genoqual"/>
        </section>

            <section name="discovery">
                <param name="mapqual" value="2"/>
                <param name="qualtra" value="19"/>
                <param name="madcutoff" value="8"/>
                <param name="minclip" value="24"/>
                <param name="minrefsep" value="24"/>
                <param name="maxreadsep" value="39"/>
            </section>
            <section name="genotyping">
                <param name="genoqual" value="4"/>

**What it does**

@WID@

Delly *lr* uses the long-read SV discovery mode.

**Input**

Delly *lr* needs a sorted, indexed and duplicate marked BAM file for every input sample. An indexed reference genome is required to identify split-reads. Additionally a VCF/BCF file for genotyping can be applied.

**Output**

The output is available in BCF and VCF format. Additionally a output file for SV-reads is provided.

.. class:: infomark

**References**

1:d5124d5c8131

<?xml version="1.0"?>
<tool id="delly_lr" name="Delly long-read (lr)" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="18.01">
    <description>optimized calling and genotyping of structural variants</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro="requirements"/>
    <expand macro="version_command"/>
    <command detect_errors="exit_code"><![CDATA[
## initialize
@BAM@
        
## run
delly lr
## generic options
--svtype $generic.svtype

#end if
--outfile 'result.bcf'
## discovery options
--mapqual $discovery.mapqual
--minclip $discovery.minclip
--min-clique-size $discovery.mincliquesize
--minrefsep $discovery.minrefsep
--maxreadsep $discovery.maxreadsep
## genotyping options
#if $genotyping.vcffile
    --vcffile '$genotyping.vcffile'

#for $i, $current in enumerate($samples)
    'sample_${i}.bam'
#end for

## postprocessing
@LOG@
@VCF@
@DUMP@
    ]]></command>
    <inputs>
        <expand macro="samples"/>
        <section name="generic" title="Generic options" expanded="true">
            <expand macro="genome"/>

            </param>
        </section>
        <section name="discovery" title="Discovery options" expanded="true">
            <param argument="--mapqual" type="integer" value="1" label="Set minimum mapping quality"/>
            <expand macro="minclip"/>
            <expand macro="mincliquesize"/>
            <expand macro="minrefsep" defaut="30"/>
            <expand macro="maxreadsep" defaut="75"/>
        </section>
        <section name="genotyping" title="Genotyping options" expanded="true">
            <expand macro="vcffile"/>
            <expand macro="genoqual"/>
        </section>

            <section name="discovery">
                <param name="mapqual" value="2"/>
                <param name="qualtra" value="19"/>
                <param name="madcutoff" value="8"/>
                <param name="minclip" value="24"/>
                <param name="mincliquesize" value="1"/>
                <param name="minrefsep" value="24"/>
                <param name="maxreadsep" value="39"/>
            </section>
            <section name="genotyping">
                <param name="genoqual" value="4"/>

**What it does**

@WID@

Delly *long-read (lr)* uses the long-read SV discovery mode.

**Input**

Delly *long-read (lr)* needs a sorted, indexed and duplicate marked BAM file for every input sample. An indexed reference genome is required to identify split-reads. Additionally a VCF/BCF file for genotyping can be applied.

**Output**

The output is available in BCF and VCF format. Additionally an output file for SV-reads is provided.

.. class:: infomark

**References**