diff bam_to_bigwig.xml @ 0:dd444ce346ae draft default tip

planemo upload commit 54155cb43a07d563ee6861c2410a167b9c879ff4-dirty
author rmarenco
date Thu, 14 Jul 2016 18:55:27 -0400
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/bam_to_bigwig.xml	Thu Jul 14 18:55:27 2016 -0400
@@ -0,0 +1,30 @@
+<tool id="bam_to_bigwig" name="Convert Bam to BigWig" version="0.1.0">
+    <requirements>
+        <requirement type="package" version="312">ucsc_tools</requirement>
+        <requirement type="package" version="2.24">bedtools</requirement>
+    </requirements>
+    <stdio>
+        <exit_code range="1:" />
+    </stdio>
+    <command><![CDATA[
+        faToTwoBit $input_reference_fasta two_bit.2bit;
+        twoBitInfo two_bit.2bit stdout | sort -k2nr > chrom.sizes;
+        bedtools genomecov -bg -split -ibam $input_bam -g chrom.sizes > temp.bg;
+        bedGraphToBigWig temp.bg chrom.sizes $output_bigwig;
+    ]]></command>
+    <inputs>
+        <param format="bam" name="input_bam" type="data" label="Choose BAM file" />
+        <param format="fasta" name="input_reference_fasta"
+               type="data" label="Choose Reference Genome (Fasta)" />
+    </inputs>
+    <outputs>
+        <data format="bigwig" name="output_bigwig" />
+    </outputs>
+    <help>
+        This tool converts a Bam file to a BigWig file.
+
+        It needs:
+            - A Bam file
+            - The reference genome as a Fasta file
+    </help>
+</tool>
\ No newline at end of file