Mercurial > repos > iuc > bedtools

changeset 24:33c3ddea63c5 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bedtools commit 09721a12bdd393441b162056316241aba3895a4f
author iuc
date Sun, 18 Feb 2018 15:07:08 -0500
parents 13400f3c3ec5
children e36671a0f646
files coverageBed.xml
diffstat 1 files changed, 9 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/coverageBed.xml	Fri Feb 09 09:00:06 2018 -0500
+++ b/coverageBed.xml	Sun Feb 18 15:07:08 2018 -0500
@@ -1,4 +1,4 @@
-<tool id="bedtools_coveragebed" name="Compute both the depth and breadth of coverage" version="@WRAPPER_VERSION@.1">
+<tool id="bedtools_coveragebed" name="Compute both the depth and breadth of coverage" version="@WRAPPER_VERSION@.2">
     <description>of features in file B on the features in file A (bedtools coverage)</description>
     <macros>
         <import>macros.xml</import>
@@ -23,11 +23,18 @@
         -a '$inputA'
         #if str($reduce_or_iterate.reduce_or_iterate_selector) == 'iterate':
             -b '$reduce_or_iterate.inputB'
+            #if $reduce_or_iterate.inputB.is_of_type('bam'):
+                -sorted  -g <(samtools view -H $reduce_or_iterate.inputB | tr ':' '\t' | grep SN | cut -f 3,5)
+            #end if
         #else:
             #set inputBs = "' '".join([str($file) for $file in $reduce_or_iterate.inputB])
             -b '$inputBs'
         #end if
-        | sort -k1,1 -k2,2n
+        #if $inputA.is_of_type('gff') or $inputA.is_of_type('gff3'):
+            | sort -k1,1 -k4,2n
+        #else:
+            | sort -k1,1 -k2,2n
+        #end if
         > '$output'
 ]]>
     </command>