view quickmerge.xml @ 0:436c1d3c990a draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/quickmerge commit c171c7bde4063ec25fff4a284f5d0330d46ca290

<tool
  id="quickmerge"
  name="QuickMerge"
  version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@"
  profile="21.05"
>
  <description>merge long-read and hybrid genome assemblies to increase contiguity</description>

  <macros>
    <token name="@TOOL_VERSION@">0.3</token>
    <token name="@VERSION_SUFFIX@">0</token>
  </macros>

  <xrefs>
    <xref type="bio.tools">quickmerge</xref>
  </xrefs>

  <requirements>
    <requirement type="package" version="@TOOL_VERSION@">quickmerge</requirement>
    <requirement type="package" version="0.1.4">fastkit</requirement>
  </requirements>

  <command detect_errors="exit_code"><![CDATA[

## Strip spaces from FASTA header
fastkit format --strip-header-space '$input_ref' > input.ref.fasta
&& fastkit format --strip-header-space '$input_query' > input.query.fasta

## Nucmer step
## ----------------------------------------------------------------------------

&& nucmer
  $advanced.nucmer_match_protocol
  -l $advanced.nucmer_min_len
  input.ref.fasta
  input.query.fasta


## Delta-filter step
## ----------------------------------------------------------------------------

&& delta-filter
  $advanced.delta_map_ref
  $advanced.delta_map_query
  -l $advanced.delta_min_len

out.delta > out.rq.delta


## Quickmerge step
## ----------------------------------------------------------------------------

&& quickmerge
  -d out.rq.delta
  -q input.query.fasta
  -r input.ref.fasta
  -p out
  -hco $advanced.hco
  -c $advanced.c
  -l $advanced.l
  -ml $advanced.ml

  ]]></command>

  <inputs>
    <param
      name="input_query"
      type="data"
      format="fasta,fasta.gz"
      label="Query (hybrid) assembly"
      help="A hydrid or long-read assembly."
    />
    <param
      name="input_ref"
      type="data"
      format="fasta,fasta.gz"
      label="Reference assembly"
      help="The self assembly - this can also be a long-read assembly."
    />

    <section name="advanced" title="Advanced parameters" expanded="false">
      <param
        name="nucmer_match_protocol"
        type="select"
        display="radio"
        label="Anchor matching protocol"
        help="By default, nucmer will include anchor matches that are unique
          in the reference but not necessarily unique in the query. You can
          change the protocol to include all matches, or only those unique to
          both the reference and query sequence."
      >
        <option value="" selected="true">
          Default (unique in reference only)
        </option>
        <option value="--mum">
          Include matches that are unique to both reference and query
        </option>
        <option value="--maxmatch">
          Include all matches
        </option>
      </param>

      <param
        name="nucmer_min_len"
        type="integer"
        label="Nucmer -l"
        help="Minimum length for a single match (nt)"
        value="20"
        min="0"
      />

      <param
        name="delta_min_len"
        type="integer"
        label="Delta-filter -l"
        help="Minimum alignment length (nt)"
        value="0"
        min="0"
      />

      <param
        name="delta_map_ref"
        type="boolean"
        truevalue="-r"
        falsevalue=""
        checked="true"
        label="Delta-filter -r"
        help="Maps each position of each reference to its best hit in the
          query, allowing for query overlaps."
      />

      <param
        name="delta_map_query"
        type="boolean"
        truevalue="-q"
        falsevalue=""
        checked="true"
        label="Delta-filter -q"
        help="Maps each position of each query to its best hit in the
          reference, allowing for reference overlaps."
      />

      <param
        type="integer"
        argument="-l"
        help="Minimum seed contig length to be merged (nt). A good rule of
          thumb is to start with the N50 of the self assembly."
        value="0"
        min="0"
      />

      <param
        type="integer"
        argument="-ml"
        help="Set the merging length cutoff (nt). This is especially helpful
          for repeat-rich genomes, where a higher minimum length can reduce
          errors arising from long repeats."
        value="5000"
        min="0"
      />

      <param
        type="float"
        argument="-hco"
        help="Controls the overlap cutoff used in selection of anchor contigs.
          Increasing this will raise stringency, resulting in fewer joins but
          higher confidence."
        value="5.0"
        min="0.0"
      />

      <param
        type="float"
        argument="-c"
        help="Controls the overlap cutoff for contigs used for extension of
          the anchor contig. Increasing this will raise stringency,
          resulting in fewer joins but higher confidence."
        value="1.5"
        min="0.0"
      />
    </section>

    <section name="output" title="Optional outputs" expanded="true">
      <param
        name="check_align_summary"
        type="boolean"
        label="Alignment summary table"
      />
      <param
        name="check_anchor_summary"
        type="boolean"
        label="Anchor summary table"
      />
      <param
        name="check_param_summary"
        type="boolean"
        label="Parameter summary table"
      />
    </section>
  </inputs>

  <outputs>
    <data
      name="merged_fasta"
      format="fasta"
      from_work_dir="merged_out.fasta"
      label="${tool.name} on ${on_string}: Merged assembly"
    />

    <data
      name="align_summary"
      format="tabular"
      from_work_dir="aln_summary_out.tsv"
      label="${tool.name} on ${on_string}: Alignment summary"
    >
      <filter>output['check_align_summary']</filter>
    </data>

    <data
      name="anchor_summary"
      format="tabular"
      from_work_dir="anchor_summary_out.txt"
      label="${tool.name} on ${on_string}: Anchor summary"
    >
      <filter>output['check_anchor_summary']</filter>
    </data>

    <data
      name="param_summary"
      format="tabular"
      from_work_dir="param_summary_out.txt"
      label="${tool.name} on ${on_string}: Parameter summary"
    >
      <filter>output['check_param_summary']</filter>
    </data>
  </outputs>

  <tests>
    <!-- Standard run -->
    <test expect_num_outputs="1">
      <param name="input_query" value="ecoli_nanopore.fasta" ftype="fasta" />
      <param name="input_ref" value="ecoli_illumina.fasta" ftype="fasta" />
      <output name="merged_fasta" ftype="fasta" file="merged.fasta" />
    </test>

    <!-- Standard run with additional outputs -->
    <test expect_num_outputs="4">
      <param name="input_query" value="ecoli_nanopore.fasta" ftype="fasta" />
      <param name="input_ref" value="ecoli_illumina.fasta" ftype="fasta" />
      <param name="check_align_summary" value="true" />
      <param name="check_anchor_summary" value="true" />
      <param name="check_param_summary" value="true" />
      <output name="merged_fasta" ftype="fasta" file="merged.fasta" />
      <output name="align_summary" ftype="tabular" file="2_aln_summary_out.tsv" />
      <output name="anchor_summary" ftype="tabular" file="2_anchor_summary_out.txt" />
      <output name="param_summary" ftype="tabular" file="2_param_summary_out.txt" />
    </test>

    <!-- Standard run with advanced params set -->
    <test expect_num_outputs="1">
      <param name="input_query" value="ecoli_nanopore.fasta" ftype="fasta" />
      <param name="input_ref" value="ecoli_illumina.fasta" ftype="fasta" />
      <param name="advanced.nucmer_min_len" value="20" />
      <param name="advanced.delta_min_len" value="0" />
      <param name="advanced.hco" value="0" />
      <param name="advanced.c" value="5000" />
      <param name="advanced.l" value="1.5" />
      <param name="advanced.ml" value="5.0" />
      <output name="merged_fasta" ftype="fasta" file="merged.fasta" />
    </test>
  </tests>

  <help><![CDATA[

.. class:: infomark

**What it does**

The program uses complementary information from genomes assembled with long reads in order to improve contiguity, and works with assemblies derived from both Pacific Biosciences or Oxford Nanopore. Quickmerge will even work with hybrid assemblies made by combining long reads and Illumina short reads. This can be useful when long read coverage is limiting. For more details, please see the paper: (`10.1093/nar/gkw654 <https://doi.org/10.1093/nar/gkw654>`_) that describes it.

Although this program was written to merge a hybrid assembly (e.g. as generated by DBG2OLC) and a PacBio or ONP only assembly, it can also be used to merge two different long molecule only assemblies (e.g. one generated with PBcR or canu and another generated with FALCON).

Please consult the `QuickMerge Wiki <https://github.com/mahulchak/quickmerge/wiki>`_ for extended usage information.

**Why use quickmerge?**

- More bang for your buck by using complementary information from different assemblers.
- It is fast. Takes less than a minute to run on most genomes. You run nucmer once (nucmer is the most time consuming step) and then you can run quickmerge over a large number of parameters in a very short time.
- Requires only FASTA files and does not depend on any special data or computational resources.
- Saves money. When long read coverage is limiting, there are some hybrid approaches that lead to good results (e.g. DBG2OLC). So quickmerge allows you to cut your long molecule requirement by as much as half by replacing the same with Illumina short reads. E.g. if you think you would get a N50 of 8Mb from 75X long reads (ONP or PacBio), try sequencing 45X long and 70X Illumina reads instead of 75X long reads. You may not need that extra 35X long reads.
- Allows dramatic improvements with reanalysis of legacy data collected when long reads were prohibitively expensive. See an excellent example from `this paper by Thomas Mather <https://academic.oup.com/g3journal/article/10/3/899/6026189>`_ on the Soybean Aphid genome.


**The process**

- **Nucmer** aligns the two assemblies so that the merger can find the correct splice sites
- **Delta-filter** filters out alignments with repeats and duplicates
- **Quickmerge** merges the aligned assemblies

**Inputs**

- **Query assembly**: a hydrid or long-read assembly (see `QuickMerge Wiki <https://github.com/mahulchak/quickmerge/wiki>`_ for more details)

- **Reference assembly**: "Self" assembly. Can also be a hybrid assembly

- **Nucmer -l**: the minimum length of a single match

- **Delta-filter -l**: the minimum length of a single alignment

- **hco**: controls the overlap cutoff used in selection of anchor contigs (default 5.0)

- **c**: controls the overlap cutoff for contigs used for extension of the anchor contig (default 1.5)

- **l**: controls the length cutoff for anchor contigs. A good rule of thumb is to start with the N50 of the self assembly. E.g. if the N50 of your self assembly is 2Mb then use 2000000 as your cutoff. Lowering this value may lead to more merging but may increase the probability of mis-joins.

- **ml**: controls the minimum alignment length to be considered for merging. This is especially helpful for repeat-rich genomes (default 5000; higher recommended).

For both ``hco`` and ``c``, bigger the number, more stringent is the criteria for contig selection (which will lead to fewer contigs being merged). If they are too small (<1), chances of spurious merging will increase. It is better to be conservative while merging contigs!


**Outputs**

- **Merged assembly**: your merged assembly in FASTA format

- **Alignment summary**: tabular summary of alignments used for merging

- **Anchor summary**: tabular summary of anchor matches used for merging

- **Parameter summary**: summary of merge parameters


**Helpful tips**

- For optimal merging results, identify the major misassemblies (especially translocations and inversions) in the component assemblies and break the contigs at such misassembly boundaries. Alignment of the component assemblies to the merged assembly may help to identify such assembly errors because a specific error typically occurs in only one of the assemblies.

- The FASTA sequence headers *should not have white spaces* in them. In case they do, as might happen for assemblies obtained from FALCON assembler, the white space needs to be removed before running.

- You can run Ka-kit's `finisherSC <https://github.com/kakitone/finishingTool>`_ after running quickmerge to improve the contiguity even further.

- Assembly polishing with Quiver and pilon before and after assembly merging is strongly recommended. However, if you are running finisherSC, you may perform the quiver polishing after the finisher step.

- Check the merged assembly by aligning the hybrid and/or non-hybrid assembly to the merged assembly. You can use ``nucmer`` for alignment (with default anchoring) and ``mummerplot`` for dot plot visualization.


Quickmerge was wrapped by the Galaxy Australia team.

  ]]></help>
  <citations>
    <citation type="doi">https://doi.org/10.1093/nar/gkw654</citation>
  </citations>
</tool>