# HG changeset patch
# User bgruening
# Date 1375551732 14400
# Node ID 09b65c12a75aa7ed424e43c88c5671952d7a9086
Uploaded
diff -r 000000000000 -r 09b65c12a75a bamCompare.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/bamCompare.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,210 @@
+
+ Normalize and compare two BAM files to output ratio, log2ratio or difference.
+
+ deepTools
+ numpy
+ ucsc_tools
+
+
+ bamCompare
+ --bamfile1 '$bamFile1'
+ -bai1 '${bamFile1.metadata.bam_index}'
+ --bamfile2 '$bamFile2'
+ -bai2 '${bamFile2.metadata.bam_index}'
+
+ --outFileName '$outFileName'
+ --outFileFormat '$outFileFormat'
+
+ --fragmentLength $fragmentLength
+ --binSize $binSize
+
+ #if $scaling.method == 'SES':
+ --scaleFactorsMethod SES
+ --sampleLength $scaling.sampleLength
+ #elif $scaling.method == 'readCount':
+ --scaleFactorsMethod readCount
+ #elif $scaling.method == 'own':
+ --scaleFactors '$scaling.scaleFactor1:$scaling.scaleFactor2'
+ #end if
+
+ --ratio $comparison.type
+
+
+ #if $comparison.type=='subtract':
+ #if $comparison.normalization.type=='rpkm':
+ --normalizeUsingRPKM
+ #elif $comparison.normalization.type=='1x':
+ --normalizeTo1x $comparison.normalization.normalizeTo1x
+ #end if
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if $advancedOpt.smoothLength:
+ --smoothLength '$advancedOpt.smoothLength'
+ #end if
+
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+ $advancedOpt.doNotExtendPairedEnds
+ $advancedOpt.ignoreDuplicates
+
+ #if $advancedOpt.minMappingQuality:
+ --minMappingQuality '$advancedOpt.minMappingQuality'
+ #end if
+
+ --missingDataAsZero $advancedOpt.missingDataAsZero
+
+ #end if
+ --numberOfProcessors 4
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+This tool compares two BAM files based on the number of mapped reads. To
+compare the BAM files the genome is partitioned into bins of equal size, then
+the number of reads found in each BAM file are counted for such bins and
+finally a summarizing value is reported. This vaule can be the ratio of the
+number of reads per bin, the log2 of the ratio or the difference. This tool
+can normalize the number of reads on each BAM file using the SES method
+proposed by Diaz et al. (2012). "Normalization, bias correction, and peak
+calling for ChIP-seq". Statistical applications in genetics and molecular
+biology, 11(3). Normalization based on read counts is also available. The
+output is either a bedgraph or a bigwig file containing the bin location and
+the resulting comparison values. By default if reads are mated the fragment
+length reported in the BAM file is used.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a bamCorrelate.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/bamCorrelate.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,167 @@
+
+ corrlates pairs of bam files
+
+ deepTools
+
+
+ #set files=[]
+ #set labels=[]
+ #for $i in $inputs
+ #set $files += [str($i.bamfile)]
+ #if str($i.label.value) != "":
+ #set $labels += ["\"%s\"" % ($i.label.value)]
+ #else
+ #set $labels += ["\"%s\"" % ($i.bamfile.name)]
+ #end if
+ #end for
+ bamCorrelate
+ --bamfiles #echo " ".join($files)
+ --labels #echo " ".join($labels)
+
+ --fragmentLength $fragmentLength
+ --corMethod $corMethod
+
+ #set newoutFileName=str($outFileName)+".png"
+ --plotFile $newoutFileName
+
+ #if $outputOpt.showOutputOpt == "yes"
+ #if $outputOpt.outFileRawCounts:
+ --outRawCounts '$outputOpt.outFileRawCounts'
+ #end if
+ #if $outputOpt.outFileCorMatrix:
+ --outFileCorMatrix '$outputOpt.outFileCorMatrix'
+ #end if
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if $advancedOpt.smoothLength:
+ --smoothLength '$advancedOpt.smoothLength'
+ #end if
+
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+
+ --binSize '$advancedOpt.binSize'
+ --numberOfSamples '$advancedOpt.numberOfSamples'
+
+ $advancedOpt.doNotExtendPairedEnds
+ $advancedOpt.ignoreDuplicates
+ $advancedOpt.includeZeros
+
+ #if $advancedOpt.minMappingQuality:
+ --minMappingQuality '$advancedOpt.minMappingQuality'
+ #end if
+ #end if
+
+ --numberOfProcessors 4; mv $newoutFileName $outFileName
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ (outputOpt['showOutputOpt'] == 'yes' and outputOpt['saveRawCounts'] == True)
+
+
+ (outputOpt['showOutputOpt'] == 'yes' and outputOpt['saveCorMatrix'] == True)
+
+
+
+
+**What it does**
+
+Genomes are split into bins of given length. For each bin the number of reads
+found for each of the bam files is counted. A correlation is computed for all
+pairs of bam files.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a bamCoverage.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/bamCoverage.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,143 @@
+
+ Given a BAM file, generates a coverage bigwig file. Multiple options available to count reads and normalize coverage.
+
+ deepTools
+ ucsc_tools
+ numpy
+
+ bamCoverage
+ --bam '$bamInput'
+ --bamIndex ${bamInput.metadata.bam_index}
+ --outFileName '$outFileName'
+ --outFileFormat '$outFileFormat'
+
+ --fragmentLength $fragmentLength
+ --binSize $binSize
+
+
+ #if $scaling.type=='rpkm':
+ --normalizeUsingRPKM
+ #elif $scaling.type=='1x':
+ --normalizeTo1x $scaling.normalizeTo1x
+ #elif $scaling.type=='own':
+ --scaleFactor $scaling.scaleFactor
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if $advancedOpt.smoothLength:
+ --smoothLength '$advancedOpt.smoothLength'
+ #end if
+
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+ $advancedOpt.doNotExtendPairedEnds
+ $advancedOpt.ignoreDuplicates
+
+ #if $advancedOpt.minMappingQuality:
+ --minMappingQuality '$advancedOpt.minMappingQuality'
+ #end if
+
+ #end if
+ --numberOfProcessors 4
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+Given a BAM file, this tool generates a bigWig or bedGraph file of fragment or read coverages. The way the method works is by first calculating all the number of reads (either extended to match the fragment length or not) that overlap each bin in the genome. Bins with zero counts are skipped, i.e. not added to the output file. The resulting read counts can be normalized using either a given scaling factor, the RPKM formula or to get a 1x depth of coverage (RPGC).
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a bamFingerprint.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/bamFingerprint.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,151 @@
+
+ plots profiles of bam files
+
+ deepTools
+
+
+ #set files=[]
+ #set labels=[]
+ #for $i in $inputs
+ #set $files += [str($i.bamfile)]
+ #if str($i.label.value) != "":
+ #set $labels += ["\"%s\"" % ($i.label.value)]
+ #else
+ #set $labels += ["\"%s\"" % ($i.bamfile.name)]
+ #end if
+ #end for
+ bamFingerprint
+ --bamfiles #echo " ".join($files)
+ --labels #echo " ".join($labels)
+
+ --fragmentLength $fragmentLength
+
+ #set newoutFileName=str($outFileName)+".png"
+ --plotFile $newoutFileName
+
+ #if $outputOpt.showOutputOpt == "yes"
+ #if $outputOpt.saveRawCounts:
+ --outRawCounts '$outFileRawCounts'
+ #end if
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if $advancedOpt.smoothLength:
+ --smoothLength '$advancedOpt.smoothLength'
+ #end if
+
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+
+ --binSize '$advancedOpt.binSize'
+ --numberOfSamples '$advancedOpt.numberOfSamples'
+
+ $advancedOpt.doNotExtendPairedEnds
+ $advancedOpt.ignoreDuplicates
+ $advancedOpt.skipZeros
+
+ #if $advancedOpt.minMappingQuality:
+ --minMappingQuality '$advancedOpt.minMappingQuality'
+ #end if
+ #end if
+
+ --numberOfProcessors 4; mv $newoutFileName $outFileName
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ (outputOpt['showOutputOpt'] == 'yes' and outputOpt['saveRawCounts'] == True)
+
+
+
+
+**What it does**
+
+Samples indexed bam files and plots a profile for each bam file. At each
+sample position all reads overlaping a window (bin) of specified length are
+counted. This counts are then sorted and the cumulative sum plotted
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a bigwigCompare.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/bigwigCompare.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,116 @@
+
+ compares two bigwig files based on the number of mapped reads
+
+ deepTools
+ ucsc_tools
+ numpy
+
+
+ bigwigCompare
+ --bigwig1 '$bigwigFile1'
+ --bigwig2 '$bigwigFile2'
+
+ --outFileName '$outFileName'
+ --outFileFormat '$outFileFormat'
+
+ --ratio $comparison_type
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+
+ --missingDataAsZero $advancedOpt.missingDataAsZero
+ --scaleFactors '$advancedOpt.scaleFactor1:$advancedOpt.scaleFactor2'
+ --pseudocount '$advancedOpt.pseudocount'
+ --binSize $advancedOpt.binSize
+
+ #end if
+ --numberOfProcessors 4
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+This tool compares two bigwig files based on the number of mapped reads. To
+compare the bigwig files the genome is partitioned into bins of equal size,
+then the number of reads found in each BAM file are counted for such bins and
+finally a summarizing value is reported. This vaule can be the ratio of the
+number of readsper bin, the log2 of the ratio, the sum or the difference.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a computeGCBias.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/computeGCBias.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,147 @@
+
+
+
+ deepTools
+
+
+
+
+
+ computeGCBias
+ --bamfile '$bamInput'
+ --species '$species'
+ --GCbiasFrequenciesFile $outFileName
+ --fragmentLength $fragmentLength
+
+ #if $source.ref_source=="history":
+ --genome $source.input1
+ #else:
+ --genome "${source.input1_2bit.fields.path}"
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+
+ --binSize '$advancedOpt.binSize'
+ --sampleSize '$advancedOpt.sampleSize'
+ --regionSize '$advancedOpt.regionSize'
+
+ #if $advancedOpt.filterOut:
+ --filterOut $advancedOpt.filterOut
+ #end if
+
+ #if $advancedOpt.extraSampling:
+ --extraSampling $advancedOpt.extraSampling
+ #end if
+
+ #end if
+
+ #set move=""
+ #if $output.showOutputSettings == "yes"
+ #if $output.saveBiasPlot:
+ --biasPlot biasPlot.png
+ #set move="mv biasPlot.png $biasPlot"
+ #end if
+ #end if
+ ; $move
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveBiasPlot'] == True)
+
+
+
+
+**What it does**
+
+Computes the GC bias ussing Benjamini's method [citation]. The resulting GC
+bias can later be used to plot the bias or to correct the bias.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a computeMatrix.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/computeMatrix.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,208 @@
+
+ summarizes and prepares an intermediary file containing scores associated with genomic regions that can be used afterwards to plot a heatmap or a profile
+
+ deepTools
+
+
+ computeMatrix
+ $mode.mode_select
+ --regionsFileName '$regionsFile'
+ --scoreFileName '$scoreFile'
+ --outFileName '$outFileName'
+
+ #if $output.showOutputSettings == "yes"
+ #if $output.saveData:
+ --outFileNameData '$outFileNameData'
+ #end if
+ #if $output.saveMatrix:
+ --outFileNameMatrix '$outFileNameMatrix'
+ #end if
+
+ #if $output.saveSortedRegions:
+ --outFileSortedRegions '$outFileSortedRegions'
+ #end if
+ #end if
+
+ #if $mode.mode_select == "reference-point":
+ --referencePoint $mode.referencePoint
+ $mode.nanAfterEnd
+ --beforeRegionStartLength $mode.beforeRegionStartLength
+ --afterRegionStartLength $mode.afterRegionStartLength
+ #else
+ --regionBodyLength $mode.regionBodyLength
+ --startLabel $mode.startLabel
+ --endLabel $mode.endLabel
+ #if $mode.regionStartLength.regionStartLength_select == "yes":
+ --beforeRegionStartLength $mode.regionStartLength.beforeRegionStartLength
+ --afterRegionStartLength $mode.regionStartLength.afterRegionStartLength
+ #end if
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ --sortRegions '$advancedOpt.sortRegions'
+ --sortUsing '$advancedOpt.sortUsing'
+ --averageTypeBins '$advancedOpt.averageTypeBins'
+ $advancedOpt.missingDataAsZero
+ $advancedOpt.skipZeros
+ -bs $advancedOpt.binSize
+
+ #if $advancedOpt.minThreshold:
+ --minThreshold $advancedOpt.minThreshold
+ #end if
+ #if $advancedOpt.maxThreshold:
+ --maxThreshold $advancedOpt.maxThreshold
+ #end if
+ #if $advancedOpt.scale:
+ --scale $advancedOpt.scale
+ #end if
+
+ #end if
+ --numberOfProcessors 4
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveData'] == True)
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveMatrix'] == True)
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveSortedRegions'] == True)
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+This tool summarizes and prepares an intermediary file containing scores associated with genomic regions that can be used afterwards to plot a heatmap or a profile. Typically, these genomic regions are genes, but any other regions defined in a BED or GFF format can be used. This tool can also be used to filter and sort regions according to their score.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
diff -r 000000000000 -r 09b65c12a75a correctGCBias.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/correctGCBias.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,112 @@
+
+
+
+
+ deepTools
+ ucsc_tools
+
+
+ correctGCBias
+ --bamfile '$bamInput'
+ --species '$species'
+ --GCbiasFrequenciesFile $GCbiasFrequenciesFile
+
+ #if $source.ref_source=="history":
+ --genome $source.input1
+ #else:
+ --genome "${source.input1_2bit.fields.path}"
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes":
+ #if str($advancedOpt.region.value) != '':
+ --region '$advancedOpt.region'
+ #end if
+
+ --binSize '$advancedOpt.binSize'
+ #end if
+
+ #set newoutFileName="corrected."+str($outFileFormat)
+
+ --correctedFile $newoutFileName; mv $newoutFileName $outFileName
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+Computes the GC bias ussing Benjamini's method [citation]. The resulting GC
+bias can later be used to plot the bias or to correct the bias.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a heatmapper.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/heatmapper.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,370 @@
+
+ creates a heatmap for a score associated to genomic regions
+
+
+ samtools
+ numpy
+ matplotlib
+ scipy
+ ucsc_tools
+ deepTools
+
+
+
+ heatmapper
+ --matrixFile $matrixFile
+ #if $output.showOutputSettings == "yes"
+ #set newoutFileName=str($outFileName)+"."+str($output.outFileFormat)
+ --outFileName $newoutFileName
+ #if $outFileNameData:
+ --outFileNameData '$outFileNameData'
+ #end if
+
+ #if $outFileNameMatrix:
+ --outFileNameMatrix '$outFileNameMatrix'
+ #end if
+
+ #if $outFileSortedRegions:
+ --outFileSortedRegions '$outFileSortedRegions'
+ #end if
+ #else
+ #set newoutFileName=str($outFileName)+".png"
+ --outFileName $newoutFileName
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes"
+ #if $advancedOpt.sortRegions:
+ --sortRegions '$advancedOpt.sortRegions'
+ #end if
+
+ #if $advancedOpt.sortUsing:
+ --sortUsing '$advancedOpt.sortUsing'
+ #end if
+
+ #if $advancedOpt.averageTypeSummaryPlot:
+ --averageTypeSummaryPlot '$advancedOpt.averageTypeSummaryPlot'
+ #end if
+
+ #if str($advancedOpt.missingDataColor.value) != "None":
+ --missingDataColor '$advancedOpt.missingDataColor'
+ #end if
+
+ --colorMap '$advancedOpt.colorMap'
+
+ #if $advancedOpt.zMin:
+ --zMin $advancedOpt.zMin
+ #end if
+ #if $advancedOpt.zMax:
+ --zMax $advancedOpt.zMax
+ #end if
+
+ #if $advancedOpt.yMin:
+ --yMin $advancedOpt.yMin
+ #end if
+ #if $advancedOpt.yMax:
+ --yMax $advancedOpt.yMax
+ #end if
+
+ --xAxisLabel '$advancedOpt.xAxisLabel'
+ --yAxisLabel '$advancedOpt.yAxisLabel'
+
+ --heatmapWidth $advancedOpt.heatmapWidth
+ --heatmapHeight $advancedOpt.heatmapHeight
+
+ --whatToShow '$advancedOpt.whatToShow'
+
+ --startLabel '$advancedOpt.startLabel'
+ --endLabel '$advancedOpt.endLabel'
+ --refPointLabel '$advancedOpt.referencePointLabel'
+ --regionsLabel '$advancedOpt.regionsLabel'
+
+ #if str($advancedOpt.plotTitle.value) != "None":
+ --plotTitle '$advancedOpt.plotTitle'
+ #end if
+
+ $advancedOpt.onePlotPerGroup
+ #end if
+ ; mv $newoutFileName $outFileName
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveData'] == True)
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveMatrix'] == True)
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveSortedRegions'] == True)
+
+
+
+
+
+
+
+
+
+
+**What it does**
+
+HeatMapper visualizes scores associated with genomic regions, for example log2 fold change values obtained from ChIP-seq experiments. Those values can be visualized individually along each of the regions provided by the user.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
diff -r 000000000000 -r 09b65c12a75a profiler.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/profiler.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,189 @@
+
+
+ creates a profile plot for a score associated to genomic regions
+
+
+ deepTools
+
+
+ profiler
+ --matrixFile $matrixFile
+
+ #if $output.showOutputSettings == "yes"
+ #set newoutFileName=str($outFileName)+"."+str($output.outFileFormat)
+ --outFileName $newoutFilename
+ #if $output.outFileNameData:
+ --outFileNameData '$output.outFileNameData'
+ #end if
+
+ #if $output.outFileNameMatrix:
+ --outFileNameMatrix '$output.outFileNameMatrix'
+ #end if
+
+ #if $output.outFileSortedRegions:
+ --outFileSortedRegions '$output.outFileSortedRegions'
+ #end if
+ #else
+ #set newoutFileName=str($outFileName)+".png"
+ --outFileName $newoutFileName
+ #end if
+
+ #if $scaleRegions.showScaleRegionsOpt == "yes":
+ --startLabel $scaleRegions.startLabel
+ --endLabel $scaleRegions.endLabel
+ --refPointLabel $scaleRegions.refPointLabel
+ #end if
+
+ #if $advancedOpt.showAdvancedOpt == "yes"
+ #if $advancedOpt.averageType:
+ --averageType '$advancedOpt.averageType'
+ #end if
+ --plotHeight $advancedOpt.plotHeight
+ --plotWidth $advancedOpt.plotWidth
+ --plotType $advancedOpt.plotType
+
+ --regionsLabel '$advancedOpt.regionsLabel'
+
+ #if str($advancedOpt.plotTitle.value) != "None":
+ --plotTitle '$advancedOpt.plotTitle'
+ #end if
+
+ $advancedOpt.onePlotPerGroup
+
+ #if $advancedOpt.yMin:
+ --yMin $advancedOpt.yMin
+ #end if
+ #if $advancedOpt.yMax:
+ --yMax $advancedOpt.yMax
+ #end if
+
+ --xAxisLabel '$advancedOpt.xAxisLabel'
+ #if str($advancedOpt.yAxisLabel.value) != "None":
+ --yAxisLabel '$advancedOpt.yAxisLabel'
+ #end if
+ #end if
+ ; mv $newoutFileName $outFileName
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveData'] == True)
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveMatrix'] == True)
+
+
+ (output['showOutputSettings'] == 'yes' and output['saveSortedRegions'] == True)
+
+
+
+
+**What it does**
+
+This tool creates a profile plot for a score associated to genomic regions.
+Typically, these regions are genes, but any other regions defined in a BED or
+GFF format will work. A preprocessed matrix generated by the tool
+computeMatrix is required.
+
+-----
+
+.. class:: infomark
+
+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`_.
+
+This tool is developed by the `Bioinformatics and Deep-Sequencing Unit`_ at the `Max Planck Institute for Immunobiology and Epigenetics`_.
+
+
+.. _Bioinformatics and Deep-Sequencing Unit: http://www3.ie-freiburg.mpg.de/facilities/research-facilities/bioinformatics-and-deep-sequencing-unit/
+.. _Max Planck Institute for Immunobiology and Epigenetics: http://www3.ie-freiburg.mpg.de
+.. _Fidel Ramirez: ramirez@ie-freiburg.mpg.de
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a test-data/master.mat.gz
Binary file test-data/master.mat.gz has changed
diff -r 000000000000 -r 09b65c12a75a test-data/master.png
Binary file test-data/master.png has changed
diff -r 000000000000 -r 09b65c12a75a test-data/test.bw
Binary file test-data/test.bw has changed
diff -r 000000000000 -r 09b65c12a75a test-data/test2.bed
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test2.bed Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,8 @@
+ch1 100 150 CG11023 0 +
+ch2 150 175 cda5 0 -
+ch3 100 125 cda8 0 +
+#Group 1
+ch1 75 125 C11023 0 +
+ch2 125 150 ca5 0 -
+ch3 75 100 ca8 0 +
+#Group 2
diff -r 000000000000 -r 09b65c12a75a tool-data/deepTools_seqs.loc.sample
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool-data/deepTools_seqs.loc.sample Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,30 @@
+#This is a sample file distributed with Galaxy that enables tools
+#to use a directory of 2bit genome files for use with deepTools. You will
+#need to supply these files and then create a deepTools_seqs.loc file
+#similar to this one (store it in this directory) that points to
+#the directories in which those files are stored. The deepTools_seqs.loc
+#file has this format (white space characters are TAB characters):
+#
+#
+#
+#So, for example, if your deepTools_seqs.loc began like this:
+#
+#hg18Human (Homo sapiens): hg18/depot/data2/galaxy/twobit/hg18.2bit
+#hg19Human (Homo sapiens): hg19/depot/data2/galaxy/twobit/hg19.2bit
+#mm9Mouse (Mus musculus): mm9/depot/data2/galaxy/twobit/mm9.2bit
+#
+#then your /depot/data2/galaxy/twobit/ directory
+#would need to contain the following 2bit files:
+#
+#-rw-r--r-- 1 james universe 830134 2005-09-13 10:12 hg18.2bit
+#-rw-r--r-- 1 james universe 527388 2005-09-13 10:12 hg19.2bit
+#-rw-r--r-- 1 james universe 269808 2005-09-13 10:12 mm9.2bit
+#
+#Your deepTools_seqs.loc file should include an entry per line for
+#each file you have stored that you want to be available. Note that
+#your files should all have the extension '2bit'.
+#
+#Note that for backwards compatibility with workflows, the unique ID of
+#an entry must be the path that was in the original loc file, because that
+#is the value stored in the workflow for that parameter.
+#
diff -r 000000000000 -r 09b65c12a75a tool_data_table_conf.xml.sample
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_data_table_conf.xml.sample Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,7 @@
+
+
+
+
diff -r 000000000000 -r 09b65c12a75a tool_dependencies.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_dependencies.xml Sat Aug 03 13:42:12 2013 -0400
@@ -0,0 +1,67 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bedGraphToBigWig
+ http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.i386/bedGraphToBigWig
+ http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.i386/bedGraphToBigWig
+ http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.x86_64/bedGraphToBigWig
+
+
+ $INSTALL_DIR/bedGraphToBigWig
+
+
+ http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/bigWigInfo
+ http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.i386/bigWigInfo
+ http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.i386/bigWigInfo
+ http://hgdownload.cse.ucsc.edu/admin/exe/macOSX.x86_64/bigWigInfo
+
+
+ $INSTALL_DIR/bigWigInfo
+
+
+ $INSTALL_DIR
+
+
+
+ The tools downloaded by this dependency definition are free for academic use. TODO: UCSC tools are only available with their latest version. That is not good for reproducibility.
+
+
+
+
+
+
+
+ git clone --recursive https://github.com/fidelram/deepTools.git
+ git reset --hard 98e5d8a61431ea8605c0643d991a1a5d8999b4dc
+ $INSTALL_DIR/lib/python
+
+ export PYTHONPATH=$PYTHONPATH:$INSTALL_DIR/lib/python &&
+ python setup.py install --install-lib $INSTALL_DIR/lib/python --install-scripts $INSTALL_DIR/bin
+
+
+
+ $INSTALL_DIR/bin
+ $INSTALL_DIR/lib/python
+
+
+
+ Installation of deepTools from Fidel Ramirez.
+
+