Mercurial > repos > jtilman > flexbar
annotate flexbar.xml @ 32:c554e6d251fd draft
Uploaded
| author | jtilman |
|---|---|
| date | Tue, 19 Jun 2018 10:00:20 -0400 |
| parents | d43903772a88 |
| children | 4a095317caca |
| rev | line source |
|---|---|
| 19 | 1 |
| 31 | 2 <!-- Flexbar tool definition for Galaxy, version 3.4 --> |
| 19 | 3 <!-- Author: Johannes Roehr --> |
| 4 | |
| 5 | |
| 31 | 6 <tool id="flexbar" name="Flexbar" version="3.4" force_history_refresh="True"> |
| 22 | 7 |
| 19 | 8 <description>flexible barcode and adapter removal</description> |
| 9 | |
| 10 <requirements> | |
| 31 | 11 <requirement type="binary" version="3.4">flexbar</requirement> |
| 19 | 12 </requirements> |
| 13 | |
| 31 | 14 <version_command>flexbar --versions</version_command> |
| 19 | 15 |
| 16 <command interpreter="perl"> | |
| 17 | |
| 18 flexbar.pl flexbar | |
| 19 | |
|
23
0ef7736004bb
Use of galaxy variable to control number of threads.
jtilman
parents:
22
diff
changeset
|
20 --threads \${GALAXY_SLOTS:-1} |
| 19 | 21 |
| 22 --reads $reads | |
| 23 | |
| 24 #if $cReads2.select == "on": | |
| 25 #if $cReads2.reads2.ext == $reads.ext: | |
| 26 --reads2 $cReads2.reads2 | |
| 27 #end if | |
| 28 #end if | |
| 29 | |
| 30 --max-uncalled $maxUncalled | |
| 22 | 31 --min-read-length $minReadLen |
| 19 | 32 |
| 33 #if $trimEnds.select == "on": | |
| 34 --pre-trim-left $trimEnds.trimLeft | |
| 35 --pre-trim-right $trimEnds.trimRight | |
| 36 #end if | |
| 37 | |
| 38 #if $cTrimPhred.select == "on": | |
| 31 | 39 --qtrim TAIL --qtrim-threshold $cTrimPhred.trimPhred |
| 40 | |
| 41 #if $reads.ext == "fastqsanger": | |
| 42 --qtrim-format sanger | |
| 43 #end if | |
| 44 #if $reads.ext == "fastqsolexa": | |
| 45 --qtrim-format solexa | |
| 46 #end if | |
| 47 #if $reads.ext == "fastqillumina": | |
| 48 --qtrim-format i1.3 | |
| 49 #end if | |
| 19 | 50 #end if |
| 51 | |
| 52 #if $cTrimLen.select == "on": | |
| 53 --post-trim-length $cTrimLen.trimLen | |
| 54 #end if | |
| 55 | |
| 56 | |
| 57 #if $cBarcodes.select == "on": | |
| 58 --barcodes $cBarcodes.barcodes | |
| 59 | |
| 60 #if $cBarcodes.cbReads.select == "yes": | |
| 61 --barcode-reads $cBarcodes.cbReads.bReads | |
| 62 #end if | |
| 63 | |
| 64 #if $cBarcodes.cbReads.select == "no": | |
| 65 $cBarcodes.cbReads.bKeep | |
| 66 #end if | |
| 67 | |
| 22 | 68 $cBarcodes.bUnassigned |
| 69 | |
| 19 | 70 --barcode-trim-end $cBarcodes.bTrimEnd |
| 71 | |
| 72 #if $cBarcodes.cbTailLen.select == "yes": | |
| 73 --barcode-tail-length $cBarcodes.cbTailLen.bTailLen | |
| 74 #end if | |
| 75 | |
| 76 #if $cBarcodes.cbMinOverlap.select == "yes": | |
| 77 --barcode-min-overlap $cBarcodes.cbMinOverlap.bMinOverlap | |
| 78 #end if | |
| 79 | |
| 31 | 80 --barcode-error-rate $cBarcodes.bThresh |
| 19 | 81 |
| 82 #if $cBarcodes.cbAlignScores.select == "yes": | |
| 83 --barcode-match $bMatch | |
| 84 --barcode-mismatch $bMismatch | |
| 85 --barcode-gap $bGap | |
| 86 #end if | |
| 87 #end if | |
| 88 | |
| 89 | |
| 90 #if $cAdapters.select == "on": | |
| 91 | |
| 92 #if $cAdapters.ccAdapters.select == "data": | |
| 93 --adapters $cAdapters.ccAdapters.adaptersData | |
| 94 #end if | |
| 95 | |
| 96 #if $cAdapters.ccAdapters.select == "seq": | |
| 97 --adapter-seq $cAdapters.ccAdapters.adapterSeq | |
| 98 #end if | |
| 99 | |
| 100 --adapter-trim-end $cAdapters.aTrimEnd | |
| 101 | |
| 102 #if $cAdapters.caTailLen.select == "yes": | |
| 103 --adapter-tail-length $cAdapters.caTailLen.aTailLen | |
| 104 #end if | |
| 105 | |
| 22 | 106 $cAdapters.aReadSet |
| 107 | |
| 19 | 108 --adapter-min-overlap $cAdapters.aMinOverlap |
| 31 | 109 --adapter-error-rate $cAdapters.aThresh |
| 19 | 110 |
| 111 #if $cAdapters.caAlignScores.select == "yes": | |
| 112 --adapter-match $aMatch | |
| 113 --adapter-mismatch $aMismatch | |
| 114 --adapter-gap $aGap | |
| 115 #end if | |
| 116 #end if | |
| 117 | |
| 118 | |
| 22 | 119 #if $cOutput.select == "show": |
| 120 $cOutput.fastaOutput | |
| 121 $cOutput.lenDist | |
| 122 $cOutput.singleReads | |
| 123 #end if | |
| 124 | |
| 19 | 125 #if $cLogging.select == "show": |
| 126 $cLogging.logLevel | |
| 127 $cLogging.numTags | |
| 128 $cLogging.remTags | |
| 129 $cLogging.rndTags | |
| 130 #end if | |
| 131 | |
| 132 | |
| 133 $output $output.id $__new_file_path__ $reads.ext | |
| 134 | |
| 135 </command> | |
| 136 | |
| 22 | 137 |
| 19 | 138 <inputs> |
| 139 | |
| 31 | 140 <param format="fastq,fastqsanger,fastqsolexa,fastqillumina" name="reads" type="data" label="Reads fastq file" optional="false"/> |
| 19 | 141 |
| 142 | |
| 143 <conditional name="cReads2"> | |
| 31 | 144 <param name="select" type="select" label="Paired reads file"> |
| 19 | 145 <option value="off" selected="true">Off</option> |
| 146 <option value="on">On</option> | |
| 147 </param> | |
| 20 | 148 <when value="off"> |
| 149 </when> | |
| 19 | 150 <when value="on"> |
| 31 | 151 <param format="fastq,fastqsanger,fastqsolexa,fastqillumina" name="reads2" type="data" label="Reads 2" optional="false" help="same format as first read set"/> |
| 19 | 152 </when> |
| 153 </conditional> | |
| 154 | |
| 155 | |
| 156 <param name="maxUncalled" size="4" type="integer" value="0" label="1) Max uncalled" optional="false" help="allowed uncalled bases per read"/> | |
| 157 | |
| 158 <conditional name="trimEnds"> | |
| 159 <param name="select" type="select" label="2) Trimming of ends"> | |
| 160 <option value="off" selected="true">Off</option> | |
| 161 <option value="on">On</option> | |
| 162 </param> | |
| 20 | 163 <when value="off"> |
| 164 </when> | |
| 19 | 165 <when value="on"> |
| 166 <param name="trimLeft" size="4" type="integer" value="0" label="Left" optional="false"/> | |
| 167 <param name="trimRight" size="4" type="integer" value="0" label="Right" optional="false" help="trims specified number of bases from read ends"/> | |
| 168 </when> | |
| 169 </conditional> | |
| 170 | |
| 171 <conditional name="cTrimPhred"> | |
| 31 | 172 <param name="select" type="select" label="3) Quality-based trimming"> |
| 19 | 173 <option value="off" selected="true">Off</option> |
| 174 <option value="on">On</option> | |
| 175 </param> | |
| 20 | 176 <when value="off"> |
| 177 </when> | |
| 19 | 178 <when value="on"> |
| 179 <param name="trimPhred" size="4" type="integer" value="10" label="Threshold" optional="false" help="trim right end until specified or higher quality reached"/> | |
| 180 </when> | |
| 181 </conditional> | |
| 182 | |
| 183 | |
| 184 <conditional name="cBarcodes"> | |
| 185 <param name="select" type="select" label="4) Barcode detection"> | |
| 186 <option value="off" selected="true">Off</option> | |
| 187 <option value="on">On</option> | |
| 188 </param> | |
| 20 | 189 <when value="off"> |
| 190 </when> | |
| 19 | 191 <when value="on"> |
| 192 <param format="fasta" name="barcodes" type="data" label="Barcodes" optional="false"/> | |
| 22 | 193 |
| 19 | 194 <conditional name="cbReads"> |
| 195 <param name="select" type="select" label="Separate barcode reads"> | |
| 196 <option value="no" selected="true">No</option> | |
| 197 <option value="yes">Yes</option> | |
| 198 </param> | |
| 199 <when value="yes"> | |
| 31 | 200 <param format="fastq,fastqsanger,fastqsolexa,fastqillumina" name="bReads" type="data" label="Separate barcode reads" optional="false"/> |
| 19 | 201 </when> |
| 202 <when value="no"> | |
| 203 <param name="bKeep" type="select" label="Remove barcodes within reads"> | |
| 204 <option value="" selected="true">Yes</option> | |
| 205 <option value="--barcode-keep">No</option> | |
| 206 </param> | |
| 207 </when> | |
| 208 </conditional> | |
| 22 | 209 |
| 210 <param name="bUnassigned" type="select" label="Include unassigned reads"> | |
| 211 <option value="" selected="true">No</option> | |
| 212 <option value="--barcode-unassigned">Yes</option> | |
| 213 </param> | |
| 214 | |
| 19 | 215 <param name="bTrimEnd" type="select" label="Trim-end mode" optional="false"> |
| 216 <option value="ANY" selected="true">Any</option> | |
| 217 <option value="RIGHT">Right</option> | |
| 31 | 218 <option value="RTAIL">Right tail</option> |
| 19 | 219 <option value="LEFT">Left</option> |
| 31 | 220 <option value="LTAIL">Left tail</option> |
| 19 | 221 </param> |
| 22 | 222 |
| 19 | 223 <conditional name="cbTailLen"> |
| 224 <param name="select" type="select" label="Change tail length"> | |
| 225 <option value="no" selected="true">No</option> | |
| 226 <option value="yes">Yes</option> | |
| 227 </param> | |
| 20 | 228 <when value="no"> |
| 229 </when> | |
| 19 | 230 <when value="yes"> |
| 231 <param name="bTailLen" size="4" type="integer" value="10" label="Tail length" optional="false"/> | |
| 232 </when> | |
| 233 </conditional> | |
| 234 | |
| 235 <conditional name="cbMinOverlap"> | |
| 236 <param name="select" type="select" label="Change min-overlap" help="default: barcode length"> | |
| 237 <option value="no" selected="true">No</option> | |
| 238 <option value="yes">Yes</option> | |
| 239 </param> | |
| 20 | 240 <when value="no"> |
| 241 </when> | |
| 19 | 242 <when value="yes"> |
| 243 <param name="bMinOverlap" size="4" type="integer" value="8" label="Min-overlap" optional="false"/> | |
| 244 </when> | |
| 245 </conditional> | |
| 246 | |
| 247 <param name="bThresh" size="4" type="integer" value="1" label="Threshold" optional="false" help="allowed mismatches and indels per 10 bases"/> | |
| 248 | |
| 249 <conditional name="cbAlignScores"> | |
| 250 <param name="select" type="select" label="Modify alignment scores"> | |
| 251 <option value="no" selected="true">No</option> | |
| 252 <option value="yes">Yes</option> | |
| 253 </param> | |
| 20 | 254 <when value="no"> |
| 255 </when> | |
| 19 | 256 <when value="yes"> |
| 257 <param name="bMatch" size="3" type="integer" value="1" label="Match" optional="false"/> | |
| 258 <param name="bMismatch" size="3" type="integer" value="-1" label="Mismatch" optional="false"/> | |
| 22 | 259 <param name="bGap" size="3" type="integer" value="-9" label="Gap" optional="false"/> |
| 19 | 260 </when> |
| 261 </conditional> | |
| 262 </when> | |
| 263 </conditional> | |
| 264 | |
| 265 | |
| 266 <conditional name="cAdapters"> | |
| 267 <param name="select" type="select" label="5) Adapter removal"> | |
| 268 <option value="off" selected="true">Off</option> | |
| 269 <option value="on">On</option> | |
| 270 </param> | |
| 271 | |
| 20 | 272 <when value="off"> |
| 273 </when> | |
| 19 | 274 <when value="on"> |
| 275 <conditional name="ccAdapters"> | |
| 276 <param name="select" type="select" label="Adapter source"> | |
| 277 <option value="data" selected="true">Fasta</option> | |
| 278 <option value="seq">Sequence</option> | |
| 279 <!-- <option value="file">File</option> --> | |
| 280 </param> | |
| 281 <when value="data"> | |
| 282 <param format="fasta" name="adaptersData" type="data" label="Adapters" optional="false"/> | |
| 283 </when> | |
| 284 <when value="seq"> | |
| 285 <param name="adapterSeq" size="40" label="Adapter" type="text" value="AAAAAAAAAAAAAA" optional="false"/> | |
| 286 </when> | |
| 287 <!-- <when value="file"> | |
| 288 <param name="adaptersFile" type="file" label="Adapters file" optional="false"/> | |
| 289 </when> --> | |
| 290 </conditional> | |
| 291 | |
| 292 <param name="aTrimEnd" type="select" label="Trim-end mode" optional="false"> | |
| 293 <option value="ANY">Any</option> | |
| 294 <option value="RIGHT" selected="true">Right</option> | |
| 31 | 295 <option value="RTAIL">Right tail</option> |
| 19 | 296 <option value="LEFT">Left</option> |
| 31 | 297 <option value="LTAIL">Left tail</option> |
| 19 | 298 </param> |
| 299 | |
| 300 <conditional name="caTailLen"> | |
| 301 <param name="select" type="select" label="Change tail length"> | |
| 302 <option value="no" selected="true">No</option> | |
| 303 <option value="yes">Yes</option> | |
| 304 </param> | |
| 20 | 305 <when value="no"> |
| 306 </when> | |
| 19 | 307 <when value="yes"> |
| 308 <param name="aTailLen" size="4" type="integer" value="10" label="Tail length" optional="false"/> | |
| 309 </when> | |
| 310 </conditional> | |
| 311 | |
| 22 | 312 <param name="aReadSet" type="select" label="Removal for single read set"> |
| 313 <option value="" selected="true">No</option> | |
| 314 <option value="--adapter-read-set 1">1st</option> | |
| 315 <option value="--adapter-read-set 2">2nd</option> | |
| 316 </param> | |
| 317 | |
| 19 | 318 <param name="aMinOverlap" size="4" type="integer" value="1" label="Min-overlap" optional="false"/> |
| 319 <param name="aThresh" size="4" type="integer" value="3" label="Threshold" optional="false" help="allowed mismatches and indels per 10 bases"/> | |
| 320 | |
| 321 <conditional name="caAlignScores"> | |
| 322 <param name="select" type="select" label="Modify alignment scores"> | |
| 323 <option value="no" selected="true">No</option> | |
| 324 <option value="yes">Yes</option> | |
| 325 </param> | |
| 20 | 326 <when value="no"> |
| 327 </when> | |
| 19 | 328 <when value="yes"> |
| 329 <param name="aMatch" size="3" type="integer" value="1" label="Match" optional="false"/> | |
| 330 <param name="aMismatch" size="3" type="integer" value="-1" label="Mismatch" optional="false"/> | |
| 31 | 331 <param name="aGap" size="3" type="integer" value="-6" label="Gap" optional="false"/> |
| 19 | 332 </when> |
| 333 </conditional> | |
| 334 </when> | |
| 335 </conditional> | |
| 336 | |
| 337 | |
| 338 <conditional name="cTrimLen"> | |
| 339 <param name="select" type="select" label="6) Trimming to length"> | |
| 340 <option value="off" selected="true">Off</option> | |
| 341 <option value="on">On</option> | |
| 342 </param> | |
| 20 | 343 <when value="off"> |
| 344 </when> | |
| 19 | 345 <when value="on"> |
| 346 <param name="trimLen" size="4" type="integer" value="30" label="Length" optional="false" help="trim reads to certain length from right"/> | |
| 347 </when> | |
| 348 </conditional> | |
| 349 | |
| 350 <param name="minReadLen" size="4" type="integer" value="18" label="7) Minimum read length" optional="false" help="shorter reads are discarded"/> | |
| 351 | |
| 22 | 352 <conditional name="cOutput"> |
| 353 <param name="select" type="select" label="Output selection"> | |
| 354 <option value="off" selected="true">Off</option> | |
| 355 <option value="show">Show</option> | |
| 356 </param> | |
| 357 <when value="off"> | |
| 358 </when> | |
| 359 <when value="show"> | |
| 360 <param name="fastaOutput" type="select" label="Fasta output"> | |
| 361 <option value="" selected="true">Off</option> | |
| 362 <option value="--fasta-output">Always</option> | |
| 363 </param> | |
| 364 | |
| 365 <param name="lenDist" type="select" label="Read length distribution"> | |
| 366 <option value="" selected="true">Off</option> | |
| 367 <option value="--length-dist">On</option> | |
| 368 </param> | |
| 369 | |
| 370 <param name="singleReads" type="select" label="Single reads"> | |
| 371 <option value="" selected="true">Off</option> | |
| 372 <option value="--single-reads">On</option> | |
| 373 </param> | |
| 374 </when> | |
| 375 </conditional> | |
| 19 | 376 |
| 377 <conditional name="cLogging"> | |
| 378 <param name="select" type="select" label="Logging and tagging options"> | |
| 379 <option value="off" selected="true">Off</option> | |
| 380 <option value="show">Show</option> | |
| 381 </param> | |
| 20 | 382 <when value="off"> |
| 383 </when> | |
| 19 | 384 <when value="show"> |
| 385 <param name="logLevel" type="select" label="Alignment logging"> | |
| 386 <option value="" selected="true">Off</option> | |
| 31 | 387 <option value="--align-log ALL">All</option> |
| 388 <option value="--align-log MOD">Modified</option> | |
| 389 <option value="--align-log TAB">Tabular</option> | |
| 19 | 390 </param> |
| 391 | |
| 392 <param name="numTags" type="select" label="Number tags"> | |
| 393 <option value="" selected="true">Off</option> | |
| 394 <option value="--number-tags">On</option> | |
| 395 </param> | |
| 396 | |
| 397 <param name="remTags" type="select" label="Removal tags"> | |
| 398 <option value="" selected="true">Off</option> | |
| 399 <option value="--removal-tags">On</option> | |
| 400 </param> | |
| 401 | |
| 31 | 402 <param name="rndTags" type="select" label="UMI tags"> |
| 19 | 403 <option value="" selected="true">Off</option> |
| 31 | 404 <option value="--umi-tags">On</option> |
| 19 | 405 </param> |
| 406 </when> | |
| 407 </conditional> | |
| 408 | |
| 409 </inputs> | |
|
23
0ef7736004bb
Use of galaxy variable to control number of threads.
jtilman
parents:
22
diff
changeset
|
410 |
| 19 | 411 <stdio> |
| 412 <exit_code range="1:" level="fatal" description="Error!" /> | |
| 413 </stdio> | |
| 414 | |
| 415 <outputs> | |
| 416 <data format="txt" name="output" metadata_source="reads"/> | |
| 417 </outputs> | |
| 418 | |
| 419 | |
| 420 <help> | |
| 421 | |
| 422 **Description** | |
| 423 | |
| 31 | 424 The program Flexbar preprocesses high-throughput sequencing data efficiently. It demultiplexes barcoded runs and removes adapter sequences. Moreover, trimming and filtering features are provided. Flexbar increases read mapping rates and improves genome as well as transcriptome assemblies. Unique molecular identifiers can be extracted in a flexible way. The program supports sequencing data in fastq format, e.g. from the Illumina platform. Flexbar is available on the project_ page. |
| 19 | 425 |
| 24 | 426 .. _project: https://github.com/seqan/flexbar |
| 19 | 427 |
| 428 ------ | |
| 429 | |
| 430 **Trim-end modes** | |
| 431 | |
| 22 | 432 **Any:** longer side of read remains after overlap removal |
| 19 | 433 |
| 22 | 434 **Left:** right side remains after removal, align before or at read end |
| 19 | 435 |
| 22 | 436 **Right:** left part remains after removal, align after or at read start |
| 19 | 437 |
| 22 | 438 **Left tail:** consider first n bases of reads in alignment |
| 19 | 439 |
| 22 | 440 **Right tail:** use only last n bases, see tail-length options |
| 19 | 441 |
| 442 ------ | |
| 443 | |
| 444 **Documentation** | |
| 445 | |
| 31 | 446 Further documentation is available on the `manual`__ page and via the command line help screen. |
| 19 | 447 |
| 24 | 448 .. __: https://github.com/seqan/flexbar/wiki |
| 19 | 449 |
| 450 ------ | |
| 451 | |
| 31 | 452 **References** |
| 19 | 453 |
| 31 | 454 Johannes T. Roehr, Christoph Dieterich, Knut Reinert: |
| 455 Flexbar 3.0 – SIMD and multicore parallelization. Bioinformatics 2017. | |
| 456 | |
| 457 Matthias Dodt, Johannes T. Roehr, Rina Ahmed, Christoph Dieterich: | |
| 458 Flexbar – flexible barcode and adapter processing for next-generation sequencing platforms. Biology 2012. | |
| 19 | 459 |
| 460 </help> | |
| 461 | |
| 462 </tool> | |
| 463 |