Mercurial > repos > devteam > samtools_sort
comparison macros.xml @ 9:17bed26ad17e draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tool_collections/samtools/samtools_sort commit 9f6dd28ae31897068c9f8b5d842750d5d7cd600c
author | iuc |
---|---|
date | Wed, 19 Sep 2018 09:53:20 -0400 |
parents | 71d5c34fef4e |
children | c15c63771494 |
comparison
equal
deleted
inserted
replaced
8:71d5c34fef4e | 9:17bed26ad17e |
---|---|
1 <macros> | 1 <macros> |
2 <xml name="requirements"> | 2 <xml name="requirements"> |
3 <requirements> | 3 <requirements> |
4 <requirement type="package" version="1.3.1">samtools</requirement> | 4 <requirement type="package" version="@TOOL_VERSION@">samtools</requirement> |
5 <yield/> | 5 <yield/> |
6 </requirements> | 6 </requirements> |
7 </xml> | 7 </xml> |
8 <token name="@TOOL_VERSION@">1.3.1</token> | 8 <token name="@TOOL_VERSION@">1.9</token> |
9 <token name="@FLAGS@">#set $flags = sum(map(int, str($filter).split(',')))</token> | |
10 <token name="@PREPARE_IDX@"><![CDATA[ | |
11 ##prepare input and indices | |
12 ln -s '$input' infile && | |
13 #if $input.is_of_type('bam'): | |
14 #if str( $input.metadata.bam_index ) != "None": | |
15 ln -s '${input.metadata.bam_index}' infile.bai && | |
16 #else: | |
17 samtools index infile infile.bai && | |
18 #end if | |
19 #elif $input.is_of_type('cram'): | |
20 #if str( $input.metadata.cram_index ) != "None": | |
21 ln -s '${input.metadata.cram_index}' infile.crai && | |
22 #else: | |
23 samtools index infile infile.crai && | |
24 #end if | |
25 #end if | |
26 ]]></token> | |
27 <token name="@PREPARE_IDX_MULTIPLE@"><![CDATA[ | |
28 ##prepare input and indices | |
29 #for $i, $bam in enumerate( $input_bams ): | |
30 ln -s '$bam' '${i}' && | |
31 #if $bam.is_of_type('bam'): | |
32 #if str( $bam.metadata.bam_index ) != "None": | |
33 ln -s '${bam.metadata.bam_index}' '${i}.bai' && | |
34 #else: | |
35 samtools index '${i}' '${i}.bai' && | |
36 #end if | |
37 #elif $bam.is_of_type('cram'): | |
38 #if str( $bam.metadata.cram_index ) != "None": | |
39 ln -s '${bam.metadata.cram_index}' '${i}.crai' && | |
40 #else: | |
41 samtools index '${i}' '${i}.crai' && | |
42 #end if | |
43 #end if | |
44 #end for | |
45 ]]></token> | |
46 <token name="@PREPARE_FASTA_IDX@"><![CDATA[ | |
47 ##checks for reference data ($addref_cond.addref_select=="history" or =="cached") | |
48 ##and sets the -t/-T parameters accordingly: | |
49 ##- in case of history a symbolic link is used because samtools (view) will generate | |
50 ## the index which might not be possible in the directory containing the fasta file | |
51 ##- in case of cached the absolute path is used which allows to read the cram file | |
52 ## without specifying the reference | |
53 #if $addref_cond.addref_select == "history": | |
54 ln -s '${addref_cond.ref}' reference.fa && | |
55 samtools faidx reference.fa && | |
56 #set reffa=str($addref_cond.ref) | |
57 #set reffai="reference.fa.fai" | |
58 #elif $addref_cond.addref_select == "cached": | |
59 #set reffa=str($addref_cond.ref.fields.path) | |
60 #set reffai=str($addref_cond.ref.fields.path) | |
61 #else | |
62 #set reffa=None | |
63 #set reffai=None | |
64 #end if | |
65 ]]></token> | |
66 <token name="@ADDTHREADS@"><![CDATA[ | |
67 ##compute the number of ADDITIONAL threads to be used by samtools (-@) | |
68 addthreads=\${GALAXY_SLOTS:-1} && (( addthreads-- )) && | |
69 ]]></token> | |
70 | |
71 <token name="@ADDMEMORY@"><![CDATA[ | |
72 ##compute the number of memory available to samtools sort (-m) | |
73 ##use only 75% of available: https://github.com/samtools/samtools/issues/831 | |
74 addmemory=\${GALAXY_MEMORY_MB_PER_SLOT:-768} && | |
75 ((addmemory=addmemory*75/100)) && | |
76 ]]></token> | |
77 <xml name="flag_options"> | |
78 <option value="1">read is paired</option> | |
79 <option value="2">read is mapped in a proper pair</option> | |
80 <option value="4">read is unmapped</option> | |
81 <option value="8">mate is unmapped</option> | |
82 <option value="16">read reverse strand</option> | |
83 <option value="32">mate reverse strand</option> | |
84 <option value="64">read is the first in a pair</option> | |
85 <option value="128">read is the second in a pair</option> | |
86 <option value="256">alignment or read is not primary</option> | |
87 <option value="512">read fails platform/vendor quality checks</option> | |
88 <option value="1024">read is a PCR or optical duplicate</option> | |
89 <option value="2048">supplementary alignment</option> | |
90 </xml> | |
9 <xml name="citations"> | 91 <xml name="citations"> |
10 <citations> | 92 <citations> |
11 <citation type="bibtex"> | 93 <citation type="bibtex"> |
12 @misc{SAM_def, | 94 @misc{SAM_def, |
13 title={Definition of SAM/BAM format}, | 95 title={Definition of SAM/BAM format}, |
47 <xml name="stdio"> | 129 <xml name="stdio"> |
48 <stdio> | 130 <stdio> |
49 <exit_code range="1:" level="fatal" description="Error" /> | 131 <exit_code range="1:" level="fatal" description="Error" /> |
50 </stdio> | 132 </stdio> |
51 </xml> | 133 </xml> |
52 <token name="@no-chrom-options@"> | |
53 ----- | |
54 | |
55 .. class:: warningmark | |
56 | |
57 **No options available? How to re-detect metadata** | |
58 | |
59 If you see a "No options available" within the "**Select references (chromosomes and contigs) you would like to restrict bam to**" drop down, you need to re-detect metadata for the dataset you are trying to process. To do this follow these steps: | |
60 | |
61 1. Click on the **pencil** icon adjacent to the dataset in the history | |
62 2. A new menu will appear in the center pane of the interface | |
63 3. Click **Datatype** tab | |
64 4. Set **New Type** to **BAM** | |
65 5. Click **Save** | |
66 | |
67 The medatada will be re-detected and you will be able to see the list of reference sequences in the "**Select references (chromosomes and contigs) you would like to restrict bam to**" drop-down. | |
68 </token> | |
69 </macros> | 134 </macros> |