Mercurial > repos > jjohnson > mothur_toolsuite
annotate mothur/tools/mothur/shhh.flows.xml @ 31:a3eed59297ea
Patches courtesy of Peter Briggs, Bioinformatics Core Facility University of Manchester
make.contigs.xml.patch:# make.contigs.xml.patch
make.contigs.xml.patch:#
make.contigs.xml.patch:# 1. Fix cosmetic typo in <description> (forard -> forward)
make.contigs.xml.patch:# 2. Address error due to having 'mismatch' as the name for both an input and an output parameter:
make.contigs.xml.patch:# rename output parameter to 'cmismatch'
make.contigs.xml.patch:# 3. Remove 'threshold' parameter: make.contigs in mothur doesn't support a 'threshold' parameter
metagenomics.py.patch:# metagenomics.py.patch
metagenomics.py.patch:#
metagenomics.py.patch:# 1. Groups class: names were being taken from the wrong field (affected shhh.flows tool)
metagenomics.py.patch:# 2. Axes class: make 'sniff' method more sensitive to try and restrict arbitrary tabular
metagenomics.py.patch:# data uploads being sniffed as this type
mothur_wrapper.py.patch:# mothur_wrapper.py.patch
mothur_wrapper.py.patch:#
mothur_wrapper.py.patch:# 1. Update 'cmd_dict' settings for shhh.flows and shhh.seqs (otherwise these functions will
mothur_wrapper.py.patch:# fail on execution)
mothur_wrapper.py.patch:# 2. Fix add_option calls defining '--match' and '--mismatch' command line options (otherwise
mothur_wrapper.py.patch:# syntax error causes immediate failure)
screen.seqs.xml.patch:# screen.seqs.xml.patch
screen.seqs.xml.patch:#
screen.seqs.xml.patch:# Replace pattern for align.report output file in definiting of 'results' parameter in
screen.seqs.xml.patch:# <command> section (otherwise output_alignreport data item is empty).
shhh.flows.xml.patch:# shhh.flows.xml.patch
shhh.flows.xml.patch:#
shhh.flows.xml.patch:# Replace 'format_source' with 'format' for output parameters (otherwise formats are not
shhh.flows.xml.patch:# correctly assigned to output datasets)
shhh.seqs.xml.patch:# shhh.seqs.xml.patch
shhh.seqs.xml.patch:#
shhh.seqs.xml.patch:# 1. Fix patterns in --result (in <command> section) for shhh_seqs.fasta and shhh_seqs.names
shhh.seqs.xml.patch:# output files (otherwise files are not collected and associated data items are empty)
shhh.seqs.xml.patch:# 2. Replace 'format_source' with 'format' for output parameters (otherwise formats are not
shhh.seqs.xml.patch:# correctly assigned to output datasets)
trim.flows.xml.patch:# trim.flows.xml.patch
trim.flows.xml.patch:#
trim.flows.xml.patch:# Remove erroneous space from --result definition in <command> section (otherwise causes tool
trim.flows.xml.patch:# failure)
trim.seqs.xml.patch:# trim.seqs.xml.patch
trim.seqs.xml.patch:#
trim.seqs.xml.patch:# 1. Remove reference to undefined 'oligo.allvalues' varible in <command> section (otherwise
trim.seqs.xml.patch:# causes failure on execution)
trim.seqs.xml.patch:# 2. Fix format for input parameter 'names' (format should be 'names' not 'name')
trim.seqs.xml.patch:# 3. Add output parameter 'scrap_names' (to ensure consistent collection of all outputs)
trim.seqs.xml.patch:# 4. Update --result definition in <command> section to collect both trim.names and scrap.names
author | Jim Johnson <jj@umn.edu> |
---|---|
date | Tue, 30 Jul 2013 09:26:31 -0500 |
parents | 49058b1f8d3f |
children | ec8df51e841a |
rev | line source |
---|---|
27
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
1 <tool id="mothur_shhh_flows" name="Shhh.flows" version="1.26.0" force_history_refresh="True"> |
15
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
2 <description>Denoise flowgrams (PyroNoise algorithm)</description> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
3 <command interpreter="python"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
4 mothur_wrapper.py |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
5 #import re, os.path |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
6 --cmd='shhh.flows' |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
7 --result='^mothur.\S+\.logfile$:'$logfile,'^\S+\.shhh\.fasta$:'$shhh_fasta,'^\S+\.shhh\.qual$:'$shhh_qual,'^\S+\.shhh\.names$:'$shhh_names,'^\S+\.shhh\.groups$:'$shhh_groups,'^\S+\.shhh\.counts$:'$shhh_counts |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
8 --outputdir='$logfile.extra_files_path' |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
9 --flow=$flow |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
10 --lookup=$prob.lookup |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
11 #if $maxiter.__str__ != '': |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
12 --maxiter=$maxiter |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
13 #end if |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
14 #if $mindelta.__str__ != '': |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
15 --mindelta=$mindelta |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
16 #end if |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
17 #if $cutoff.__str__ != '': |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
18 --cutoff=$cutoff |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
19 #end if |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
20 #if $sigma.__str__ != '': |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
21 --sigma=$sigma |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
22 #end if |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
23 #if $order.__str__.strip() != '': |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
24 --order=$order |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
25 #end if |
27
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
26 #if $large.__str__ != '': |
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
27 --large=$large |
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
28 #end if |
15
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
29 --processors=8 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
30 </command> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
31 <inputs> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
32 <param name="flow" type="data" format="sff.flow" label="flow - flowgram data" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
33 help="Use sffinfo to generate flow data from an sff file and usually trimmed by trim.flows"/> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
34 <conditional name="prob"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
35 <param name="source" type="select" label="Select Taxonomy from" help=""> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
36 <option value="ref">Cached Reference</option> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
37 <option value="hist">History</option> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
38 </param> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
39 <when value="ref"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
40 <param name="lookup" type="select" format="tabular" label="lookup - intensity value per homopolymer length" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
41 help="table of the probability of observing an intensity value for a given homopolymer length"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
42 <options from_file="mothur_lookup.loc"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
43 <column name="name" index="0" /> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
44 <column name="value" index="1" /> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
45 </options> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
46 </param> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
47 </when> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
48 <when value="hist"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
49 <param name="lookup" type="data" format="tabular" label="lookup - intensity value per homopolymer length" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
50 help="from http://www.mothur.org/wiki/Lookup_files"/> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
51 </when> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
52 </conditional> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
53 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
54 <param name="maxiter" type="integer" value="1000" optional="true" label="maxiter - maximum iterations to run (default 1000)" help="if the delta value does not first drop below the mindelta value. "> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
55 <validator type="in_range" message="Number of differences can't be negative" min="0"/> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
56 </param> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
57 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
58 <param name="mindelta" type="float" value="" optional="true" label="mindelta - threshold for determining how much change in the flowgram correction is allowed" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
59 help="default .0000001 (10^-6)"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
60 <validator type="in_range" message="mindelta between 0. and .1" min="0.0" max="0.1"/> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
61 </param> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
62 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
63 <param name="cutoff" type="float" value="" optional="true" label="cutoff - seed the expectation-maximizaton step" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
64 help="default .01 (usually doesn't need to be changed)"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
65 <validator type="in_range" message="cutoff between 0. and 1." min="0.0" max="1.0"/> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
66 </param> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
67 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
68 <param name="sigma" type="float" value="" optional="true" label="sigma - the dispersion of the data in the expectation-maximization step of the algorithm" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
69 help="default .06 (usually doesn't need to be changed)"> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
70 <validator type="in_range" message="sigma between 0. and 1." min="0.0" max="1.0"/> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
71 </param> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
72 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
73 <param name="order" type="text" value="" label="order - flow order for nucleotides in the sequencer" |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
74 help="default is TACG"/> |
27
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
75 |
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
76 <param name="large" type="integer" value="10000" optional="true" label="large - split your flow file and process the pieces separately (default 10000)" help=""> |
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
77 <validator type="in_range" message="large value must be positive" min="1"/> |
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
78 </param> |
15
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
79 </inputs> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
80 <outputs> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
81 <data format="html" name="logfile" label="${tool.name} on ${on_string}: logfile" /> |
31
a3eed59297ea
Patches courtesy of Peter Briggs, Bioinformatics Core Facility University of Manchester
Jim Johnson <jj@umn.edu>
parents:
27
diff
changeset
|
82 <data format="fasta" name="shhh_fasta" label="${tool.name} on ${on_string}: shhh.fasta"/> |
a3eed59297ea
Patches courtesy of Peter Briggs, Bioinformatics Core Facility University of Manchester
Jim Johnson <jj@umn.edu>
parents:
27
diff
changeset
|
83 <data format="qual454" name="shhh_qual" label="${tool.name} on ${on_string}: shhh.qual"/> |
a3eed59297ea
Patches courtesy of Peter Briggs, Bioinformatics Core Facility University of Manchester
Jim Johnson <jj@umn.edu>
parents:
27
diff
changeset
|
84 <data format="names" name="shhh_names" label="${tool.name} on ${on_string}: shhh.names"/> |
a3eed59297ea
Patches courtesy of Peter Briggs, Bioinformatics Core Facility University of Manchester
Jim Johnson <jj@umn.edu>
parents:
27
diff
changeset
|
85 <data format="groups" name="shhh_groups" label="${tool.name} on ${on_string}: shhh.groups"/> |
a3eed59297ea
Patches courtesy of Peter Briggs, Bioinformatics Core Facility University of Manchester
Jim Johnson <jj@umn.edu>
parents:
27
diff
changeset
|
86 <data format="tabular" name="shhh_counts" label="${tool.name} on ${on_string}: shhh.counts"/> |
15
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
87 </outputs> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
88 <requirements> |
27
49058b1f8d3f
Update to mothur version 1.27 and add tool_dependencies.xml to automatically install mothur
Jim Johnson <jj@umn.edu>
parents:
15
diff
changeset
|
89 <requirement type="package" version="1.27">mothur</requirement> |
15
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
90 </requirements> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
91 <tests> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
92 </tests> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
93 <help> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
94 **mothur overview** |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
95 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
96 Mothur_, initiated by Dr. Patrick Schloss and his software development team |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
97 in the Department of Microbiology and Immunology at The University of Michigan, |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
98 provides bioinformatics for the microbial ecology community. |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
99 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
100 .. _Mothur: http://www.mothur.org/wiki/Main_Page |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
101 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
102 **Command Documenation** |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
103 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
104 The shhh.flows_ command is Pat Schloss's translation of Chris Quince's PyroNoise algorithm [1] from C to C++ with the incorporation of mothur's bells and whistles. Based on processing of test datasets provided by Quince, shhh.flows gives the same/similar output to AmpliconNoise. shhh.flows uses a expectation-maximization algorithm to correct flowgrams to identify the idealized form of each flowgram and translate that flowgram to a DNA sequence. Our testing has shown that when Titanium data are trimmed to 450 flows using trim.flows, shhh.flows provides the highest quality data for any other method available. In contrast, when we use the min/max number of flows suggested by Quince of 360/720, the error rate is not that great. This much improved error rate does come at a computational cost. Whereas the features in trim.seqs take on the order of minutes, shhh.flows can take on the order of hours. You will also need a lookup file that tells shhh.flows the probability of observing an intensity value for a given homopolymer length. You can get mothur-compatible files at: http://www.mothur.org/wiki/Lookup_files |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
105 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
106 .. _shhh.flows: http://www.mothur.org/wiki/Shhh.flows |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
107 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
108 |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
109 </help> |
a6189f58fedb
Mothur - updated for Mothur version 1.22.0
Jim Johnson <jj@umn.edu>
parents:
diff
changeset
|
110 </tool> |