# HG changeset patch
# User jbrayet
# Date 1448439807 18000
# Node ID 96f453b29673af2f7007ad4db3833f88b7c020ef
# Parent  bb21f59022f8e936ac058b7f9f310746360e5d6e
Uploaded

diff -r bb21f59022f8 -r 96f453b29673 samtools_rmdup.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/samtools_rmdup.xml	Wed Nov 25 03:23:27 2015 -0500
@@ -0,0 +1,54 @@
+<tool id="samtools_rmdup" name="RmDup" version="2.0">
+  <description>remove PCR duplicates</description>
+  <macros>
+    <import>macros.xml</import>
+  </macros>
+  <expand macro="requirements"></expand>
+  <expand macro="stdio"></expand>
+  <expand macro="version_command"></expand>
+  <command>samtools rmdup 
+  #if str( $bam_paired_end_type.bam_paired_end_type_selector ) == "PE"
+      ${bam_paired_end_type.force_se}
+  #else:
+      -s
+  #end if
+  "$input1" "$output1"
+  </command>
+  <inputs>
+    <param name="input1" type="data" format="bam" label="BAM File" />
+    
+    <conditional name="bam_paired_end_type">
+      <param name="bam_paired_end_type_selector" type="select" label="Is this paired-end or single end data">
+        <option value="PE" selected="True">BAM is paired-end</option>
+        <option value="SE">BAM is single-end (-s)</option>
+      </param>
+      <when value="PE">
+        <param name="force_se" type="boolean" label="Treat as single-end" help="-S" truevalue="-S" falsevalue="" checked="False"/>
+      </when>
+      <when value="SE" /> <!-- No extra parameters here -->
+    </conditional>
+    
+  </inputs>
+  <outputs>
+    <data name="output1" format="bam" />
+  </outputs>
+  <tests>
+    <test>
+      <param name="input1" value="samtools-rmdup-input1.bam" ftype="bam" />
+      <param name="bam_paired_end_type_selector" value="PE" />
+      <param name="force_se" />
+      <output name="output1" file="samtools-rmdup-test1.bam" ftype="bam" sort="True" />
+    </test>
+  </tests>
+  <help>
+
+**What it does**
+
+Remove potential PCR duplicates: if multiple read pairs have identical external coordinates, only retain the pair with highest mapping quality. In the paired-end mode, this command ONLY works with FR orientation and requires ISIZE is correctly set. It does not work for unpaired reads (e.g. two ends mapped to different chromosomes or orphan reads). This tool has the following parameters::
+
+  -s    rmdup for SE reads
+  -S    treat PE reads as SE in rmdup (force -s)
+
+  </help>
+    <expand macro="citations"></expand>
+</tool>