|
0
|
1 <?xml version="1.0"?>
|
|
|
2 <tool id="translate_bed_sequences" name="Translate BED Sequences" version="0.0.1">
|
|
|
3 <description>3 frame translation of BED augmented with a sequence column</description>
|
|
|
4 <command interpreter="python">translate_bed_sequences.py --input "$input"
|
|
|
5 #if $reference:
|
|
|
6 --reference $reference
|
|
|
7 #else:
|
|
|
8 --reference ${input.metadata.dbkey}
|
|
|
9 #end if
|
|
|
10 #if $seqtype:
|
|
|
11 --seqtype $seqtype
|
|
|
12 #end if
|
|
|
13 #if $score_name:
|
|
|
14 --score_name $score_name
|
|
|
15 #end if
|
|
|
16 #if $filter.filterseqs == 'yes':
|
|
|
17 #if $filter.leading_bp:
|
|
|
18 --leading_bp $filter.leading_bp
|
|
|
19 #end if
|
|
|
20 #if $filter.trailing_bp:
|
|
|
21 --trailing_bp $filter.trailing_bp
|
|
|
22 #end if
|
|
|
23 #else:
|
|
|
24 --unfiltered
|
|
|
25 #end if
|
|
|
26 #if $trim.trimseqs == 'no':
|
|
|
27 --untrimmed
|
|
|
28 #if $trim.max_stop_codons != None:
|
|
|
29 --max_stop_codons $trim.max_stop_codons
|
|
|
30 #end if
|
|
|
31 #end if
|
|
|
32 #if $min_length:
|
|
|
33 --min_length $min_length
|
|
|
34 #end if
|
|
|
35 --output "$output"
|
|
|
36 </command>
|
|
|
37 <inputs>
|
|
|
38 <param name="input" type="data" format="bed" label="BED file with added sequence column"
|
|
|
39 help="Output from 'Extract Genomic DNA' run on tophat junctions.bed "/>
|
|
|
40 <param name="reference" type="text" value="" optional="true" label="Genome reference name"
|
|
|
41 help="By default, the database metadata will be used."/>
|
|
|
42 <param name="seqtype" type="text" value="" optional="true" label="The SEQTYPE:STATUS to include in the fasta ID lines"
|
|
|
43 help="For example: pep:novel"/>
|
|
|
44 <param name="score_name" type="text" value="" optional="true" label="Add the bed score field fasta ID line with this tag name"
|
|
|
45 help="For example: with the tag name 'depth' and bed score 12: depth:12"/>
|
|
|
46 <conditional name="filter">
|
|
|
47 <param name="filterseqs" type="select" label="Filter out translations with stop codons before the splice site">
|
|
|
48 <option value="yes" selected="true">Yes</option>
|
|
|
49 <option value="no">No</option>
|
|
|
50 </param>
|
|
|
51 <when value="yes">
|
|
|
52 <param name="leading_bp" type="integer" value="" min="0" optional="true" label="Stop codon filtering start position base pairs"
|
|
|
53 help="Do not reject translation is stop_codons are within base pairs of the BED start position for positive strand"/>
|
|
|
54 <param name="trailing_bp" type="integer" value="" min="0" optional="true" label="Stop codon filtering end position base pairs"
|
|
|
55 help="Do not reject translation is stop_codons are within base pairs of the BED end position for negative strand"/>
|
|
|
56 </when>
|
|
|
57 <when value="no"/>
|
|
|
58 </conditional>
|
|
|
59 <conditional name="trim">
|
|
|
60 <param name="trimseqs" type="select" label="Trim translations to stop codons">
|
|
|
61 <option value="yes" selected="true">Yes</option>
|
|
|
62 <option value="no">No</option>
|
|
|
63 </param>
|
|
|
64 <when value="no">
|
|
|
65 <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"/>
|
|
|
66 </when>
|
|
|
67 </conditional>
|
|
|
68 <param name="min_length" type="integer" value="" min="0" optional="true" label="Minimum length of a translation to be reported"/>
|
|
|
69 </inputs>
|
|
|
70 <stdio>
|
|
|
71 <exit_code range="1:" level="fatal" description="Error" />
|
|
|
72 </stdio>
|
|
|
73 <outputs>
|
|
|
74 <data name="output" metadata_source="input" format="fasta" label="${tool.name} on ${on_string}">
|
|
|
75 <filter>'found' in str(outputs)</filter>
|
|
|
76 </data>
|
|
|
77 </outputs>
|
|
|
78 <tests>
|
|
|
79 <test>
|
|
|
80 <param name="input" value="Extract_Genomic_DNA.bed" ftype="bed" dbkey="hg19"/>
|
|
|
81 <param name="reference" value="GRCh37"/>
|
|
|
82 <param name="seqtype" value="pep:novel"/>
|
|
|
83 <param name="score_name" value="depth"/>
|
|
|
84 <output name="output" file="translated_bed_sequences.fa"/>
|
|
|
85 </test>
|
|
|
86 </tests>
|
|
|
87 <help>
|
|
|
88 **Translate BED Sequences**
|
|
|
89
|
|
|
90 This tool takes a BED input file that has been processed
|
|
|
91 by the Galaxy tool "Extract Genomic DNA" to add a 13th column with the transcript sequence.
|
|
|
92
|
|
|
93 It generates a peptide fasta file with the 3-frame translations of the spliced sequence
|
|
|
94 defined by each entry in the input BED file.
|
|
|
95
|
|
|
96
|
|
|
97 </help>
|
|
|
98 </tool>
|
