diff --git a/public/shiva/src/main/scala/nl/lumc/sasc/biopet/pipelines/shiva/ShivaTrait.scala b/public/shiva/src/main/scala/nl/lumc/sasc/biopet/pipelines/shiva/ShivaTrait.scala
index 2b13c8d41f42589c190edd500d9f49118ed7dd2b..19c11ca883cd6b8597ec41e3be1ed4bdb0e33c06 100644
--- a/public/shiva/src/main/scala/nl/lumc/sasc/biopet/pipelines/shiva/ShivaTrait.scala
+++ b/public/shiva/src/main/scala/nl/lumc/sasc/biopet/pipelines/shiva/ShivaTrait.scala
@@ -197,6 +197,16 @@ trait ShivaTrait extends MultiSampleQScript with Reference {
                   val bamLn = Ln(qscript, oldBamFile, bamFile.get)
                   bamLn.deps :+= baiLn.output
                   add(bamLn)
+
+                  val bamMetrics = new BamMetrics(qscript)
+                  bamMetrics.sampleId = Some(sampleId)
+                  bamMetrics.libId = Some(libId)
+                  bamMetrics.inputBam = bamFile.get
+                  bamMetrics.outputDir = new File(libDir, "metrics")
+                  bamMetrics.init()
+                  bamMetrics.biopetScript()
+                  addAll(bamMetrics.functions)
+                  addSummaryQScript(bamMetrics)
                 }
             }
           }