Mercurial > repos > bgruening > deeptools
diff computeGCBias.xml @ 5:c54d31467be4 draft
Uploaded
author | bgruening |
---|---|
date | Fri, 29 Nov 2013 08:59:04 -0500 |
parents | d957e25e18a3 |
children | c53a73b8eef9 |
line wrap: on
line diff
--- a/computeGCBias.xml Sat Nov 16 07:54:13 2013 -0500 +++ b/computeGCBias.xml Fri Nov 29 08:59:04 2013 -0500 @@ -8,63 +8,52 @@ <import>deepTools_macros.xml</import> </macros> <command> - #import tempfile - #set $temp_dir = os.path.abspath(tempfile.mkdtemp()) + ln -s $bamInput local_bamInput.bam; + ln -s $bamInput.metadata.bam_index local_bamInput.bam.bai; - #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)) + computeGCBias - computeGCBias - - @THREADS@ + @THREADS@ - --bamfile '$temp_bam_path' - --GCbiasFrequenciesFile $outFileName - --fragmentLength $fragmentLength + --bamfile 'local_bamInput.bam' + --GCbiasFrequenciesFile $outFileName + --fragmentLength $fragmentLength - @reference_genome_source@ - + @reference_genome_source@ - #if $effectiveGenomeSize.effectiveGenomeSize_opt == "specific": - --effectiveGenomeSize $effectiveGenomeSize.effectiveGenomeSize - #else: - --effectiveGenomeSize $effectiveGenomeSize.effectiveGenomeSize_opt - #end if - + #if $effectiveGenomeSize.effectiveGenomeSize_opt == "specific": + --effectiveGenomeSize $effectiveGenomeSize.effectiveGenomeSize + #else: + --effectiveGenomeSize $effectiveGenomeSize.effectiveGenomeSize_opt + #end if - #if $advancedOpt.showAdvancedOpt == "yes": - #if str($advancedOpt.region.value) != '': - --region '$advancedOpt.region' - #end if - - --sampleSize '$advancedOpt.sampleSize' - --regionSize '$advancedOpt.regionSize' + #if $advancedOpt.showAdvancedOpt == "yes": + #if str($advancedOpt.region.value) != '': + --region '$advancedOpt.region' + #end if - #if $advancedOpt.filterOut: - --filterOut $advancedOpt.filterOut - #end if + --sampleSize '$advancedOpt.sampleSize' + --regionSize '$advancedOpt.regionSize' + + #if $advancedOpt.filterOut: + --filterOut $advancedOpt.filterOut + #end if - #if $advancedOpt.extraSampling: - --extraSampling $advancedOpt.extraSampling - #end if - - #end if + #if $advancedOpt.extraSampling: + --extraSampling $advancedOpt.extraSampling + #end if + #end if - #if $saveBiasPlot: - --biasPlot $biasPlot - #end if + #if $saveBiasPlot: + --biasPlot $biasPlot + #end if -## #if $output.showOutputSettings == "yes" -## #if $output.saveBiasPlot: -## --biasPlot biasPlot.png ; -## mv biasPlot.png $biasPlot -## #end if -## #end if - - ; rm $temp_dir -rf + ## #if $output.showOutputSettings == "yes" + ## #if $output.saveBiasPlot: + ## --biasPlot biasPlot.png ; + ## mv biasPlot.png $biasPlot + ## #end if + ## #end if </command> <inputs>