|
0
|
1 <tool id="diffbind" name="DiffBind" version="2.0.1">
|
|
|
2 <description> differential binding analysis of ChIP-Seq peak data</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="binary">Rscript</requirement>
|
|
|
5 <requirement type="R-module">DESeq2</requirement>
|
|
|
6 <requirement type="package" version="3.0.1">R_3_0_1</requirement>
|
|
|
7 <requirement type="package" version="1.0.17">deseq2</requirement>
|
|
6
|
8 <requirement type="package" version="1.8.3">diffbind</requirement>
|
|
11
|
9 <requirement type="set_environment">DIFFBIND_SCRIPT_PATH</requirement>
|
|
0
|
10 </requirements>
|
|
10
|
11 <command>
|
|
|
12 ## seems that diffbind also needs file extensions to work properly
|
|
|
13 #set $counter = 1
|
|
|
14 #for $sample in $samples:
|
|
|
15 ln -s $sample.bamreads #echo str($counter) + "_bamreads.bam"# &&
|
|
|
16 ln -s ${sample.bamreads.metadata.bam_index} #echo str($counter) + "_bamreads.bai"# &&
|
|
|
17 ln -s $sample.bamcontrol #echo str($counter) + "_bamcontrol.bam"# &&
|
|
|
18 ln -s ${sample.bamcontrol.metadata.bam_index} #echo str($counter) + "_bamcontrol.bai"# &&
|
|
|
19 #set $counter = $counter + 1
|
|
|
20 #end for
|
|
|
21
|
|
12
|
22 Rscript \$DIFFBIND_SCRIPT_PATH/diffbind.R
|
|
0
|
23 -i $infile
|
|
|
24 -o $outfile
|
|
|
25 -p $plots
|
|
|
26 -f $format
|
|
|
27 </command>
|
|
|
28 <stdio>
|
|
|
29 <regex match="Execution halted"
|
|
|
30 source="both"
|
|
|
31 level="fatal"
|
|
|
32 description="Execution halted." />
|
|
|
33 <regex match="Input-Error 01"
|
|
|
34 source="both"
|
|
|
35 level="fatal"
|
|
|
36 description="Error in your input parameters: Make sure you only apply factors to selected samples." />
|
|
|
37 <regex match="Error in"
|
|
|
38 source="both"
|
|
|
39 level="fatal"
|
|
|
40 description="An undefined error occured, please check your intput carefully and contact your administrator." />
|
|
|
41 </stdio>
|
|
|
42 <inputs>
|
|
|
43 <repeat name="samples" title="Samples" min="2">
|
|
|
44 <param name="sample_id" type="text" value="Sample ID" label="Specify a sample id" help="e.g. BT474.1-" />
|
|
|
45 <param name="tissue" type="text" value="Tissue" label="Specify the tissue" help="e.g. BT474" />
|
|
|
46 <param name="factor" type="text" value="Factor Name" label="Specify a factor name" help="e.g. ER" />
|
|
|
47 <param name="condition" type="text" value="Condition" label="Specify the condition" help="e.g. Resistent" />
|
|
3
|
48 <param name="replicate" type="integer" value="1" label="Specify the replicate number" help="e.g. 1" />
|
|
0
|
49 <param format="bam" name="bamreads" type="data" label="Read BAM file" help="Specify the Read BAM file, used for Peak calling."/>
|
|
|
50 <param format="bam" name="bamcontrol" type="data" label="Control BAM file" help="Specify the Control BAM file, used for Peak calling."/>
|
|
4
|
51 <param format="bed" name="peaks" type="data" label="Peak file" help="Result of your Peak calling experiment."/>
|
|
0
|
52 </repeat>
|
|
|
53 <param name="pdf" type="boolean" truevalue="" falsevalue="" checked="true"
|
|
|
54 label="Visualising the analysis results"
|
|
|
55 help="output an additional PDF files" />
|
|
|
56 <param name="format" type="select" label="Output Format">
|
|
|
57 <option value="bed">BED</option>
|
|
|
58 <option value="gff">GFF</option>
|
|
|
59 <option value="wig">WIG</option>
|
|
|
60 </param>
|
|
|
61 </inputs>
|
|
|
62 <configfiles>
|
|
9
|
63 <configfile name="infile">
|
|
0
|
64 SampleID,Tissue,Factor,Condition,Replicate,bamReads,bamControl,Peaks
|
|
10
|
65 #set $counter = 1
|
|
0
|
66 #for $sample in $samples:
|
|
10
|
67 $sample.sample_id,$sample.tissue,$sample.factor,$sample.condition,$sample.replicate,#echo str($counter) + '_bamreads.bam'#,#echo str($counter) + '_bamcontrol.bam'#,$sample.peaks
|
|
|
68 #set $counter = $counter + 1
|
|
0
|
69 #end for
|
|
9
|
70 </configfile>
|
|
0
|
71 </configfiles>
|
|
|
72 <outputs>
|
|
|
73 <data format="bed" name="outfile" label="Differential binding sites on ${on_string}">
|
|
|
74 <change_format>
|
|
|
75 <when input="format" value="wig" format="wig" />
|
|
|
76 <when input="format" value="gff" format="gff" />
|
|
|
77 </change_format>
|
|
|
78 </data>
|
|
9
|
79 <data format="pdf" name="plots" label="Differential binding sites on ${on_string}">
|
|
0
|
80 <filter>pdf == True</filter>
|
|
|
81 </data>
|
|
|
82 </outputs>
|
|
|
83
|
|
|
84 <help>
|
|
|
85
|
|
|
86 .. class:: infomark
|
|
|
87
|
|
|
88 **What it does**
|
|
|
89
|
|
|
90
|
|
|
91 ------
|
|
|
92
|
|
|
93 **References**
|
|
|
94
|
|
|
95 DiffBind_ Authors: Rory Stark, Gordon Brown
|
|
|
96
|
|
|
97 .. _DiffBind: http://www.bioconductor.org/packages/release/bioc/html/DiffBind.html
|
|
|
98
|
|
|
99 Wrapper authors: Bjoern Gruening, Pavankumar Videm
|
|
|
100
|
|
|
101 </help>
|
|
|
102 </tool>
|
