34
|
1 <tool id="get_chrom_sizes" name="Get Chromosome Sizes" version="2.8">
|
26
|
2 <description>Generates a TAB Delimited chrom.sizes File from an Inputted FASTA</description>
|
1
|
3 <command>
|
33
|
4
|
|
5 #if $refGenomeSource.genomeSource == "builtin":
|
35
|
6 cat "#chrom\tlength" > $output1
|
|
7 cat ${refGenomeSource.input1_builtin.fields.path} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf ${input2} substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat >> $output1
|
33
|
8 #else:
|
34
|
9 cat ${refGenomeSource.input1_file} | awk '$0 ~ ">" {if (NR > 1) {print c;} c=0; printf ${input2} substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END {print c; }' | cat > $output1
|
29
|
10 #end if
|
|
11
|
1
|
12 </command>
|
0
|
13 <inputs>
|
4
|
14 <conditional name="refGenomeSource">
|
|
15 <param name="genomeSource" type="select" label="Will you select a reference from your history or use a built-in FASTA?">
|
|
16 <option value="builtin">Use a built-in FASTA</option>
|
|
17 <option value="history">Use one from the history</option>
|
|
18 </param>
|
|
19 <when value="builtin">
|
12
|
20 <param name="input1_builtin" type="select" format="fasta" label="Select a reference" help="if your reference of interest is not listed - contact RiboGalaxy team">
|
4
|
21 <options from_data_table="builtin_fastas">
|
|
22 <filter type="sort_by" column="2" />
|
|
23 <validator type="no_options" message="No built-ins are available" />
|
|
24 </options>
|
|
25 </param>
|
|
26 </when>
|
|
27 <when value="history">
|
6
|
28 <param name="input1_file" type="data" format="fasta" label="FASTA File" />
|
4
|
29 </when>
|
|
30 </conditional>
|
|
31 <param name="input2" type="select" label="Chromosome Column Prefix (eg. chr)">
|
|
32 <option value="chr">'chr' (required for upload to GWIPS-Viz)</option>
|
33
|
33 <option value="" selected='true'>None - do not add a prefix</option>
|
4
|
34 </param>
|
3
|
35
|
0
|
36 </inputs>
|
|
37 <outputs>
|
|
38 <data name="output1" format="tabular"/>
|
|
39 </outputs>
|
|
40 <tests>
|
|
41 <test>
|
1
|
42 <param name="input1" value="test.fasta" ftype="fasta" />
|
3
|
43 <param name="input2" value="chr"/>
|
0
|
44 <output name="output1" file="test.fasta.sizes" ftype="tabular" lines_diff="4" />
|
|
45 </test>
|
|
46 </tests>
|
|
47 <help>
|
|
48 **What it does**
|
|
49
|
25
|
50 Creates a chromosome sizes file from a fasta file. This is needed for creating ribosome profiles.
|
|
51
|
|
52 Some commonly used genomes are provided to save on storage.
|
0
|
53 </help>
|
|
54 <citations/>
|
|
55 </tool>
|