From 5d18d6297a1eac7270e94efce2269ec0add76363 Mon Sep 17 00:00:00 2001
From: Wai Yi Leung <w.y.leung@lumc.nl>
Date: Mon, 19 Oct 2015 18:04:45 +0200
Subject: [PATCH] Adjusting the Sambamba view filter, exclude "split reads" in
RNA-seq experiments
---
.../main/scala/nl/lumc/sasc/biopet/pipelines/gears/Gears.scala | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/public/gears/src/main/scala/nl/lumc/sasc/biopet/pipelines/gears/Gears.scala b/public/gears/src/main/scala/nl/lumc/sasc/biopet/pipelines/gears/Gears.scala
index 4ae307a12..f67b1b8fc 100644
--- a/public/gears/src/main/scala/nl/lumc/sasc/biopet/pipelines/gears/Gears.scala
+++ b/public/gears/src/main/scala/nl/lumc/sasc/biopet/pipelines/gears/Gears.scala
@@ -69,7 +69,7 @@ class Gears(val root: Configurable) extends QScript with SummaryQScript {
// sambamba view -f bam -F "unmapped or mate_is_unmapped" <alnFile> > <extracted.bam>
val samFilterUnmapped = new SambambaView(qscript)
samFilterUnmapped.input = bamfile
- samFilterUnmapped.filter = Some("(unmapped or mate_is_unmapped) and not (secondary_alignment)")
+ samFilterUnmapped.filter = Some("(unmapped or mate_is_unmapped) and not (secondary_alignment) and [XH] == null")
samFilterUnmapped.output = new File(outputDir, s"$outputName.unmapped.bam")
samFilterUnmapped.isIntermediate = false
add(samFilterUnmapped)
--
GitLab