|
0
|
1 <tool id="correctGCBias" name="correctGCBias" version="1.0">
|
|
|
2 <description>
|
|
|
3 </description>
|
|
|
4 <requirements>
|
|
|
5 <requirement type="package" version="1.5.1_98e5d8a61431ea8605c0643d991a1a5d8999b4dc">deepTools</requirement>
|
|
|
6 <requirement type="package" version="0.1">ucsc_tools</requirement>
|
|
|
7 </requirements>
|
|
|
8 <command>
|
|
|
9 correctGCBias
|
|
|
10 --bamfile '$bamInput'
|
|
|
11 --species '$species'
|
|
|
12 --GCbiasFrequenciesFile $GCbiasFrequenciesFile
|
|
|
13
|
|
|
14 #if $source.ref_source=="history":
|
|
|
15 --genome $source.input1
|
|
|
16 #else:
|
|
|
17 --genome "${source.input1_2bit.fields.path}"
|
|
|
18 #end if
|
|
|
19
|
|
|
20 #if $advancedOpt.showAdvancedOpt == "yes":
|
|
|
21 #if str($advancedOpt.region.value) != '':
|
|
|
22 --region '$advancedOpt.region'
|
|
|
23 #end if
|
|
|
24
|
|
|
25 --binSize '$advancedOpt.binSize'
|
|
|
26 #end if
|
|
|
27
|
|
|
28 #set newoutFileName="corrected."+str($outFileFormat)
|
|
|
29
|
|
|
30 --correctedFile $newoutFileName; mv $newoutFileName $outFileName
|
|
|
31
|
|
|
32 </command>
|
|
|
33
|
|
|
34 <inputs>
|
|
|
35
|
|
|
36 <param name="GCbiasFrequenciesFile" type="data" format="tabular" label="Output of computeGCBias" />
|
|
|
37
|
|
|
38 <param name="bamInput" format="bam" type="data" label="Input BAM file" help="The BAM file must be sorted and indexed."/>
|
|
|
39
|
|
|
40 <param name="species" type="text" value="" label="Species name abbreviation" />
|
|
|
41
|
|
|
42 <conditional name="source">
|
|
|
43 <param name="ref_source" type="select" label="Reference genome">
|
|
|
44 <option value="cached">locally cached</option>
|
|
|
45 <option value="history">in your history</option>
|
|
|
46 </param>
|
|
|
47 <when value="cached">
|
|
|
48 <param name="input1_2bit" type="select" label="Using reference genome" help="If your genome of interest is not listed, contact your Galaxy team">
|
|
|
49 <options from_data_table="deepTools_seqs" />
|
|
|
50 </param>
|
|
|
51 </when>
|
|
|
52 <when value="history">
|
|
|
53 <param name="input1" type="data" format="twobit" label="Select a reference dataset in 2bit format" />
|
|
|
54 </when>
|
|
|
55 </conditional>
|
|
|
56
|
|
|
57 <param name="outFileFormat" type="select" label="File format of the output">
|
|
|
58 <option value="bam">bam</option>
|
|
|
59 <option value="bw">bigwig</option>
|
|
|
60 <option value="bg">bedgraph</option>
|
|
|
61 </param>
|
|
|
62
|
|
|
63 <conditional name="advancedOpt">
|
|
|
64 <param name="showAdvancedOpt" type="select" label="Show advanced options" >
|
|
|
65 <option value="no" selected="true">no</option>
|
|
|
66 <option value="yes">yes</option>
|
|
|
67 </param>
|
|
|
68 <when value="no" />
|
|
|
69 <when value="yes">
|
|
|
70 <param name="region" type="text" value=""
|
|
|
71 label="Region of the genome to limit the operation to"
|
|
|
72 help="This is useful when testing parameters to reduce the computing time. The format is chr:start:end, for example "chr10" or "chr10:456700:891000"" />
|
|
|
73
|
|
|
74 <param name="binSize" type="integer" value="50" min="1"
|
|
|
75 label="Bin size in bp"
|
|
|
76 help="Size of the bins in bp for the ouput of the bigwig/bedgraph file."/>
|
|
|
77 </when>
|
|
|
78 </conditional>
|
|
|
79 </inputs>
|
|
|
80
|
|
|
81 <outputs>
|
|
|
82 <data format="bam" name="outFileName">
|
|
|
83 <change_format>
|
|
|
84 <when input="outFileFormat" value="bw" format="bigwig" />
|
|
|
85 <when input="outFileFormat" value="bam" format="bam" />
|
|
|
86 <when input="outFileFormat" value="bg" format="bedgraph" />
|
|
|
87 </change_format>
|
|
|
88 </data>
|
|
|
89 </outputs>
|
|
|
90 <help>
|
|
|
91
|
|
|
92 **What it does**
|
|
|
93
|
|
|
94 Computes the GC bias ussing Benjamini's method [citation]. The resulting GC
|
|
|
95 bias can later be used to plot the bias or to correct the bias.
|
|
|
96
|
|
|
97 -----
|
|
|
98
|
|
|
99 .. class:: infomark
|
|
|
100
|
|
|
101 Please acknowledge that this tool **is still in development** and we will be very happy to receive feedback from the users. If you run into any trouble please sent an email to `Fidel Ramirez`_.
|
|
|
102
|
|
|
103 This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
|
|
|
104
|
|
|
105
|
|
|
106 .. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
|
|
|
107 .. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
|
|
|
108 .. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
|
|
|
109
|
|
|
110 </help>
|
|
|
111
|
|
|
112 </tool>
|
