|
4
|
1 <tool id="blockclust" name="BlockClust" version="0.1">
|
|
|
2 <description>Non-coding RNA clustering from deep sequencing read profiles</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="1.0">blockclust</requirement>
|
|
5
|
5 <requirement type="package" version="3.0.1">R_3_0_1</requirement>
|
|
|
6 <requirement type="package" version="0.1.19 ">samtools</requirement>
|
|
|
7 <requirement type="R-module">ggplot2</requirement>
|
|
|
8 <requirement type="R-module">ape</requirement>
|
|
9
|
9 <requirement type="set_environment">BLOCKCLUST_DATA_PATH</requirement>
|
|
4
|
10 </requirements>
|
|
|
11 <!--<version_command> -version</version_command>-->
|
|
|
12 <command>
|
|
|
13 #if str($mode.operation) == "pre":
|
|
|
14 #set $outputdir = $tags_bed.extra_files_path
|
|
|
15 BlockClustPipeLine.pl -m PRE -bam $mode.reads_bam -tbed $tags_bed;
|
|
|
16 #end if
|
|
|
17
|
|
|
18 #if str($mode.operation) == "clust":
|
|
|
19 #set $outputdir = $clusters.extra_files_path
|
|
|
20 #set $accept_bed=list()
|
|
|
21 #set $reject_bed=list()
|
|
|
22 ## prepare annotations
|
|
|
23 #if str($mode.reference) == "hg19":
|
|
9
|
24 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/hg19/hg19.accept.bed")
|
|
|
25 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/hg19/hg19.reject.bed")
|
|
4
|
26 #elif str($mode.reference) == "mm10":
|
|
9
|
27 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/mm10/mm10.accept.bed")
|
|
|
28 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/mm10/mm10.reject.bed")
|
|
4
|
29 #elif str($mode.reference) == "dm3":
|
|
9
|
30 $accept_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/dm3/dm3.accept.bed")
|
|
|
31 $reject_bed.append("\$BLOCKCLUST_DATA_PATH/annotations/dm3/dm3.reject.bed")
|
|
4
|
32 #end if
|
|
9
|
33 BlockClustPipeLine.pl -m TEST -f SEQUENCE -c \$BLOCKCLUST_DATA_PATH/blockclust.config
|
|
4
|
34 -t $mode.input_bbo
|
|
|
35 -a #echo ''.join( $accept_bed )
|
|
|
36 -r #echo ''.join( $reject_bed )
|
|
|
37 -o $outputdir;
|
|
|
38 cp #echo os.path.join($outputdir, 'mcl_clusters','all_clusters.bed')# $clusters;
|
|
|
39 cp #echo os.path.join($outputdir, 'hclust_tree.pdf')# $hclust_plot;
|
|
|
40 cp #echo os.path.join($outputdir, 'discretized.gspan.tab')# $sim_tab_out;
|
|
|
41 #end if
|
|
|
42
|
|
|
43 #if str($mode.operation) == "post":
|
|
|
44 #set $outputdir = $clusters_bed.extra_files_path
|
|
9
|
45 BlockClustPipeLine.pl -m POST -cbed $mode.clusters_bed -cm $mode.cmsearch_out -tab $mode.sim_tab_in -o $outputdir -rfam \$BLOCKCLUST_DATA_PATH/rfam_map.txt;
|
|
4
|
46 cp #echo os.path.join($outputdir, 'cluster_distribution.pdf')# $cluster_dist;
|
|
|
47 #end if
|
|
|
48 </command>
|
|
|
49 <inputs>
|
|
|
50 <conditional name="mode">
|
|
|
51 <param name="operation" type="select" label="Select mode of operation">
|
|
|
52 <option value="pre">Pre-processing </option>
|
|
|
53 <option value="clust">Clustering</option>
|
|
|
54 <option value="post">Post-processing</option>
|
|
|
55 </param>
|
|
|
56 <when value="pre">
|
|
|
57 <param name="reads_bam" type="data" format="bam" label="BAM file containing alignments" />
|
|
|
58 </when>
|
|
|
59 <when value="clust">
|
|
|
60 <param name="input_bbo" type="data" format="tabular" label="Input blockgroups file" />
|
|
|
61 <param name="reference" type="select" label="Select reference genome">
|
|
|
62 <option value="hg19">Human (hg19)</option>
|
|
|
63 <option value="mm10">Mouse (mm10)</option>
|
|
|
64 <option value="dm3">Fly (dm3)</option>
|
|
|
65 </param>
|
|
|
66 </when>
|
|
|
67 <when value="post">
|
|
5
|
68 <param name="cmsearch_out" type="data" format="tabular" label="Output of cmsearch tool" />
|
|
4
|
69 <param name="clusters_bed" type="data" format="bed" label="BED file containing clusters (output of BlockClust)" />
|
|
|
70 <param name="sim_tab_in" type="data" format="tabular" label="Pairwise similarities file" />
|
|
|
71 </when>
|
|
|
72 </conditional>
|
|
|
73 </inputs>
|
|
|
74
|
|
|
75 <outputs>
|
|
|
76 <data format="bed" name="tags_bed" label="BlockClust: BAM to BED on ${on_string}">
|
|
|
77 <filter> mode["operation"]=="pre"</filter>
|
|
|
78 </data>
|
|
5
|
79 <data format="pdf" name="hclust_plot" label="BlockClust: Hierarchical clustering plot on ${on_string}" >
|
|
4
|
80 <filter> mode["operation"]=="clust"</filter>
|
|
|
81 </data>
|
|
5
|
82 <data format="bed" name="clusters" label="BlockClust: Clustering BED on ${on_string}">
|
|
4
|
83 <filter> mode["operation"]=="clust"</filter>
|
|
|
84 </data>
|
|
|
85 <data format="tabular" name="sim_tab_out" label="BlockClust: Pairwise similarities on ${on_string}">
|
|
|
86 <filter> mode["operation"]=="clust"</filter>
|
|
|
87 </data>
|
|
|
88 <data format="tabular" name="analysis" label="BlockClust: Cluster analysis on ${on_string}" >
|
|
|
89 <filter> mode["operation"]=="post"</filter>
|
|
|
90 </data>
|
|
|
91 <data format="pdf" name="cluster_dist" label="BlockClust: Cluster distribution on ${on_string}" >
|
|
|
92 <filter> mode["operation"]=="post"</filter>
|
|
|
93 </data>
|
|
|
94 </outputs>
|
|
|
95 <help>
|
|
|
96
|
|
|
97 .. class:: infomark
|
|
|
98
|
|
|
99 **What it does**
|
|
|
100
|
|
|
101 Clusters the read profiles (i.e., blockgroups) from the blockbuster tool.
|
|
|
102
|
|
|
103 **Inputs**
|
|
|
104
|
|
|
105 BlockClust needs output of tool blockbuster as input
|
|
|
106
|
|
|
107 **Output**
|
|
|
108
|
|
|
109 BlockClust produces a fasta file containing clusters.
|
|
|
110
|
|
|
111 ------
|
|
|
112
|
|
|
113 **Licenses**
|
|
|
114
|
|
|
115 If **BlockClust** is used to obtain results for scientific publications it
|
|
|
116 should be cited as [1]_.
|
|
|
117
|
|
|
118 **References**
|
|
|
119
|
|
|
120 .. [1]
|
|
|
121
|
|
|
122
|
|
|
123 </help>
|
|
|
124 </tool>
|
