Mercurial > repos > bgruening > deeptools
diff correctGCBias.xml @ 26:2ad3b027dfcd draft
Uploaded
author | bgruening |
---|---|
date | Thu, 06 Feb 2014 06:40:05 -0500 |
parents | |
children | bf1b1dcdd67b |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/correctGCBias.xml Thu Feb 06 06:40:05 2014 -0500 @@ -0,0 +1,108 @@ +<tool id="deeptools_correctGCBias" name="correctGCBias" version="1.0.4"> + <description>uses the output from computeGCBias to generate corrected BAM files</description> + <expand macro="requirements" /> + <expand macro="stdio" /> + <macros> + <token name="@BINARY@">correctGCBias</token> + <import>deepTools_macros.xml</import> + </macros> + <command> + #import tempfile + #set $temp_dir = os.path.abspath(tempfile.mkdtemp()) + + #set $temp_bam_handle = tempfile.NamedTemporaryFile( dir=$temp_dir ) + #set $temp_bam_path = $temp_bam_handle.name + '.bam' + #silent $temp_bam_handle.close() + #silent os.system("ln -s %s %s" % (str($bamInput), $temp_bam_path)) + #silent os.system("ln -s %s %s.bai" % (str($bamInput.metadata.bam_index), $temp_bam_path)) + + + correctGCBias + + @THREADS@ + + --bamfile '$temp_bam_path' + --GCbiasFrequenciesFile $GCbiasFrequenciesFile + + @reference_genome_source@ + + + #if $effectiveGenomeSize.effectiveGenomeSize_opt == "specific": + --effectiveGenomeSize $effectiveGenomeSize.effectiveGenomeSize + #else: + --effectiveGenomeSize $effectiveGenomeSize.effectiveGenomeSize_opt + #end if + + #if str($region).strip() != '': + --region '$region' + #end if + + #if $advancedOpt.showAdvancedOpt == "yes": + --binSize '$advancedOpt.binSize' + #end if + + ###set newoutFileName="corrected."+str($outFileFormat) + ##--correctedFile $newoutFileName; + --correctedFile "corrected.bam"; + + ##mv $newoutFileName $outFileName + mv "corrected.bam" $outFileName + </command> + <inputs> + <param name="GCbiasFrequenciesFile" type="data" format="tabular" label="Output of computeGCBias" /> + <param name="bamInput" format="bam" type="data" label="BAM file" help="This should be same file that was used for computeGCbias. The BAM file must be sorted."/> + <expand macro="reference_genome_source" /> + <expand macro="effectiveGenomeSize" /> + + <!-- + <param name="outFileFormat" type="select" label="File format of the output"> + <option value="bam">bam</option> + <option value="bw">bigwig</option> + <option value="bg">bedgraph</option> + </param> + --> + <expand macro="region_limit_operation" /> + + <conditional name="advancedOpt"> + <param name="showAdvancedOpt" type="select" label="Show advanced options" > + <option value="no" selected="true">no</option> + <option value="yes">yes</option> + </param> + <when value="no" /> + <when value="yes"> + <param name="binSize" type="integer" value="50" min="1" + label="Bin size in bp" + help="Size of the bins in bp for the output of the bigwig/bedgraph file."/> + </when> + </conditional> + </inputs> + <outputs> + <data format="bam" name="outFileName"> + <!--<change_format> + <when input="outFileFormat" value="bw" format="bigwig" /> + <when input="outFileFormat" value="bam" format="bam" /> + <when input="outFileFormat" value="bg" format="bedgraph" /> + </change_format>--> + </data> + </outputs> + <help> + +**What it does** + +This tool requires the output from computeGCBias to correct a given BAM file according to the method proposed by +Benjamini and Speed (2012) Nucleic Acids Res. +The resulting BAM file can be used in any downstream analyses, but be aware that you should not filter out duplicates from here on. + +You can find more details on the correctGCBias wiki page: https://github.com/fidelram/deepTools/wiki/QC#wiki-correctGCbias + + +**Output files**: + +- GC-normalized BAM file + +----- + +@REFERENCES@ + + </help> +</tool>