annotate translate_bed_sequences.xml @ 5:c626a939eef7 draft default tip

Uploaded
author jjohnson
date Tue, 12 Jan 2016 14:38:03 -0500
parents 3b526a780849
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
1 <?xml version="1.0"?>
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
2 <tool id="translate_bed_sequences" name="Translate BED Sequences" version="0.1.0">
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
3 <description>3 frame translation of BED augmented with a sequence column</description>
1
639ee511d552 Add biopython requirement
Jim Johnson <jj@umn.edu>
parents: 0
diff changeset
4 <requirements>
639ee511d552 Add biopython requirement
Jim Johnson <jj@umn.edu>
parents: 0
diff changeset
5 <requirement type="package" version="1.62">biopython</requirement>
639ee511d552 Add biopython requirement
Jim Johnson <jj@umn.edu>
parents: 0
diff changeset
6 <requirement type="python-module">Bio</requirement>
639ee511d552 Add biopython requirement
Jim Johnson <jj@umn.edu>
parents: 0
diff changeset
7 </requirements>
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
8 <command interpreter="python">
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
9 translate_bed_sequences.py --input "$input"
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
10 #if $fa_db:
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
11 --fa_db='$fa_db'
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
12 #end if
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
13 #if $fa_sep:
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
14 --fa_sep='$fa_sep'
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
15 #end if
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
16 #if $id_prefix:
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
17 --id_prefix='$id_prefix'
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
18 #end if
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
19 #if $reference:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
20 --reference $reference
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
21 #else:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
22 --reference ${input.metadata.dbkey}
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
23 #end if
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
24 #if $refsource:
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
25 --refsource $refsource
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
26 #end if
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
27 #if $seqtype:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
28 --seqtype $seqtype
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
29 #end if
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
30 #if $score_name:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
31 --score_name $score_name
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
32 #end if
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
33 #if $filter.filterseqs == 'yes':
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
34 #if $filter.leading_bp:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
35 --leading_bp $filter.leading_bp
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
36 #end if
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
37 #if $filter.trailing_bp:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
38 --trailing_bp $filter.trailing_bp
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
39 #end if
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
40 #else:
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
41 --unfiltered
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
42 #end if
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
43 #if $trim.trimseqs == 'no':
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
44 --untrimmed
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
45 #if str($trim.max_stop_codons) != '':
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
46 --max_stop_codons $trim.max_stop_codons
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
47 #end if
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
48 #end if
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
49 #if str($min_length) != '':
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
50 --min_length $min_length
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
51 #end if
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
52 --bed $translated_bed
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
53 --output "$output"
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
54 </command>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
55 <inputs>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
56 <param name="input" type="data" format="bed" label="BED file with added sequence column"
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
57 help="Output from 'Extract Genomic DNA' run on tophat junctions.bed "/>
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
58 <param name="fa_db" type="text" value="" optional="true" label="fasta ID source, e.g. generic"
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
59 help="Any Compomics application such as PeptideShaker, requires a source">
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
60 </param>
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
61 <param name="fa_sep" type="text" value="" optional="true" label="fasta ID source, e.g. generic"
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
62 help="Only used when a fasta ID source is given, default to the pipe character">
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
63 </param>
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
64 <param name="id_prefix" type="text" value="" optional="true" label="ID prefix for generated IDs"
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
65 help="Can be used to distinguish samples">
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
66 <validator type="regex" message="Allowed chars:a-z A-Z 0-9 _ - |">^[a-zA-Z0-9_-|]*$</validator>
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
67 </param>
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
68 <param name="refsource" type="text" value="Ensembl" optional="true" label="Genome reference source"
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
69 help=""/>
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
70 <param name="reference" type="text" value="" optional="true" label="Genome reference name"
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
71 help="By default, the database metadata will be used."/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
72 <param name="seqtype" type="text" value="" optional="true" label="The SEQTYPE:STATUS to include in the fasta ID lines"
3
3b526a780849 Change default seqtype from pep:novel to pep:splice
Jim Johnson <jj@umn.edu>
parents: 2
diff changeset
73 help="For example: pep:splice"/>
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
74 <param name="score_name" type="text" value="" optional="true" label="Add the bed score field fasta ID line with this tag name"
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
75 help="For example: with the tag name 'depth' and bed score 12: depth:12"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
76 <conditional name="filter">
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
77 <param name="filterseqs" type="select" label="Filter out translations with stop codons before the splice site">
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
78 <option value="yes" selected="true">Yes</option>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
79 <option value="no">No</option>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
80 </param>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
81 <when value="yes">
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
82 <param name="leading_bp" type="integer" value="" min="0" optional="true" label="Stop codon filtering start position base pairs"
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
83 help="Do not reject translation is stop_codons are within base pairs of the BED start position for positive strand"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
84 <param name="trailing_bp" type="integer" value="" min="0" optional="true" label="Stop codon filtering end position base pairs"
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
85 help="Do not reject translation is stop_codons are within base pairs of the BED end position for negative strand"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
86 </when>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
87 <when value="no"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
88 </conditional>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
89 <conditional name="trim">
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
90 <param name="trimseqs" type="select" label="Trim translations to stop codons">
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
91 <option value="yes" selected="true">Yes</option>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
92 <option value="no">No</option>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
93 </param>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
94 <when value="no">
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
95 <param name="max_stop_codons" type="integer" value="" min="0" optional="true" label="Maximum number of stop codons allowed in a translation to be reported"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
96 </when>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
97 </conditional>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
98 <param name="min_length" type="integer" value="" min="0" optional="true" label="Minimum length of a translation to be reported"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
99 </inputs>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
100 <stdio>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
101 <exit_code range="1:" level="fatal" description="Error" />
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
102 </stdio>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
103 <outputs>
5
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
104 <data name="translated_bed" metadata_source="input" format="bed" label="${tool.name} on ${on_string} bed">
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
105 </data>
c626a939eef7 Uploaded
jjohnson
parents: 3
diff changeset
106 <data name="output" metadata_source="input" format="fasta" label="${tool.name} on ${on_string} fasta">
0
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
107 </data>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
108 </outputs>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
109 <tests>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
110 <test>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
111 <param name="input" value="Extract_Genomic_DNA.bed" ftype="bed" dbkey="hg19"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
112 <param name="reference" value="GRCh37"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
113 <param name="seqtype" value="pep:novel"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
114 <param name="score_name" value="depth"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
115 <output name="output" file="translated_bed_sequences.fa"/>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
116 </test>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
117 </tests>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
118 <help>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
119 **Translate BED Sequences**
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
120
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
121 This tool takes a BED input file that has been processed
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
122 by the Galaxy tool "Extract Genomic DNA" to add a 13th column with the transcript sequence.
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
123
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
124 It generates a peptide fasta file with the 3-frame translations of the spliced sequence
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
125 defined by each entry in the input BED file.
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
126
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
127 </help>
57e586ee821e Uploaded
jjohnson
parents:
diff changeset
128 </tool>