# HG changeset patch # User Jim Johnson # Date 1312907967 18000 # Node ID 370b3fc4e7d397adbba71220a194302f1dfe1fa2 # Parent 11867a45a43d5c8c87e3b67c13074c99d7065299 Mothur galaxy tool updates for Mothur version 1.21.0, fix help section warnings diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/chimera.slayer.xml --- a/mothur/tools/mothur/chimera.slayer.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/chimera.slayer.xml Tue Aug 09 11:39:27 2011 -0500 @@ -1,4 +1,4 @@ - + Find putative chimeras using slayer mothur_wrapper.py @@ -32,6 +32,7 @@ $options.trim $options.split #end if + ## --blastlocation=path_to_blast --processors=2 @@ -115,7 +116,7 @@ (A) the ends of a query sequence are searched against an included database of reference chimera-free 16S sequences to identify potential parents of a chimera; (B) candidate parents of a chimera are selected as those that form a branched best scoring alignment to the NAST-formatted query sequence; - (C) the NAST alignment of the query sequence is improved in a `chimera-aware' profile-based NAST realignment to the selected reference parent sequences; and + (C) the NAST alignment of the query sequence is improved in a 'chimera-aware' profile-based NAST realignment to the selected reference parent sequences; and (D) an evolutionary framework is used to flag query sequences found to exhibit greater sequence homology to an in silico chimera formed between any two of the selected reference parent sequences. Note: diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/classify.seqs.xml --- a/mothur/tools/mothur/classify.seqs.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/classify.seqs.xml Tue Aug 09 11:39:27 2011 -0500 @@ -1,4 +1,4 @@ - + Assign sequences to taxonomy mothur_wrapper.py @@ -97,6 +97,8 @@ + + diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/get.groups.xml --- a/mothur/tools/mothur/get.groups.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/get.groups.xml Tue Aug 09 11:39:27 2011 -0500 @@ -1,4 +1,4 @@ - + Select groups mothur_wrapper.py @@ -9,43 +9,61 @@ --cmd='get.groups' --outputdir='$logfile.extra_files_path' --group=$group_in - #if $groups.__str__ != "None" and len($groups.__str__) > 0: - --groups=$groups - #end if - #if $accnos.__str__ != "None" and len($accnos.__str__) > 0: - --accnos=$accnos + #if $groupnames.source == 'groups': + #if $groupnames.groups.__str__ != "None" and len($groupnames.groups.__str__) > 0: + --groups=$groupnames.groups + #end if + #else + #if $groupnames.accnos.__str__ != "None" and len($groupnames.accnos.__str__) > 0: + --accnos=$groupnames.accnos + #end if #end if #if $fasta_in.__str__ != "None" and len($fasta_in.__str__) > 0: --fasta=$fasta_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($fasta_in.__str__)) + ":'" + $fasta_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($fasta_in.__str__)) + ":'" + $fasta_out.__str__] #end if #if $name_in.__str__ != "None" and len($name_in.__str__) > 0: --name=$name_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($name_in.__str__)) + ":'" + $name_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($name_in.__str__)) + ":'" + $name_out.__str__] #end if #if $list_in.__str__ != "None" and len($list_in.__str__) > 0: --list=$list_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($list_in.__str__)) + ":'" + $list_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($list_in.__str__)) + ":'" + $list_out.__str__] + #end if + #if $shared_in.__str__ != "None" and len($shared_in.__str__) > 0: + --shared=$shared_in + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.*\\.pick\\.\2',$os.path.basename($shared_in.__str__)) + ":'" + $shared_out.__str__] #end if #if $taxonomy_in.__str__ != "None" and len($taxonomy_in.__str__) > 0: --taxonomy=$taxonomy_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($taxonomy_in.__str__)) + ":'" + $taxonomy_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($taxonomy_in.__str__)) + ":'" + $taxonomy_out.__str__] #end if --result=#echo ','.join($results) - - - - - - - - - + + + + + + + + + + + + + + + + + + + + @@ -61,6 +79,9 @@ list_in != None + + shared_in != None + taxonomy_in != None @@ -81,8 +102,12 @@ **Command Documenation** -The get.groups_ command selects sequences from a specific group or set of groups from the following file types: fasta, name, group, list, taxonomy. +The get.groups_ command selects sequences from a specific group or set of groups from the following file types: fasta, fasta, name_, group_, list_, taxonomy_. +.. _name: http://www.mothur.org/wiki/Name_file +.. _group: http://www.mothur.org/wiki/Group_file +.. _list: http://www.mothur.org/wiki/List_file +.. _taxonomy: http://www.mothur.org/wiki/Taxonomy_outline .. _get.groups: http://www.mothur.org/wiki/Get.groups diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/get.seqs.xml --- a/mothur/tools/mothur/get.seqs.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/get.seqs.xml Tue Aug 09 11:39:27 2011 -0500 @@ -102,6 +102,5 @@ .. _list.seqs: http://www.mothur.org/wiki/list.seqs .. _get.seqs: http://www.mothur.org/wiki/Get.seqs - diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/get.sharedseqs.xml --- a/mothur/tools/mothur/get.sharedseqs.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/get.sharedseqs.xml Tue Aug 09 11:39:27 2011 -0500 @@ -98,7 +98,7 @@ **Command Documenation** -The get.sharedseqs_ command takes a list and group file and outputs a *.shared.seqs file for each distance. This is useful for those cases where you might be interested in identifying sequences that are either unique or shared by specific groups, which you could then classify. +The get.sharedseqs_ command takes a list and group file and outputs a .shared.seqs file for each distance. This is useful for those cases where you might be interested in identifying sequences that are either unique or shared by specific groups, which you could then classify. .. _get.sharedseqs: http://www.mothur.org/wiki/Get.sharedseqs diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/indicator.xml --- a/mothur/tools/mothur/indicator.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/indicator.xml Tue Aug 09 11:39:27 2011 -0500 @@ -1,4 +1,4 @@ - + Identify indicator "species" for nodes on a tree mothur_wrapper.py @@ -6,7 +6,9 @@ #import re, os.path --cmd='indicator' --outputdir='$logfile.extra_files_path' - --tree=$tree + #if $tree.__str__ != "None" and len($tree.__str__) > 0: + --tree=$tree + #end if #if isinstance($otu.datatype, $__app__.datatypes_registry.get_datatype_by_extension('shared').__class__): --shared=$otu #elif isinstance($otu.datatype, $__app__.datatypes_registry.get_datatype_by_extension('relabund').__class__): @@ -24,7 +26,6 @@ --result='^mothur.\S+\.logfile$:'$logfile,'^\S+\.indicator.tre$:'$tree_out,'^\S+\.indicator.summary$:'$summary - @@ -42,13 +43,17 @@ + - - + + tree != None + + mothur @@ -66,9 +71,10 @@ **Command Documenation** -The indicator_ command reads a shared_ or relabund file and a tree file, and outputs a .indicator.tre and .indicator.summary file. The new tree contains labels at each internal node. The label is the node number so you can relate the tree to the summary file. The summary file lists the indicator value for each OTU for each node. +The indicator_ command reads a shared_ or relabund_ file and a tree file, and outputs a .indicator.summary file and when a tree file is given a .indicator.tre file. The summary file lists the indicator value for each OTU for each node. The new tree contains labels at each internal node. The label is the node number so you can relate the tree to the summary file. .. _shared: http://www.mothur.org/wiki/Shared_file +.. _relabund: http://www.mothur.org/wiki/Get.relabund .. _indicator: http://www.mothur.org/wiki/Indicator diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/mothur_wrapper.py --- a/mothur/tools/mothur/mothur_wrapper.py Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/mothur_wrapper.py Tue Aug 09 11:39:27 2011 -0500 @@ -176,21 +176,22 @@ """ cmd_dict = dict() cmd_dict['align.check'] = dict({'required' : ['fasta','map']}) - cmd_dict['align.seqs'] = dict({'required' : ['fasta','reference',], 'optional' : ['search','ksize','align','match','mismatch','gapopen','gapextend','flip','threshold','processors']}) + cmd_dict['align.seqs'] = dict({'required' : ['fasta','reference',], 'optional' : ['search','ksize','align','match','mismatch','gapopen','gapextend','flip','threshold','save','processors']}) cmd_dict['amova'] = dict({'required' : ['phylip','design'] , 'optional' : ['alpha','iters']}) cmd_dict['anosim'] = dict({'required' : ['phylip','design'] , 'optional' : ['alpha','iters']}) cmd_dict['bin.seqs'] = dict({'required' : ['list','fasta'], 'optional' : ['name','label','group']}) cmd_dict['bootstrap.shared'] = dict({'required' : ['shared'], 'optional' : ['calc','groups','iters','label']}) #catchall cmd_dict['chimera.bellerophon'] = dict({'required' : ['fasta'], 'optional' : ['filter','correction','window','increment','processors']}) - cmd_dict['chimera.ccode'] = dict({'required' : ['fasta','reference'], 'optional' : ['filter','mask','window','numwanted','processors']}) - cmd_dict['chimera.check'] = dict({'required' : ['fasta','reference'], 'optional' : ['ksize','svg','name','increment','processors']}) - cmd_dict['chimera.pintail'] = dict({'required' : ['fasta','reference'], 'optional' : ['conservation','quantile','filter','mask','window','increment','processors']}) - cmd_dict['chimera.slayer'] = dict({'required' : ['fasta','reference'], 'optional' : ['name','search','window','increment','match','mismatch','numwanted','parents','minsim','mincov','iters','minbs','minsnp','divergence','realign','split','processors']}) + cmd_dict['chimera.ccode'] = dict({'required' : ['fasta','reference'], 'optional' : ['filter','mask','window','numwanted','save','processors']}) + cmd_dict['chimera.check'] = dict({'required' : ['fasta','reference'], 'optional' : ['ksize','svg','name','increment','save','processors']}) + cmd_dict['chimera.pintail'] = dict({'required' : ['fasta','reference'], 'optional' : ['conservation','quantile','filter','mask','window','increment','save','processors']}) + cmd_dict['chimera.slayer'] = dict({'required' : ['fasta','reference'], 'optional' : ['name','search','window','increment','match','mismatch','numwanted','parents','minsim','mincov','iters','minbs','minsnp','divergence','realign','split','blastlocation','save','processors']}) cmd_dict['chimera.uchime'] = dict({'required' : ['fasta'], 'optional' : ['name','reference','abskew','chimealns','minh','mindiv','xn','dn','xa','chunks','minchunk','idsmoothwindow','minsmoothid','maxp','skipgaps','skipgaps2','minlen','maxlen','ucl','queryfract','processors']}) cmd_dict['chop.seqs'] = dict({'required' : ['fasta','numbases'], 'optional' : ['countgaps','keep','short']}) cmd_dict['classify.otu'] = dict({'required' : ['list','taxonomy'],'optional' : ['name','cutoff','label','group','probs','basis','reftaxonomy']}) - cmd_dict['classify.seqs'] = dict({'required' : ['fasta','reference','taxonomy'],'optional' : ['name','search','ksize','method','match','mismatch','gapopen','gapextend','numwanted','probs','processors']}) + cmd_dict['classify.seqs'] = dict({'required' : ['fasta','reference','taxonomy'],'optional' : ['name','search','ksize','method','match','mismatch','gapopen','gapextend','numwanted','probs','save','processors']}) + #clear.memory ## not needed in galaxy framework cmd_dict['clearcut'] = dict({'required' : [['phylip','fasta']],'optional' : ['seed','norandom','shuffle','neighbor','expblen','expdist','ntrees','matrixout','kimura','jukes','protein','DNA']}) cmd_dict['cluster'] = dict({'required' : [['phylip','column']] , 'optional' : ['name','method','cutoff','hard','precision','sim','showabund','timing']}) cmd_dict['cluster.classic'] = dict({'required' : ['phylip'] , 'optional' : ['name','method','cutoff','hard','sim','precision']}) @@ -209,7 +210,7 @@ cmd_dict['fastq.info'] = dict({'required' : ['fastq'], 'optional' : []}) cmd_dict['filter.seqs'] = dict({'required' : ['fasta'], 'optional' : ['vertical','trump','soft','hard','processors']}) cmd_dict['get.group'] = dict({'required' : ['shared'], 'optional' : []}) - cmd_dict['get.groups'] = dict({'required' : ['group'], 'optional' : ['groups','accnos','fasta','name','list','taxonomy']}) + cmd_dict['get.groups'] = dict({'required' : ['group'], 'optional' : ['groups','accnos','fasta','name','list','shared','taxonomy']}) cmd_dict['get.lineage'] = dict({'required' : ['taxonomy','taxon'],'optional' : ['fasta','name','group','list','alignreport','dups']}) cmd_dict['get.otulist'] = dict({'required' : ['list'], 'optional' : ['label','sort']}) cmd_dict['get.oturep'] = dict({'required' : ['fasta','list',['phylip','column']], 'optional' : ['name','label','group','groups','sorted','precision','cutoff','large','weighted']}) @@ -223,7 +224,7 @@ cmd_dict['heatmap.bin'] = dict({'required' : [['list', 'sabund', 'rabund', 'shared']], 'optional' : ['label','groups','scale','sorted','numotu','fontsize']}) cmd_dict['heatmap.sim'] = dict({'required' : [['shared','phylip','column']], 'optional' : ['calc','name','label','groups']}) cmd_dict['homova'] = dict({'required' : ['phylip','design'] , 'optional' : ['alpha','iters']}) - cmd_dict['indicator'] = dict({'required' : ['tree',['shared','relabund']], 'optional' : ['groups','label','design']}) + cmd_dict['indicator'] = dict({'required' : [['tree','design'],['shared','relabund']], 'optional' : ['groups','label']}) cmd_dict['libshuff'] = dict({'required' : ['phylip','group'],'optional' : ['groups','iters','form','sim','step','cutoff']}) cmd_dict['list.seqs'] = dict({'required' : [['fasta','name','group','list','alignreport','taxonomy']]}) cmd_dict['make,fastq'] = dict({'required' : ['fasta','qfile'] , 'optional' : []}) @@ -246,7 +247,7 @@ cmd_dict['pre.cluster'] = dict({'required' : ['fasta'], 'optional' : ['name','diffs']}) cmd_dict['rarefaction.shared'] = dict({'required' : ['shared'], 'optional' : ['calc','label','iters','groups','jumble']}) cmd_dict['rarefaction.single'] = dict({'required' : [['list', 'sabund', 'rabund', 'shared']], 'optional' : ['calc','abund','iters','label','freq','processors']}) - cmd_dict['remove.groups'] = dict({'required' : ['group'], 'optional' : ['groups','accnos','fasta','name','list','taxonomy']}) + cmd_dict['remove.groups'] = dict({'required' : ['group'], 'optional' : ['groups','accnos','fasta','name','list','shared','taxonomy']}) cmd_dict['remove.lineage'] = dict({'required' : ['taxonomy','taxon'],'optional' : ['fasta','name','group','list','alignreport','dups']}) cmd_dict['remove.otus'] = dict({'required' : ['group','list','label'], 'optional' : ['groups','accnos']}) cmd_dict['remove.rare'] = dict({'required' : [['list','sabund','rabund','shared'],'nseqs'], 'optional' : ['group','groups','label','bygroup']}) @@ -469,6 +470,7 @@ parser.add_option( '--sabund', dest='sabund', help='' ) parser.add_option( '--shared', dest='shared', help='' ) parser.add_option( '--relabund', dest='relabund', help='' ) + parser.add_option( '--makerelabund', dest='makerelabund', help='Whether to convert to a relabund file' ) parser.add_option( '--ordergroup', dest='ordergroup', help='') # include read.tree options parser.add_option( '--tree', dest='tree', help='' ) diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/otu.hierarchy.xml --- a/mothur/tools/mothur/otu.hierarchy.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/otu.hierarchy.xml Tue Aug 09 11:39:27 2011 -0500 @@ -53,7 +53,7 @@ The otu.hierarchy_ command relates OTUs from a list_ at different distances. -.. _list_file: http://www.mothur.org/wiki/List_file +.. _list: http://www.mothur.org/wiki/List_file .. _otu.hierarchy: http://www.mothur.org/wiki/Otu.hierarchy diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/remove.groups.xml --- a/mothur/tools/mothur/remove.groups.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/remove.groups.xml Tue Aug 09 11:39:27 2011 -0500 @@ -1,4 +1,4 @@ - + Remove groups from groups,fasta,names,list,taxonomy mothur_wrapper.py @@ -20,22 +20,25 @@ #end if #if $fasta_in.__str__ != "None" and len($fasta_in.__str__) > 0: --fasta=$fasta_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($fasta_in.__str__)) + ":'" + $fasta_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($fasta_in.__str__)) + ":'" + $fasta_out.__str__] #end if #if $name_in.__str__ != "None" and len($name_in.__str__) > 0: --name=$name_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($name_in.__str__)) + ":'" + $name_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($name_in.__str__)) + ":'" + $name_out.__str__] #end if #if $list_in.__str__ != "None" and len($list_in.__str__) > 0: --list=$list_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($list_in.__str__)) + ":'" + $list_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($list_in.__str__)) + ":'" + $list_out.__str__] + #end if + #if $shared_in.__str__ != "None" and len($shared_in.__str__) > 0: + --shared=$shared_in + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.*\\.pick\\.\2',$os.path.basename($shared_in.__str__)) + ":'" + $shared_out.__str__] #end if #if $taxonomy_in.__str__ != "None" and len($taxonomy_in.__str__) > 0: --taxonomy=$taxonomy_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.\2',$os.path.basename($taxonomy_in.__str__)) + ":'" + $taxonomy_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.pick\\.\2',$os.path.basename($taxonomy_in.__str__)) + ":'" + $taxonomy_out.__str__] #end if --result=#echo ','.join($results) - @@ -60,6 +63,7 @@ + @@ -75,6 +79,9 @@ list_in != None + + shared_in != None + taxonomy_in != None diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/remove.lineage.xml --- a/mothur/tools/mothur/remove.lineage.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/remove.lineage.xml Tue Aug 09 11:39:27 2011 -0500 @@ -36,11 +36,12 @@ @@ -53,7 +54,6 @@ - @@ -88,7 +88,7 @@ **Command Documenation** -The remove.lineage_ command reads a taxonomy_ file and a taxon and generates a new file that contains only the sequences in the that are not from that taxon. You may also include either a fasta_, name_, group_, list_, or align.report_ file to this command and mothur will generate new files for each of those containing only the selected sequences. +The remove.lineage_ command reads a taxonomy_ file and a taxon and generates a new file that contains only the sequences in the that are not from that taxon. You may also include either a fasta, name_, group_, list_, or align.report_ file to this command and mothur will generate new files for each of those containing only the selected sequences. .. _taxonomy: http://www.mothur.org/wiki/Taxonomy_outline .. _name: http://www.mothur.org/wiki/Name_file diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/remove.otus.xml --- a/mothur/tools/mothur/remove.otus.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/remove.otus.xml Tue Aug 09 11:39:27 2011 -0500 @@ -19,8 +19,8 @@ #end if #end if #set results = ["'^mothur.\S+\.logfile$:'" + $logfile.__str__] - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.' + $label.__str__ + '.\2',$os.path.basename($group_in.__str__)) + ":'" + $group_out.__str__] - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.pick.'+ $label.__str__ + '.\2',$os.path.basename($list_in.__str__)) + ":'" + $list_out.__str__] + #set results = $results + ["'" + $re.sub(r'^(.*)\.(.*?)',r'\1.pick.' + $label.__str__ + '.\2',$os.path.basename($group_in.__str__)) + ":'" + $group_out.__str__] + #set results = $results + ["'" + $re.sub(r'^(.*)\.(.*?)',r'\1.pick.'+ $label.__str__ + '.\2',$os.path.basename($list_in.__str__)) + ":'" + $list_out.__str__] --result=#echo ','.join($results) @@ -38,7 +38,7 @@ - + At least one group must be selected diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/remove.seqs.xml --- a/mothur/tools/mothur/remove.seqs.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/remove.seqs.xml Tue Aug 09 11:39:27 2011 -0500 @@ -41,8 +41,8 @@ - - + + @@ -53,23 +53,12 @@ - - - fasta_in != None - - - - qfile_in != None - - - name_in != None diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/sub.sample.xml --- a/mothur/tools/mothur/sub.sample.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/sub.sample.xml Tue Aug 09 11:39:27 2011 -0500 @@ -1,4 +1,4 @@ - + Create a sub sample mothur_wrapper.py @@ -9,15 +9,15 @@ --outputdir='$logfile.extra_files_path' #if $input.format == "fasta": --fasta=$input.fasta_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.fasta_in.__str__)) + ":'" + $fasta_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.fasta_in.__str__)) + ":'" + $fasta_out.__str__] #if $input.name_in.__str__ != "None" and len($input.name_in.__str__) > 0: --name=$input.name_in - ## #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.name_in.__str__)) + ":'" + $names_out.__str__] + ## #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.name_in.__str__)) + ":'" + $names_out.__str__] #end if #if $input.use_group.to_filter == "yes": #if $input.use_group.group_in.__str__ != "None" and len($input.use_group.group_in.__str__) > 0: --group=$input.use_group.group_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.use_group.group_in.__str__)) + ":'" + $group_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.use_group.group_in.__str__)) + ":'" + $group_out.__str__] #if $input.use_group.groups.__str__ != "None" and len($input.use_group.groups.__str__) > 0: --groups=$input.use_group.groups #end if @@ -26,11 +26,11 @@ #end if #elif $input.format == "list": --list=$input.otu_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $list_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $list_out.__str__] #if $input.use_group.to_filter == "yes": #if $input.use_group.group_in.__str__ != "None" and len($input.use_group.group_in.__str__) > 0: --group=$input.use_group.group_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.use_group.group_in.__str__)) + ":'" + $group_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.use_group.group_in.__str__)) + ":'" + $group_out.__str__] #if $input.use_group.groups.__str__ != "None" and len($input.use_group.groups.__str__) > 0: --groups=$input.use_group.groups #end if @@ -42,7 +42,7 @@ #end if #elif $input.format == "shared": --shared=$input.otu_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $shared_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.*\\.subsample\\.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $shared_out.__str__] #if $input.groups.__str__ != "None" and len($input.groups.__str__) > 0: --groups=$input.groups #end if @@ -51,13 +51,13 @@ #end if #elif $input.format == "sabund": --sabund=$input.otu_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $sabund_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $sabund_out.__str__] #if $input.label.__str__ != "None" and len($input.label.__str__) > 0: --label=$input.label #end if #elif $input.format == "rabund": --rabund=$input.otu_in - #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1.subsample.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $rabund_out.__str__] + #set results = $results + ["'" + $re.sub(r'(^.*)\.(.*?)',r'\1\\.subsample\\.\2',$os.path.basename($input.otu_in.__str__)) + ":'" + $rabund_out.__str__] #if $input.label.__str__ != "None" and len($input.label.__str__) > 0: --label=$input.label #end if @@ -126,7 +126,7 @@ - + @@ -145,7 +145,7 @@ - + @@ -154,7 +154,7 @@ - + @@ -188,7 +188,7 @@ --> - input['use_group']['group_in'] != None + ((input['format'] == 'fasta' or input['format'] == 'list') and input['use_group'] == 'yes') diff -r 11867a45a43d -r 370b3fc4e7d3 mothur/tools/mothur/summary.single.xml --- a/mothur/tools/mothur/summary.single.xml Thu Jul 14 15:21:01 2011 -0500 +++ b/mothur/tools/mothur/summary.single.xml Tue Aug 09 11:39:27 2011 -0500 @@ -81,6 +81,7 @@ The summary.single_ command produce a summary file that has the calculator value for each line in the OTU data and for all possible comparisons between the different groups in the group_ file. This can be useful if you aren't interested in generating collector's or rarefaction curves for your multi-sample data analysis. It would be worth your while, however, to look at the collector's curves for the calculators you are interested in to determine how sensitive the values are to sampling. If the values are not sensitive to sampling, then you can trust the values. Otherwise, you need to keep sampling. For calc parameter choices see: http://www.mothur.org/wiki/Calculators +.. _group: http://www.mothur.org/wiki/Group_file .. _summary.single: http://www.mothur.org/wiki/Summary.single