annotate blockclust.xml @ 1:a92d8689aa30 draft

Uploaded
author rnateam
date Tue, 08 Jul 2014 11:59:49 -0400
parents 768a67f6ae18
children f973ec6e5192
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
1 <tool id="blockclust" name="BlockClust" version="1.0.0">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
2 <description>Non-coding RNA clustering from deep sequencing read profiles</description>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
3 <requirements>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
4 <requirement type="package" version="1.0">blockclust</requirement>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
5 <requirement type="package" version="1.1">eden</requirement>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
6 <requirement type="package" version="3.0.1">R_3_0_1</requirement>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
7 <requirement type="package" version="0.1.19">samtools</requirement>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
8 <requirement type="package" version="12.135">mcl</requirement>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
9 <requirement type="package" version="1.0">blockclust_rlibs</requirement>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
10 </requirements>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
11 <version_command>echo '1.0'</version_command>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
12 <command>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
13
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
14 #if str($tool_mode.operation) == "pre":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
15 BlockClustPipeLine.pl -m PRE -bam $tool_mode.reads_bam -tbed $tags_bed
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
16 #elif str($tool_mode.operation) == "clust":
1
a92d8689aa30 Uploaded
rnateam
parents: 0
diff changeset
17 #set $outputdir = $clusters.files_path
0
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
18 #set $accept_bed=list()
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
19 #set $reject_bed=list()
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
20 ## prepare annotations
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
21 #if str($tool_mode.reference) == "hg19":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
22 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/hg19/hg19.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
23 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/hg19/hg19.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
24 #elif str($tool_mode.reference) == "mm10":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
25 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/mm10/mm10.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
26 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/mm10/mm10.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
27 #elif str($tool_mode.reference) == "dm3":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
28 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/dm3/dm3.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
29 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/dm3/dm3.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
30 #elif str($tool_mode.reference) == "rheMac3":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
31 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/rheMac3/rheMac3.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
32 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/rheMac3/rheMac3.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
33 #elif str($tool_mode.reference) == "panTro4":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
34 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/panTro3/panTro4.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
35 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/panTro3/panTro4.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
36 #elif str($tool_mode.reference) == "xenTro3":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
37 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/xenTro3/xenTro3.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
38 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/xenTro3/xenTro3.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
39 #elif str($tool_mode.reference) == "celWS235":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
40 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/celWS235/celWS235.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
41 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/celWS235/celWS235.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
42 #elif str($tool_mode.reference) == "tair10":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
43 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/tair10/tair10.accept.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
44 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/tair10/tair10.reject.bed")
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
45 #end if
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
46 BlockClustPipeLine.pl -m TEST -f SEQUENCE -c \$BLOCKCLUST_DATA_PATH/blockclust.config
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
47 -t $tool_mode.input_bbo
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
48 -a #echo ''.join( $accept_bed )
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
49 -r #echo ''.join( $reject_bed )
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
50 -o $outputdir
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
51 #if str($tool_mode.pred.enable_pred) == "yes":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
52 -p
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
53 -pm $tool_mode.pred.pred_mode
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
54 -md \$BLOCKCLUST_DATA_PATH/models
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
55 -x 3;
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
56 #if str($tool_mode.pred.pred_mode) == "nearest_neighbour":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
57 cp #echo os.path.join($outputdir,'nearest_neighbour_predictions.txt')# $nearest_neighbour_pred_bed;
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
58 #elif str($tool_mode.pred.pred_mode) == "model_based":
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
59 cp #echo os.path.join($outputdir,'model_based_predictions.txt')# $model_based_pred_bed;
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
60 #end if
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
61 #else:
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
62 ;
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
63 #end if
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
64
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
65 cp #echo os.path.join($outputdir, 'mcl_clusters','all_clusters.bed')# $clusters;
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
66 cp #echo os.path.join($outputdir, 'hclust_tree.pdf')# $hclust_plot;
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
67 cp #echo os.path.join($outputdir, 'discretized.gspan.tab')# $sim_tab_out
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
68 #elif str($tool_mode.operation) == "post":
1
a92d8689aa30 Uploaded
rnateam
parents: 0
diff changeset
69 BlockClustPipeLine.pl -m POST -cbed $tool_mode.clusters_bed -cm $tool_mode.cmsearch_out -tab $tool_mode.sim_tab_in -o ./;
0
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
70 #end if
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
71 </command>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
72 <inputs>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
73 <conditional name="tool_mode">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
74 <param name="operation" type="select" label="Select mode of operation">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
75 <option value="pre">Pre-processing </option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
76 <option value="clust">Clustering</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
77 <option value="post">Post-processing</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
78 </param>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
79 <when value="pre">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
80 <param name="reads_bam" type="data" format="bam" label="BAM file containing alignments" />
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
81 </when>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
82 <when value="clust">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
83 <param name="input_bbo" type="data" format="tabular" label="Input blockgroups file" />
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
84 <param name="reference" type="select" label="Select reference genome">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
85 <option value="hg19">Human (hg19)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
86 <option value="mm10">Mouse (mm10)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
87 <option value="dm3">Fly (dm3)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
88 <option value="rheMac3">Monkey (rheMac3)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
89 <option value="panTro4">Chimp (panTro4)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
90 <option value="xenTro3">Frog (xenTro3)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
91 <option value="celWS235">C. elegans (celWS235)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
92 <option value="tair10">Arabidopsis thaliana (tair10)</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
93 </param>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
94 <conditional name="pred">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
95 <param name="enable_pred" type="select" label="Would you like to perform classification?">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
96 <option value="no">No</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
97 <option value="yes">Yes</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
98 </param>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
99 <when value="yes">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
100 <param name="pred_mode" type="select" label="Mode of classification">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
101 <option value="model_based">Model based</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
102 <option value="nearest_neighbour">Nearest neighbour</option>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
103 </param>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
104 </when>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
105 </conditional>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
106 </when>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
107 <when value="post">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
108 <param name="cmsearch_out" type="data" format="tabular" label="Output of cmsearch tool" />
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
109 <param name="clusters_bed" type="data" format="bed" label="BED file containing clusters (output of BlockClust)" />
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
110 <param name="sim_tab_in" type="data" format="tabular" label="Pairwise similarities file" />
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
111 </when>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
112 </conditional>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
113 </inputs>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
114
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
115 <outputs>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
116 <data format="bed" name="tags_bed" label="BlockClust: BAM to BED on ${on_string}">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
117 <filter> tool_mode["operation"]=="pre"</filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
118 </data>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
119 <data format="pdf" name="hclust_plot" label="BlockClust: Hierarchical clustering plot on ${on_string}" >
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
120 <filter> tool_mode["operation"]=="clust"</filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
121 </data>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
122 <data format="bed" name="clusters" label="BlockClust: Clustering BED on ${on_string}">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
123 <filter> tool_mode["operation"]=="clust"</filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
124 </data>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
125 <data format="bed" name="model_based_pred_bed" label="BlockClust: Model based prediction BED on ${on_string}">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
126 <filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
127 ((
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
128 tool_mode["operation"] == 'clust' and
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
129 tool_mode["pred"]["enable_pred"] == "yes" and
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
130 tool_mode["pred"]["pred_mode"] == "model_based"
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
131 ))
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
132 </filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
133 </data>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
134 <data format="bed" name="nearest_neighbour_pred_bed" label="BlockClust: Nearest neighbour prediction BED on ${on_string}">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
135 <filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
136 ((
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
137 tool_mode["operation"] == 'clust' and
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
138 tool_mode["pred"]["enable_pred"] == "yes" and
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
139 tool_mode["pred"]["pred_mode"] == "nearest_neighbour"
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
140 ))
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
141 </filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
142 </data>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
143 <data format="tabular" name="sim_tab_out" label="BlockClust: Pairwise similarities on ${on_string}">
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
144 <filter> tool_mode["operation"]=="clust"</filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
145 </data>
1
a92d8689aa30 Uploaded
rnateam
parents: 0
diff changeset
146 <data format="pdf" name="cluster_dist" from_work_dir="cluster_distribution.pdf" label="BlockClust: Cluster distribution on ${on_string}" >
0
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
147 <filter> tool_mode["operation"]=="post"</filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
148 </data>
1
a92d8689aa30 Uploaded
rnateam
parents: 0
diff changeset
149 <data format="pdf" name="cluster_hclust" from_work_dir="hclust_tree_clusters.pdf" label="BlockClust: Hierarchical clustering plot of cluster centroids on ${on_string}" >
0
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
150 <filter> tool_mode["operation"]=="post"</filter>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
151 </data>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
152 </outputs>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
153 <help>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
154
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
155 .. class:: infomark
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
156
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
157 **What it does**
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
158
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
159 Clusters the read profiles (i.e., blockgroups) from the blockbuster tool.
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
160
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
161 **Inputs**
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
162
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
163 BlockClust needs output of tool blockbuster as input
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
164
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
165 **Output**
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
166
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
167 BlockClust produces a fasta file containing clusters.
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
168
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
169 ------
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
170
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
171 **Licenses**
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
172
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
173 If **BlockClust** is used to obtain results for scientific publications it
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
174 should be cited as [1]_.
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
175
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
176 **References**
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
177
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
178 .. [1]
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
179
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
180
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
181 </help>
768a67f6ae18 Uploaded
rnateam
parents:
diff changeset
182 </tool>