diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkVariantcalling.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkVariantcalling.scala index b4f8c127c9ff39df392d45bb32ba0556072fbf97..884b2d837d50c7040edd98912e1983703659e0b8 100644 --- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkVariantcalling.scala +++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkVariantcalling.scala @@ -51,7 +51,7 @@ class GatkVariantcalling(val root: Configurable) extends QScript with BiopetQScr else if (!outputDir.endsWith("/")) outputDir += "/" val baseRecalibrator = new BaseRecalibrator(this) - if (useBaseRecalibration && baseRecalibrator.knownSites.isEmpty) { + if (preProcesBams && useBaseRecalibration && baseRecalibrator.knownSites.isEmpty) { logger.warn("No Known site found, skipping base recalibration") useBaseRecalibration = false }