Repository revision
repository tip
Select a revision to inspect and download versions of Galaxy utilities from this repository.

Repository dada2_learnerrors
Owner: iuc
Synopsis: Wrapper for DADA2 function: dada2: learnErrors.
Fast and accurate sample inference from amplicon data with single-nucleotide resolution
Type: unrestricted
Revision: 7:75403243703a
This revision can be installed: True
Times cloned / installed: 583

Repository README files - may contain important installation or license information

Wrappers for the core functionality of the dada2 package https://benjjneb.github.io/dada2/index.html.

  • assignTaxonomyAddspecies
  • filterAndTrim
  • learnErrors
  • dada
  • mergePairs
  • makeSequenceTable
  • plotComplexity
  • plotQualityProfile
  • removeBimeraDenovo
  • seqCounts

Installation

A conda release >=4.6 is needed so that the bioconductor-dada2 conda environment can be installed correctly.

The dada2 wrappers can be installed via the toolshed. Since they use datatypes that have been introduced with Galaxy release 19.09 they won't work out of the box for older Galaxy releases. In order to run the tools you may either upgrade Galaxy or execute the following two steps:

  1. find GALAXY_ROOT/shed_tools/testtoolshed.g2.bx.psu.edu/repos/iuc/ -name "dada2_*xml" -exec sed -i -e 's/profile="19.09"/profile="YOUR_RELEASE"/' {} ; (replace GALAXY_ROOT and YOUR_RELEASE appropriately)

  2. insert the following lines in config/datatypes.xml (just before the line </registration>):
    <datatype extension="dada2_dada" type="galaxy.datatypes.binary:RData" subclass="true" display_in_upload="true" />
    <datatype extension="dada2_errorrates" type="galaxy.datatypes.binary:RData" subclass="true" display_in_upload="true" />
    <datatype extension="dada2_mergepairs" type="galaxy.datatypes.binary:RData" subclass="true" display_in_upload="true" />
    <datatype extension="dada2_sequencetable" type="galaxy.datatypes.tabular:Tabular" mimetype="application/text" subclass="true" display_in_upload="true" />
    <datatype extension="dada2_uniques" type="galaxy.datatypes.tabular:Tabular" mimetype="application/text" subclass="true" display_in_upload="true" />
    

Datatypes

The dada2 Galaxy wrappers use a few extra data types to ensure that only inputs of the correct type can be used, these datatypes are available from Galaxy release 19.05, for earlier releases they need to be added manually.

For the outputs of dada, learnErrors, and mergePairs the following datatypes are used that derive from Rdata (which contains the named list that is returned from the corresponding dada function):

  • dada2_dada (Rdata: named list, see docs for dada-class)
  • dada2_errorrates (Rdata: named list, see docs for learnErrors)
  • dada2_mergepairs (Rdata: named list, see docs for mergePairs)

For the outputs of makeSequenceTable and removeBimeraDenovo the following data types are used which derive from tabular:

  • dada2_uniques
    • in R a named integer vector (names are the unique sequences)
    • in Galaxy written as a table (each row corresponding to a unique sequence, column 1: the sequence, column 2: the count)
  • dada2_sequencetable
    • in R a named integer matrix (rows = samples, columns = unique sequences)
    • in Galaxy written as a table (rows = unique sequences, columns = samples)

Note the difference between the R and Galaxy representations! The main motivation is that the dada2_sequencetable is analogous to OTU tables as produced for instance by qiime (and it seemed natural to extend this to the uniques which are essentially a sequencetables of single samples).

Test data

Test data for dada2_seqCounts is generated using planemo's --update_test_data argument and manual inspection of the test files. In addition a run of the pipeline (using collections) is executed manually using planemo serve making sure that the entries of the tables are generated in a useful way.

In order to have the Collection unzip tool available use planemo s --galaxy_root GALAXY_ROOT  --extra_tools GALAXY_ROOT/lib/galaxy/tools/

All test other test data is generated using the shell script (gentest.sh) in test-data

Contents of this repository

Name Description Version Minimum Galaxy Version
Learn Error rates 1.28+galaxy0 19.09

Categories
Metagenomics - Tools enabling the study of metagenomes