Mercurial > repos > iuc > mothur_trim_seqs
comparison trim.seqs.xml @ 2:baa3df51e2b7 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mothur commit 3418f23b9768f5aafb86488f5ec1cb97530d4fb3
author | iuc |
---|---|
date | Tue, 20 Mar 2018 21:56:31 -0400 |
parents | 143f949a5def |
children |
comparison
equal
deleted
inserted
replaced
1:de357ac9db36 | 2:baa3df51e2b7 |
---|---|
5 </macros> | 5 </macros> |
6 <expand macro="requirements"/> | 6 <expand macro="requirements"/> |
7 <expand macro="stdio"/> | 7 <expand macro="stdio"/> |
8 <expand macro="version_command"/> | 8 <expand macro="version_command"/> |
9 <command><![CDATA[ | 9 <command><![CDATA[ |
10 @SHELL_OPTIONS@ | 10 @SHELL_OPTIONS@ |
11 | 11 |
12 ## create symlinks to input datasets | 12 ## create symlinks to input datasets |
13 ln -s "$fasta" fasta.dat && | 13 ln -s '$fasta' fasta.dat && |
14 ln -s "$names" names.dat && | 14 ln -s '$names' names.dat && |
15 ln -s "$count" count.dat && | 15 ln -s '$count' count.dat && |
16 #if $oligo.add == "yes": | 16 #if $oligo.add == "yes": |
17 ln -s "$oligo.oligos" oligo.oligos.dat && | 17 ln -s '$oligo.oligos' oligo.oligos.dat && |
18 #end if | 18 #end if |
19 #if $qual.add2 == "yes": | 19 #if $qual.add2 == "yes": |
20 ln -s "$qual.qfile" qual.qfile.dat && | 20 ln -s '$qual.qfile' qual.qfile.dat && |
21 #end if | 21 #end if |
22 | 22 |
23 echo 'trim.seqs( | 23 echo 'trim.seqs( |
24 fasta=fasta.dat, | 24 fasta=fasta.dat, |
25 minlength=$minlength, | 25 minlength=$minlength, |
26 maxlength=$maxlength, | 26 maxlength=$maxlength, |
27 maxambig=$maxambig, | 27 maxambig=$maxambig, |
28 maxhomop=$maxhomop, | 28 maxhomop=$maxhomop, |
29 keepfirst=$keepfirst, | 29 keepfirst=$keepfirst, |
30 removelast=$removelast, | 30 removelast=$removelast, |
31 #if $oligo.add == "yes": | 31 #if $oligo.add == "yes": |
32 oligos=oligo.oligos.dat, | 32 oligos=oligo.oligos.dat, |
33 bdiffs=$oligo.bdiffs, | 33 bdiffs=$oligo.bdiffs, |
34 pdiffs=$oligo.pdiffs, | 34 pdiffs=$oligo.pdiffs, |
35 tdiffs=$oligo.tdiffs, | 35 tdiffs=$oligo.tdiffs, |
36 ldiffs=$oligo.ldiffs, | 36 ldiffs=$oligo.ldiffs, |
37 sdiffs=$oligo.sdiffs, | 37 sdiffs=$oligo.sdiffs, |
38 keepforward=$oligo.keepforward, | 38 keepforward=$oligo.keepforward, |
39 allfiles=$oligo.allfiles, | 39 allfiles=$oligo.allfiles, |
40 #end if | 40 #end if |
41 #if $qual.add2 == "yes": | 41 #if $qual.add2 == "yes": |
42 qfile=qual.qfile.dat, | 42 qfile=qual.qfile.dat, |
43 qaverage=$qual.qaverage, | 43 qaverage=$qual.qaverage, |
44 qthreshold=$qual.qthreshold, | 44 qthreshold=$qual.qthreshold, |
45 qwindowaverage=$qual.qwindowaverage, | 45 qwindowaverage=$qual.qwindowaverage, |
46 qwindowsize=$qual.qwindowsize, | 46 qwindowsize=$qual.qwindowsize, |
47 rollaverage=$qual.rollaverage, | 47 rollaverage=$qual.rollaverage, |
48 qstepsize=$qual.qstepsize, | 48 qstepsize=$qual.qstepsize, |
49 qtrim=$qual.qtrim, | 49 qtrim=$qual.qtrim, |
50 #end if | 50 #end if |
51 flip=$flip, | 51 flip=$flip, |
52 #if $names: | 52 #if $names: |
53 name=names.dat, | 53 name=names.dat, |
54 #end if | 54 #end if |
55 logtransform=$logtransform, | 55 logtransform=$logtransform, |
56 checkorient=$checkorient, | 56 checkorient=$checkorient, |
57 #if $count: | 57 #if $count: |
58 count=count.dat, | 58 count=count.dat, |
59 #end if | 59 #end if |
60 processors='\${GALAXY_SLOTS:-8}' | 60 processors='\${GALAXY_SLOTS:-8}' |
61 )' | 61 )' |
62 | sed 's/ //g' ## mothur trips over whitespace | 62 | sed 's/ //g' ## mothur trips over whitespace |
63 | mothur | 63 | mothur |
64 | tee mothur.out.log | 64 | tee mothur.out.log |
65 ## prevent these two files from being gathered into collection | 65 |
66 && mv fasta.trim.fasta fasta.trim | 66 ## prevent these two files from being gathered into collection |
67 && mv fasta.scrap.fasta fasta.scrap | 67 && mv fasta.trim.fasta fasta.trim |
68 && mv fasta.scrap.fasta fasta.scrap | |
68 ]]></command> | 69 ]]></command> |
69 <inputs> | 70 <inputs> |
70 <param name="fasta" type="data" format="fasta" label="fasta - Sequences"/> | 71 <param name="fasta" type="data" format="fasta" label="fasta - Sequences"/> |
71 <param name="names" type="data" format="mothur.names" optional="true" label="name - Sequence representative name list"/> | 72 <param name="names" type="data" format="mothur.names" optional="true" label="name - Sequence representative name list"/> |
72 <param name="minlength" type="integer" value="0" min="0" label="minlength - Minimum Sequence Length (default 0, ignored if < 1 )"/> | 73 <param name="minlength" type="integer" value="0" min="0" label="minlength - Minimum Sequence Length (default 0, ignored if < 1 )"/> |
111 </conditional> | 112 </conditional> |
112 <param name="flip" type="boolean" truevalue="true" falsevalue="false" checked="false" label="flip - reverse complement the trimmed sequences"/> | 113 <param name="flip" type="boolean" truevalue="true" falsevalue="false" checked="false" label="flip - reverse complement the trimmed sequences"/> |
113 <param name="count" type="data" format="mothur.count_table" optional="true" label="count - a count_table" help="The count file is similar to the name file in that it is used to represent the number of duplicate sequences for a given representative sequence. If you run trim.seqs with an oligos file that contains group labels, trim.seqs will create a new *.trim.count_table with the group information included. "/> | 114 <param name="count" type="data" format="mothur.count_table" optional="true" label="count - a count_table" help="The count file is similar to the name file in that it is used to represent the number of duplicate sequences for a given representative sequence. If you run trim.seqs with an oligos file that contains group labels, trim.seqs will create a new *.trim.count_table with the group information included. "/> |
114 <param name="logtransform" type="boolean" truevalue="true" falsevalue="false" checked="false" label="logtransform" help="allows you to indicate you want the averages for the qwindowaverage, rollaverage and qaverage to be calculated using a logtransform."/> | 115 <param name="logtransform" type="boolean" truevalue="true" falsevalue="false" checked="false" label="logtransform" help="allows you to indicate you want the averages for the qwindowaverage, rollaverage and qaverage to be calculated using a logtransform."/> |
115 <param name="checkorient" type="boolean" truevalue="true" falsevalue="false" checked="false" label="checkorient - search the reverse complement?" help="If you are running the trim.seqs command with paired barcodes or primers, you can use the checkorient parameter. When checkorient=t and mothur can't find the barcodes and primers, it will search the reverse compliment. "/> | 116 <param name="checkorient" type="boolean" truevalue="true" falsevalue="false" checked="false" label="checkorient - search the reverse complement?" help="If you are running the trim.seqs command with paired barcodes or primers, you can use the checkorient parameter. When checkorient=t and mothur can't find the barcodes and primers, it will search the reverse compliment. "/> |
117 <expand macro="param-savelog"/> | |
116 </inputs> | 118 </inputs> |
117 <outputs> | 119 <outputs> |
118 <expand macro="logfile-output"/> | 120 <expand macro="logfile-output"/> |
119 <data name="trim_fasta" format_source="fasta" from_work_dir="fasta.trim" label="${tool.name} on ${on_string}: trim.fasta"/> | 121 <data name="trim_fasta" format_source="fasta" from_work_dir="fasta.trim" label="${tool.name} on ${on_string}: trim.fasta"/> |
120 <data name="scrap_fasta" format_source="fasta" from_work_dir="fasta.scrap" label="${tool.name} on ${on_string}: scrap.fasta"/> | 122 <data name="scrap_fasta" format_source="fasta" from_work_dir="fasta.scrap" label="${tool.name} on ${on_string}: scrap.fasta"/> |
155 <param name="maxhomop" value="4"/> | 157 <param name="maxhomop" value="4"/> |
156 <output name="trim_fasta" md5="14dcaa23735a3f545e7014a69b002859" ftype="fasta"/> | 158 <output name="trim_fasta" md5="14dcaa23735a3f545e7014a69b002859" ftype="fasta"/> |
157 <output name="scrap_fasta" md5="4f791b7684662f1f962970af46429e24" ftype="fasta"/> | 159 <output name="scrap_fasta" md5="4f791b7684662f1f962970af46429e24" ftype="fasta"/> |
158 <output name="trim_names" md5="1b8c6c47052bb69524ef56ebb764fb8f" ftype="mothur.names"/> | 160 <output name="trim_names" md5="1b8c6c47052bb69524ef56ebb764fb8f" ftype="mothur.names"/> |
159 <output name="scrap_names" md5="80f9252837e4b189f06ec00469b88e85" ftype="mothur.names"/> | 161 <output name="scrap_names" md5="80f9252837e4b189f06ec00469b88e85" ftype="mothur.names"/> |
162 <param name="savelog" value="true"/> | |
160 <expand macro="logfile-test"/> | 163 <expand macro="logfile-test"/> |
161 </test> | 164 </test> |
162 <test><!-- test with count table --> | 165 <test><!-- test with count table --> |
163 <param name="fasta" value="amazon.fasta" ftype="fasta"/> | 166 <param name="fasta" value="amazon.fasta" ftype="fasta"/> |
164 <param name="count" value="amazon1.count_table" ftype="mothur.count_table"/> | 167 <param name="count" value="amazon1.count_table" ftype="mothur.count_table"/> |
165 <param name="maxhomop" value="4"/> | 168 <param name="maxhomop" value="4"/> |
166 <output name="trim_fasta" md5="14dcaa23735a3f545e7014a69b002859" ftype="fasta"/> | 169 <output name="trim_fasta" md5="14dcaa23735a3f545e7014a69b002859" ftype="fasta"/> |
167 <output name="scrap_fasta" md5="4f791b7684662f1f962970af46429e24" ftype="fasta"/> | 170 <output name="scrap_fasta" md5="4f791b7684662f1f962970af46429e24" ftype="fasta"/> |
168 <output name="trim_count" md5="836b4d72a8cda3741ef435741783b384" ftype="mothur.count_table"/> | 171 <output name="trim_count" md5="836b4d72a8cda3741ef435741783b384" ftype="mothur.count_table"/> |
169 <output name="scrap_count" md5="04ae9f50c1b6f0d8d7e1ac28f845dd4c" ftype="mothur.count_table"/> | 172 <output name="scrap_count" md5="04ae9f50c1b6f0d8d7e1ac28f845dd4c" ftype="mothur.count_table"/> |
173 <param name="savelog" value="true"/> | |
170 <expand macro="logfile-test"/> | 174 <expand macro="logfile-test"/> |
171 </test> | 175 </test> |
172 <test><!-- test with oligos --> | 176 <test><!-- test with oligos --> |
173 <param name="fasta" value="amazon.fasta" ftype="fasta"/> | 177 <param name="fasta" value="amazon.fasta" ftype="fasta"/> |
174 <param name="add" value="yes"/> | 178 <param name="add" value="yes"/> |
176 <param name="bdiffs" value="100"/> | 180 <param name="bdiffs" value="100"/> |
177 <param name="pdiffs" value="100"/> | 181 <param name="pdiffs" value="100"/> |
178 <output name="trim_fasta" md5="75a8a3ae2d1fe1ff2b860480b84e9bd6" ftype="fasta"/> | 182 <output name="trim_fasta" md5="75a8a3ae2d1fe1ff2b860480b84e9bd6" ftype="fasta"/> |
179 <output name="scrap_fasta" md5="c4fd14e70ab7d1c21d238e87624829d7" ftype="fasta"/> | 183 <output name="scrap_fasta" md5="c4fd14e70ab7d1c21d238e87624829d7" ftype="fasta"/> |
180 <output name="groups_file" md5="198957282c234e825414e175d926046a" ftype="mothur.groups"/> | 184 <output name="groups_file" md5="198957282c234e825414e175d926046a" ftype="mothur.groups"/> |
185 <param name="savelog" value="true"/> | |
181 <expand macro="logfile-test"/> | 186 <expand macro="logfile-test"/> |
182 </test> | 187 </test> |
183 <test><!-- test with oligos and allfiles parameter --> | 188 <test><!-- test with oligos and allfiles parameter --> |
184 <param name="fasta" value="amazon.fasta" ftype="fasta"/> | 189 <param name="fasta" value="amazon.fasta" ftype="fasta"/> |
185 <param name="add" value="yes"/> | 190 <param name="add" value="yes"/> |
193 <element name="F003D144" md5="445124b06d0c9146ae353631794c8093" ftype="mothur.groups"/> | 198 <element name="F003D144" md5="445124b06d0c9146ae353631794c8093" ftype="mothur.groups"/> |
194 </output_collection> | 199 </output_collection> |
195 <output_collection name="fasta_allfiles" count="9"> | 200 <output_collection name="fasta_allfiles" count="9"> |
196 <element name="F003D144" md5="025ff271ac24ecb898863d7fcbfabf10" ftype="fasta"/> | 201 <element name="F003D144" md5="025ff271ac24ecb898863d7fcbfabf10" ftype="fasta"/> |
197 </output_collection> | 202 </output_collection> |
203 <param name="savelog" value="true"/> | |
198 <expand macro="logfile-test"/> | 204 <expand macro="logfile-test"/> |
199 </test> | 205 </test> |
200 <test><!-- test with qfile--> | 206 <test><!-- test with qfile--> |
201 <param name="fasta" value="Fasting_Example1.fasta" ftype="fasta"/> | 207 <param name="fasta" value="Fasting_Example1.fasta" ftype="fasta"/> |
202 <param name="add2" value="yes"/> | 208 <param name="add2" value="yes"/> |
204 <param name="maxhomop" value="4"/> | 210 <param name="maxhomop" value="4"/> |
205 <output name="trim_fasta" md5="d02f74acd6d9fb52b04a93869bb79302" ftype="fasta"/> | 211 <output name="trim_fasta" md5="d02f74acd6d9fb52b04a93869bb79302" ftype="fasta"/> |
206 <output name="scrap_fasta" md5="a4d3ef3d91b4c0146ec84bb7aad3987c" ftype="fasta"/> | 212 <output name="scrap_fasta" md5="a4d3ef3d91b4c0146ec84bb7aad3987c" ftype="fasta"/> |
207 <output name="trim_qual" md5="3d4e2d3c7dd43b90660ab9c923d9eab1" ftype="qual454"/> | 213 <output name="trim_qual" md5="3d4e2d3c7dd43b90660ab9c923d9eab1" ftype="qual454"/> |
208 <output name="scrap_qual" md5="22931236d082c2b77811bbf912c1f4b1" ftype="qual454"/> | 214 <output name="scrap_qual" md5="22931236d082c2b77811bbf912c1f4b1" ftype="qual454"/> |
215 <param name="savelog" value="true"/> | |
209 <expand macro="logfile-test"/> | 216 <expand macro="logfile-test"/> |
210 </test> | 217 </test> |
211 </tests> | 218 </tests> |
212 <help> | 219 <help><![CDATA[ |
213 <![CDATA[ | |
214 | 220 |
215 @MOTHUR_OVERVIEW@ | 221 @MOTHUR_OVERVIEW@ |
216 | 222 |
217 **Command Documentation** | 223 **Command Documentation** |
218 | 224 |
219 The trim.seqs_ command provides the preprocessing features needed to screen and sort pyrosequences. The command will enable you to trim off primer sequences and barcodes, use the barcode information to generate a group file and split a fasta file into sub-files, screen sequences based on the qual file that comes from 454 sequencers, cull sequences based on sequence length and the presence of ambiguous bases and get the reverse complement of your sequences. While this analysis is clearly geared towards pyrosequencing collections, it can also be used with traditional Sanger sequences. | 225 The trim.seqs_ command provides the preprocessing features needed to screen and sort pyrosequences. The command will enable you to trim off primer sequences and barcodes, use the barcode information to generate a group file and split a fasta file into sub-files, screen sequences based on the qual file that comes from 454 sequencers, cull sequences based on sequence length and the presence of ambiguous bases and get the reverse complement of your sequences. While this analysis is clearly geared towards pyrosequencing collections, it can also be used with traditional Sanger sequences. |
220 | 226 |
221 .. _trim.seqs: https://www.mothur.org/wiki/Trim.seqs | 227 .. _trim.seqs: https://www.mothur.org/wiki/Trim.seqs |
222 | 228 |
223 ]]> | 229 ]]></help> |
224 </help> | |
225 <expand macro="citations"/> | 230 <expand macro="citations"/> |
226 </tool> | 231 </tool> |