Mercurial > repos > iuc > dada2_makesequencetable

comparison README.rst @ 6:8470de786d47 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/dada2 commit ea6c9c638e742c097b0ef294161eeea447c09e06
author iuc
date Fri, 30 Jun 2023 07:59:55 +0000
parents
children
comparison
equal deleted inserted replaced
5:e96bbb9e24ea 6:8470de786d47
1 Wrappers for the core functionality of the dada2 package https://benjjneb.github.io/dada2/index.html.
2
3 - assignTaxonomyAddspecies
4 - filterAndTrim
5 - learnErrors
6 - dada
7 - mergePairs
8 - makeSequenceTable
9 - plotComplexity
10 - plotQualityProfile
11 - removeBimeraDenovo
12 - seqCounts
13
14 Installation
15 ============
16
17 A conda release >=4.6 is needed so that the bioconductor-dada2 conda environment can be installed correctly.
18
19 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.
20 In order to run the tools you may either upgrade Galaxy or execute the following two steps:
21
22
23 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)
24
25 2. insert the following lines in ``config/datatypes.xml`` (just before the line ``</registration>``):
26 .. code-block::
27
28 <datatype extension="dada2_dada" type="galaxy.datatypes.binary:RData" subclass="true" display_in_upload="true" />
29 <datatype extension="dada2_errorrates" type="galaxy.datatypes.binary:RData" subclass="true" display_in_upload="true" />
30 <datatype extension="dada2_mergepairs" type="galaxy.datatypes.binary:RData" subclass="true" display_in_upload="true" />
31 <datatype extension="dada2_sequencetable" type="galaxy.datatypes.tabular:Tabular" mimetype="application/text" subclass="true" display_in_upload="true" />
32 <datatype extension="dada2_uniques" type="galaxy.datatypes.tabular:Tabular" mimetype="application/text" subclass="true" display_in_upload="true" />
33
34
35 Datatypes
36 =========
37
38 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.
39
40 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):
41
42 - dada2_dada (Rdata: named list, see docs for dada-class)
43 - dada2_errorrates (Rdata: named list, see docs for learnErrors)
44 - dada2_mergepairs (Rdata: named list, see docs for mergePairs)
45
46
47 For the outputs of makeSequenceTable and removeBimeraDenovo the following data types are used which derive from tabular:
48
49 - dada2_uniques
50 - in R a named integer vector (names are the unique sequences)
51 - in Galaxy written as a table (each row corresponding to a unique sequence, column 1: the sequence, column 2: the count)
52 - dada2_sequencetable
53 - in R a named integer matrix (rows = samples, columns = unique sequences)
54 - in Galaxy written as a table (rows = unique sequences, columns = samples)
55
56
57 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).
58
59 Test data
60 =========
61
62 Test data for ``dada2_seqCounts`` is generated using planemo's ``--update_test_data`` argument and manual
63 inspection of the test files. In addition a run of the pipeline (using collections) is executed
64 manually using ``planemo serve`` making sure that the entries of the tables are generated in a useful way.
65
66 In order to have the Collection unzip tool available use ``planemo s --galaxy_root GALAXY_ROOT --extra_tools GALAXY_ROOT/lib/galaxy/tools/``
67
68 All test other test data is generated using the shell script (``gentest.sh``) in test-data