diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala index 419c5ce12fdc6b211073bdabf53f7c46faffb5f8..a8b3daf680e02d130ed10c857c27114c918a6c21 100644 --- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala +++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala @@ -8,6 +8,7 @@ import scala.collection.JavaConversions._ import java.io.File import nl.lumc.sasc.biopet.extensions.gatk.{ CombineVariants, CombineGVCFs } import nl.lumc.sasc.biopet.extensions.picard.AddOrReplaceReadGroups +import nl.lumc.sasc.biopet.pipelines.bammetrics.BamMetrics import nl.lumc.sasc.biopet.pipelines.mapping.Mapping import org.broadinstitute.gatk.queue.QScript import org.broadinstitute.gatk.utils.commandline.{ Argument } @@ -195,6 +196,7 @@ class GatkPipeline(val root: Configurable) extends QScript with MultiSampleQScri } else throw new IllegalStateException("Readgroup sample and/or library of input bamfile is not correct, file: " + bamFile + "\nPossible to set 'correct_readgroups' to true on config to automatic fix this") } + addAll(BamMetrics(this, bamFile, runDir + "metrics/").functions) libraryOutput.mappedBamFile = bamFile } else logger.error("Sample: " + sampleID + ": No R1 found for run: " + runConfig)