diff multiCov.xml @ 8:0d3aa592ce27 draft

Uploaded
author iuc
date Tue, 28 Apr 2015 22:56:34 -0400
parents
children a2d4c30ba2f9
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/multiCov.xml	Tue Apr 28 22:56:34 2015 -0400
@@ -0,0 +1,75 @@
+<tool id="bedtools_multicovtbed" name="MultiCovBed" version="@WRAPPER_VERSION@.0">
+    <description></description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements" />
+    <expand macro="stdio" />
+    <command>
+<![CDATA[
+        #for $i, $bam in enumerate( $bams ):
+            ln -s -f $bam ${i}.bam &&
+            ln -s -f $bam.metadata.bam_index ${i}.bam.bai &&
+        #end for
+
+        bedtools multicov
+            -bed $input
+            -bams
+            #for $i, $bam in enumerate( $bams ):
+                ${i}.bam
+            #end for
+            $strand
+            -f $overlap
+            $reciprocal
+            $split
+            -q $q
+            $duplicate
+            $failed
+            $proper
+            > $output
+]]>
+    </command>
+    <inputs>
+        <param name="input" format="bed" type="data" label="Sorted BED file" />
+        <!-- Additional files, if the user needs more -->
+        <param name="bams" format="bam" type="data" multiple="True" label="BAM file" />
+
+        <expand macro="strand2" />
+        <expand macro="overlap" />
+        <expand macro="reciprocal" />
+        <expand macro="split" />
+        <param name="q" type="integer" value="0" label="Minimum mapping quality (MAPQ) allowed" help="(-q)" />
+        <param name="duplicate" type="boolean" checked="False" truevalue="-D" falsevalue=""
+            label="Include duplicate reads"
+            help="Default counts non-duplicates only. (-D)" />
+        <param name="failed" type="boolean" checked="false" truevalue="-F" falsevalue=""
+            label="Include failed-QC reads" 
+            help="Default counts pass-QC reads only (-F)"/>
+        <param name="proper" type="boolean" checked="false" truevalue="-p" falsevalue=""
+            label="Only count proper pairs"
+            help="Default counts all alignments with MAPQ > -q argument, regardless of the BAM FLAG field. (-p)" />
+    </inputs>
+    <outputs>
+        <data name="output" metadata_source="input" format_source="input" />
+    </outputs>
+    <tests>
+        <test>
+            <param name="input" value="multiCov1.bed" ftype="bed" />
+            <param name="bams" value="srma_in3.bam,srma_in3.bam" ftpye="bam"/>
+            <param name="q" value="1"/>
+            <param name="split" value="False"/>
+            <output name="output" file="multiCovBed_result1.bed" ftype="bed" />
+        </test>
+    </tests>
+    <help>
+<![CDATA[
+**What it does**
+
+bedtools multicov, reports the count of alignments from multiple position-sorted and indexed BAM files that overlap intervals in a BED file. Specifically, for each BED interval provided, it reports a separate count of overlapping alignments from each BAM file.
+
+
+@REFERENCES@
+]]>
+    </help>
+    <expand macro="citations" />
+</tool>