view drep_dereplicate.xml @ 0:8dfcdbeaeed8 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/drep commit 8fa5ff35b45c2b046c7f4800410cf39cb89a299a"

<tool id="drep_dereplicate" name="dRep dereplicate" version="@VERSION@.0" python_template_version="3.5">
    <description>De-replicate a list of genomes</description>
    <macros>
        <import>macros.xml</import>
    </macros>
    <expand macro="requirements" />
    <command detect_errors="exit_code"><![CDATA[
         @PREPARE_GENOMES@
         dRep dereplicate outdir
         @FILTER_OPTIONS@
         @GENOME_COMPARISON_OPTIONS@
         @CLUSTERING_OPTIONS@
         @SCORING_OPTIONS@
         @TAXONOMY_OPTIONS@
         @WARNING_OPTIONS@        
         @GENOMES@
         --debug
         || (rc=\$?; 
             ls -ltr `find outdir -type f`;
             cat outdir/data/checkM/checkM_outdir/checkm.log;
             cat outdir/log/logger.log;
             exit \$rc)
    ]]></command>
    <inputs>
        <expand macro="genomes"/>
        <expand macro="filtering_options"/>
        <expand macro="genome_comparison_options"/>
        <expand macro="clustering_options"/>
        <expand macro="scoring_options"/>
        <expand macro="taxonomy_options"/>
        <expand macro="warning_options"/>
        <expand macro="select_drep_outputs"/>
    </inputs>
    <outputs>
        <collection name="dereplicated_genomes" type="list" label="dereplicated_genomes">
             <discover_datasets pattern="__designation__" directory="outdir/dereplicated_genomes" ext='fasta'/>
        </collection>
        <expand macro="drep_outputs" />
    </outputs>
    <tests>
        <expand macro="test_defaults_log">
            <has_text text="dRep dereplicate finished" />
        </expand>
        <test>
            <param name="genomes" ftype="fasta" value="Enterococcus_casseliflavus_EC20.fasta,Enterococcus_faecalis_T2.fna,Enterococcus_faecalis_TX0104.fa"/>
            <conditional name="filter">
                <param name="set_options" value="yes"/>
                <conditional name="quality">
                    <param name="source" value="checkm"/>
                    <param name="checkM_method" value="taxonomy_wf"/>
                </conditional>
            </conditional>
            <output name="log">
                <assert_contents>
                    <has_text text="dRep dereplicate finished" />
                </assert_contents>
            </output>
        </test>
    </tests>
    <help><![CDATA[
**dRep dereplicate**

`dRep <https://drep.readthedocs.io/en/latest/overview.html>`_ performs rapid pair-wise comparison of genome sets.





`De-replication <https://drep.readthedocs.io/en/latest/overview.html#genome-de-replication>`_ is the process of identifying sets of genomes that are the “same” in a list of genomes, and removing all but the “best” genome from each redundant set. How similar genomes need to be to be considered “same”, how to determine which genome is “best”, and other important decisions are discussed in `Choosing parameters. <https://drep.readthedocs.io/en/latest/choosing_parameters.html>`_   Detailed options for each module are described at: https://drep.readthedocs.io/en/latest/module_descriptions.html

A common use for genome de-replication is the case of individual assembly of metagenomic data. If metagenomic samples are collected in a series, a common way to assemble the short reads is with a “co-assembly”. That is, combining the reads from all samples and assembling them together. The problem with this is assembling similar strains together can severely fragment assemblies, precluding recovery of a good genome bin. An alternative option is to assemble each sample separately, and then “de-replicate” the bins from each assembly to make a final genome set.

The steps to this process are:

  - Assemble each sample separately using your favorite assembler. You can also perform a co-assembly to catch low-abundance microbes
  - Bin each assembly (and co-assembly) separately using your favorite binner
  - Pull the bins from all assemblies together and run **dRep** on them
  - Perform downstream analysis on the de-replicated genome list



**INPUTS**

  - Genome sets in fasta format.


**OUTPUTS**

  - `Figures <https://drep.readthedocs.io/en/latest/example_output.html#figures>`_ that show the relationship of the Genome inputs.
  - `Warnings <https://drep.readthedocs.io/en/latest/example_output.html#warnings>`_ report two things: de-replicated genome similarity and secondary clusters that were almost different. 
  - A Dataset collection of the “best” genome of each secondary cluster.
  - `Tables from intermediate steps <https://drep.readthedocs.io/en/latest/advanced_use.html>`_ 

    * Chdb.csv # CheckM results for Bdb
    * Widb.csv # Winning genomes' checkM information 


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