|
0
|
1 <?xml version="1.0"?>
|
|
|
2 <tool id="gbtofasta" name="gbToFasta" version="1.0">
|
|
|
3 <description>Convert GenBank records to fasta and Create table with coding regions information for each mRNA record</description>
|
|
|
4 <requirements>
|
|
|
5 <requirement type="package" version="1.0">ucsc_tools_340_for_BLAT</requirement>
|
|
|
6 </requirements>
|
|
|
7 <command detect_errors="exit_code"><![CDATA[
|
|
|
8 gbToFaRa
|
|
|
9 /dev/null
|
|
|
10 '${output_fa}'
|
|
|
11 output_ra
|
|
|
12 output_ta
|
|
|
13 '${gbfile}'
|
|
|
14 #if $cds == "yes"
|
|
|
15 && raToTab
|
|
|
16 -cols=acc,cds
|
|
|
17 output_ra
|
|
|
18 '${output_cds}'
|
|
|
19 #end if
|
|
|
20
|
|
|
21 ]]></command>
|
|
|
22 <inputs>
|
|
|
23 <param type="data" name="gbfile" format="genbank" />
|
|
|
24 <param type="select" name="cds" label="Create table with coding regions information for each mRNA record" >
|
|
|
25 <option value="yes">Yes</option>
|
|
|
26 <option value="no">No</option>
|
|
|
27 </param>
|
|
|
28
|
|
|
29 </inputs>
|
|
|
30 <outputs>
|
|
|
31 <data format="fasta" name="output_fa" label="${gbfile.name}:fasta"></data>
|
|
|
32 <data format="fasta" name="output_cds" label="${gbfile.name}:cds">
|
|
|
33 <filter>cds == "yes"</filter>
|
|
|
34 </data>
|
|
|
35 </outputs>
|
|
|
36 <tests>
|
|
|
37 <test>
|
|
|
38 <param name="gbfile" value="Gallus_gallus_RefSeq.gb" />
|
|
|
39 <param name="cds" value="yes" />
|
|
|
40 <output name="output_fa" value="Gallus_gallus_RefSeq.fa" />
|
|
|
41 <output name="output_cds" value="Gallus_gallus_RefSeq.cds" />
|
|
|
42 </test>
|
|
|
43 <test>
|
|
|
44 <param name="gbfile" value="Gallus_gallus_RefSeq.gb" />
|
|
|
45 <param name="cds" value="no" />
|
|
|
46 <output name="output_fa" value="Gallus_gallus_RefSeq.fa" />
|
|
|
47 </test>
|
|
|
48 </tests>
|
|
|
49 <help>
|
|
|
50 <![CDATA[
|
|
|
51 gbToFasta
|
|
|
52 =========
|
|
|
53
|
|
|
54 Convert GenBank records to fasta
|
|
|
55 ---------------------------------
|
|
|
56
|
|
|
57 gbToFaRa - Convert GenBank flat format file to an fa file containing
|
|
|
58 the sequence data, an ra file containing other relevant info and
|
|
|
59 a ta file containing summary statistics.
|
|
|
60 usage:
|
|
|
61 gbToFaRa filterFile faFile raFile taFile genBankFile(s)
|
|
|
62 where filterFile is definition of which records and fields
|
|
|
63 use /dev/null if you want no filtering.
|
|
|
64
|
|
|
65 gbToFaRa /dev/null ${gbfile.fa} \
|
|
|
66 ${gbfile.ra} ${gbfile.ta} ${gbfile}
|
|
|
67
|
|
|
68 Create table with coding regions information for each mRNA record
|
|
|
69 -----------------------------------------------------------------
|
|
|
70 raToTab - Convert ra file to table.
|
|
|
71 raToTab -cols=acc,cds ${gbfile.ra} ${gbfile.cds}
|
|
|
72
|
|
|
73 Source code:
|
|
|
74 ============
|
|
|
75
|
|
|
76 http://hgdownload.cse.ucsc.edu/admin/exe/
|
|
|
77
|
|
|
78 ]]></help>
|
|
|
79 <citations>
|
|
|
80 <citation type="bibtex">@article{kent2002blat,
|
|
|
81 title={BLAT—the BLAST-like alignment tool},
|
|
|
82 author={Kent, W James},
|
|
|
83 journal={Genome research},
|
|
|
84 volume={12},
|
|
|
85 number={4},
|
|
|
86 pages={656--664},
|
|
|
87 year={2002},
|
|
|
88 publisher={Cold Spring Harbor Lab}
|
|
|
89 }</citation>
|
|
|
90 </citations>
|
|
|
91 </tool>
|
|
|
92
|
|
|
93
|
|
|
94
|
|
|
95
|
|
|
96 |
