diff --git a/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala b/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala index 53024816a4e756a6314b8981b77907fb2f6b9067..90557d3e86514ab47380e710385ea5d9fb27e86d 100644 --- a/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala +++ b/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala @@ -155,6 +155,7 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { for ((chunk, fastqfile) <- chunks) { var R1 = fastqfile._1 var R2 = fastqfile._2 + val chunkDir = if (chunking) outputDir + chunk + "/" else outputDir if (!skipFlexiprep) { flexiprep.input_R1 = fastq_R1 if (paired) flexiprep.input_R2 = fastq_R2 @@ -162,14 +163,14 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { flexiprep.init flexiprep.runInitialJobs //flexiprep.biopetScript - val flexiout = flexiprep.runTrimClip(R1, R2, outputDir + chunk, chunk) + val flexiout = flexiprep.runTrimClip(R1, R2, chunkDir + "flexiprep/", chunk) logger.debug(chunk + " - " + flexiout) R1 = flexiout._1 if (paired) R2 = flexiout._2 fastq_R1_output :+= R1 fastq_R2_output :+= R2 } - val chunkDir = if (chunking) outputDir + chunk + "/" else outputDir + if (aligner == "bwa") { val bwaCommand = new Bwa(this) bwaCommand.R1 = R1