Mercurial > repos > devteam > bwa
comparison bwa.xml @ 26:2477830927ec draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bwa commit 6e9628b1d92fdb358b79959ad54a456cfa46fa33
author | iuc |
---|---|
date | Fri, 17 May 2024 21:09:07 +0000 |
parents | e188dc7a68e6 |
children |
comparison
equal
deleted
inserted
replaced
25:e188dc7a68e6 | 26:2477830927ec |
---|---|
1 <?xml version="1.0"?> | 1 <?xml version="1.0"?> |
2 <tool id="bwa" name="Map with BWA" version="@VERSION@.5"> | 2 <tool id="bwa" name="Map with BWA" version="@TOOL_VERSION@" profile="22.05"> |
3 <description>- map short reads (< 100 bp) against reference genome</description> | 3 <description>- map short reads (< 100 bp) against reference genome</description> |
4 <xrefs> | |
5 <xref type="bio.tools">bwa</xref> | |
6 </xrefs> | |
7 <macros> | 4 <macros> |
8 <import>read_group_macros.xml</import> | 5 <import>read_group_macros.xml</import> |
9 <import>bwa_macros.xml</import> | 6 <import>bwa_macros.xml</import> |
10 <token name="@command_options@"> | 7 <token name="@command_options@"> |
11 #if str( $analysis_type.analysis_type_selector ) == "full": | 8 #if str( $analysis_type.analysis_type_selector ) == "full": |
75 <when value="do_not_set"> | 72 <when value="do_not_set"> |
76 <!-- do nothing --> | 73 <!-- do nothing --> |
77 </when> | 74 </when> |
78 </xml> | 75 </xml> |
79 </macros> | 76 </macros> |
77 <expand macro="bio_tools"/> | |
80 <expand macro="requirements"/> | 78 <expand macro="requirements"/> |
81 <expand macro="stdio"/> | 79 <expand macro="stdio"/> |
82 <command> | 80 <command> |
83 <![CDATA[ | 81 <![CDATA[ |
84 @pipefail@ | 82 @pipefail@ |
312 <data format="bam" name="bam_output" label="${tool.name} on ${on_string} (mapped reads in BAM format)"> | 310 <data format="bam" name="bam_output" label="${tool.name} on ${on_string} (mapped reads in BAM format)"> |
313 <expand macro="dbKeyActionsBwaMem"/> | 311 <expand macro="dbKeyActionsBwaMem"/> |
314 </data> | 312 </data> |
315 </outputs> | 313 </outputs> |
316 <tests> | 314 <tests> |
317 <test> | 315 <!-- `samtools sort` in the new update adds PG lines to the output so the lines_diff is changed from "2" to "4" --> |
316 <test expect_num_outputs="1"> | |
318 <param name="reference_source_selector" value="history"/> | 317 <param name="reference_source_selector" value="history"/> |
319 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> | 318 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> |
320 <param name="input_type_selector" value="single"/> | 319 <param name="input_type_selector" value="single"/> |
321 <param name="fastq_input1" ftype="fasta" value="bwa-mem-fasta1.fa"/> | 320 <param name="fastq_input1" ftype="fasta" value="bwa-mem-fasta1.fa"/> |
322 <param name="analysis_type_selector" value="illumina"/> | 321 <param name="analysis_type_selector" value="illumina"/> |
323 <output name="bam_output" ftype="bam" file="bwa-aln-test1-fasta.bam" lines_diff="2"/> | 322 <output name="bam_output" ftype="bam" file="bwa-aln-test1-fasta.bam" lines_diff="4"/> |
324 </test> | 323 </test> |
325 <test> | 324 <test expect_num_outputs="1"> |
326 <param name="reference_source_selector" value="history"/> | 325 <param name="reference_source_selector" value="history"/> |
327 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> | 326 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> |
328 <param name="input_type_selector" value="paired"/> | 327 <param name="input_type_selector" value="paired"/> |
329 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/> | 328 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/> |
330 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/> | 329 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/> |
331 <param name="analysis_type_selector" value="illumina"/> | 330 <param name="analysis_type_selector" value="illumina"/> |
332 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="2"/> | 331 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="4"/> |
333 </test> | 332 </test> |
334 <test> | 333 <test expect_num_outputs="1"> |
335 <param name="reference_source_selector" value="history"/> | 334 <param name="reference_source_selector" value="history"/> |
336 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> | 335 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> |
337 <param name="input_type_selector" value="paired"/> | 336 <param name="input_type_selector" value="paired"/> |
338 <param name="fastq_input1" ftype="fastqsanger.gz" value="bwa-mem-fastq1.fq.gz"/> | 337 <param name="fastq_input1" ftype="fastqsanger.gz" value="bwa-mem-fastq1.fq.gz"/> |
339 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/> | 338 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/> |
340 <param name="analysis_type_selector" value="illumina"/> | 339 <param name="analysis_type_selector" value="illumina"/> |
341 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="2"/> | 340 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="4"/> |
342 </test> | 341 </test> |
343 <test> | 342 <test expect_num_outputs="1"> |
344 <param name="reference_source_selector" value="history"/> | 343 <param name="reference_source_selector" value="history"/> |
345 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> | 344 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> |
346 <param name="input_type_selector" value="paired_bam"/> | 345 <param name="input_type_selector" value="paired_bam"/> |
347 <param name="bam_input" ftype="bam" value="bwa-aln-bam-input.bam"/> | 346 <param name="bam_input" ftype="bam" value="bwa-aln-bam-input.bam"/> |
348 <param name="analysis_type_selector" value="illumina"/> | 347 <param name="analysis_type_selector" value="illumina"/> |
349 <output name="bam_output" ftype="bam" file="bwa-aln-test2.bam" lines_diff="2"/> | 348 <output name="bam_output" ftype="bam" file="bwa-aln-test2.bam" lines_diff="4"/> |
350 </test> | 349 </test> |
351 <test> | 350 <test expect_num_outputs="1"> |
352 <param name="reference_source_selector" value="history"/> | 351 <param name="reference_source_selector" value="history"/> |
353 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> | 352 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/> |
354 <param name="input_type_selector" value="paired"/> | 353 <param name="input_type_selector" value="paired"/> |
355 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/> | 354 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/> |
356 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/> | 355 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/> |
357 <param name="rg_selector" value="set"/> | 356 <param name="rg_selector" value="set"/> |
358 <param name="ID" value="rg1"/> | 357 <param name="ID" value="rg1"/> |
359 <param name="PL" value="CAPILLARY"/> | 358 <param name="PL" value="CAPILLARY"/> |
360 <param name="analysis_type_selector" value="illumina"/> | 359 <param name="analysis_type_selector" value="illumina"/> |
361 <output name="bam_output" ftype="bam" file="bwa-aln-test3.bam" lines_diff="2"/> | 360 <output name="bam_output" ftype="bam" file="bwa-aln-test3.bam" lines_diff="5"/> |
362 </test> | 361 </test> |
363 </tests> | 362 </tests> |
364 <help><![CDATA[ | 363 <help><![CDATA[ |
365 **What is does** | 364 **What it does** |
366 | 365 |
367 BWA is a software package for mapping low-divergent sequences against a large reference genome, such as the | 366 BWA_ is a software package for mapping low-divergent sequences against a large reference genome, such as the human genome. |
368 human genome. The bwa-aln algorithm is designed for Illumina sequence reads up to 100bp. For longer reads use | 367 The bwa-aln algorithm is designed for Illumina sequence reads up to 100bp. For longer reads use the separate BWA-MEM Galaxy tool. |
369 the separate BWA-MEM Galaxy tool. | 368 |
370 | 369 This Galaxy tool wraps the bwa-aln, bwa-samse and -sampe modules of the BWA read mapping tool: |
371 This Galaxy tool wraps bwa-aln, bwa-samse and -sampe modules of bwa read mapping tool: | |
372 | 370 |
373 - **bwa aln** - actual mapper placing reads onto the reference sequence | 371 - **bwa aln** - actual mapper placing reads onto the reference sequence |
374 - **bwa samse** - post-processor converting suffix array coordinates into genome coordinates in SAM format for | 372 - **bwa samse** - post-processor converting suffix array coordinates into genome coordinates in SAM format for |
375 single reads | 373 single reads |
376 - **bam sampe** - post-processor for paired reads | 374 - **bam sampe** - post-processor for paired reads |
377 | 375 |
376 For more details about the different modules of the BWA package see the `BWA manual`_. | |
378 | 377 |
379 The Galaxy implementation takes fastq or BAM (unaligned BAM) datasets as input and produces output in BAM format, | 378 The Galaxy implementation takes fastq or BAM (unaligned BAM) datasets as input and produces output in BAM format, |
380 which can be further processed using various BAM utilities exiting in Galaxy (BAMTools, SAMTools, Picard). | 379 which can be further processed using various BAM utilities existing in Galaxy (BAMTools, SAMTools, Picard). |
381 | 380 |
382 ----- | 381 ----- |
383 | 382 |
384 **Indices: Selecting reference genomes for BWA** | 383 @ref_genomes@ |
385 | |
386 The Galaxy wrapper for BWA allows you to select between precomputed and user-defined indices for reference genomes | |
387 using the **Will you select a reference genome from your history or use a built-in index?** select box. | |
388 | |
389 This select box has two options: | |
390 | |
391 1. **Use a built-in genome index** - when selected (this is default), Galaxy provides the user with **Select | |
392 reference genome index** dropdown. Genomes listed in this dropdown have been pre-indexed with bwa index utility | |
393 and are ready to be mapped against. | |
394 2. **Use a genome from the history and build index** - when selected, Galaxy provides the user with **Select | |
395 reference genome sequence** dropdown. This dropdown is populated by all FASTA formatted files listed in your | |
396 current history. If your genome of interest is uploaded into history it will be shown there. Selecting a genome | |
397 from this dropdown will cause Galaxy to first transparently index it using `bwa index` command, and then run | |
398 mapping with `bwa aln`. | |
399 | |
400 | |
401 If your genome of interest is not listed here you have two choices: | |
402 | |
403 1. Contact galaxy team using **Help->Support** link at the top of the interface and let us know that an index | |
404 needs to be added | |
405 2. Upload your genome of interest as a FASTA file to Galaxy history and selected **Use a genome from the history | |
406 and build index** option. | |
407 | |
408 | 384 |
409 @RG@ | 385 @RG@ |
410 | 386 |
411 @info@ | 387 @links@ |
412 ]]></help> | 388 ]]></help> |
413 <citations> | 389 <citations> |
414 <citation type="doi">10.1093/bioinformatics/btp324</citation> | 390 <citation type="doi">10.1093/bioinformatics/btp324</citation> |
415 <citation type="doi">10.1093/bioinformatics/btp698</citation> | 391 <citation type="doi">10.1093/bioinformatics/btp698</citation> |
416 </citations> | 392 </citations> |