diff --git a/public/carp/src/main/scala/nl/lumc/sasc/biopet/pipelines/carp/Carp.scala b/public/carp/src/main/scala/nl/lumc/sasc/biopet/pipelines/carp/Carp.scala index fb900410363933159cbf0da593ac9a33013f8026..fc89aada73b69515c6de101c1fa93292160d293d 100644 --- a/public/carp/src/main/scala/nl/lumc/sasc/biopet/pipelines/carp/Carp.scala +++ b/public/carp/src/main/scala/nl/lumc/sasc/biopet/pipelines/carp/Carp.scala @@ -124,11 +124,20 @@ class Carp(val root: Configurable) extends QScript with MultiSampleQScript with val bamMetrics = BamMetrics(qscript, bamFile, new File(sampleDir, "metrics"), sampleId = Some(sampleId)) addAll(bamMetrics.functions) addSummaryQScript(bamMetrics) + + val bamMetricsFilter = BamMetrics(qscript, bamFileFilter, new File(sampleDir, "metrics-filter"), sampleId = Some(sampleId)) + addAll(bamMetricsFilter.functions) + bamMetricsFilter.summaryName = "bamMetrics-filter" + addSummaryQScript(bamMetricsFilter) + addAll(Bam2Wig(qscript, bamFile).functions) + addAll(Bam2Wig(qscript, bamFileFilter).functions) val samtoolsView = new SamtoolsView(qscript) samtoolsView.input = bamFile samtoolsView.output = bamFileFilter + samtoolsView.b = true + samtoolsView.h = true add(samtoolsView) val buildBamIndex = new BuildBamIndex(qscript)