From ad0c55bcfc1a1bdd8b2ac1c7da3792fada932e59 Mon Sep 17 00:00:00 2001 From: bow <bow@bow.web.id> Date: Tue, 3 Feb 2015 16:42:13 +0100 Subject: [PATCH] Rename libraryId -> libId for consistency --- .../sasc/biopet/pipelines/basty/Basty.scala | 2 +- .../biopet/pipelines/gatk/GatkPipeline.scala | 20 +++++++++---------- .../sasc/biopet/core/MultiSampleQScript.scala | 14 ++++++------- .../MultisamplePipelineTemplate.scala | 2 +- .../sasc/biopet/pipelines/carp/Carp.scala | 6 +++--- .../pipelines/flexiprep/Flexiprep.scala | 8 ++++---- .../flexiprep/FlexiprepSummary.scala | 2 +- .../sasc/biopet/pipelines/kopisu/Kopisu.scala | 2 +- .../biopet/pipelines/mapping/Mapping.scala | 14 ++++++------- .../sasc/biopet/pipelines/sage/Sage.scala | 10 +++++----- .../sasc/biopet/pipelines/yamsvp/Yamsvp.scala | 8 ++++---- 11 files changed, 44 insertions(+), 44 deletions(-) diff --git a/protected/basty/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala b/protected/basty/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala index c260f4561..934224a23 100644 --- a/protected/basty/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala +++ b/protected/basty/src/main/scala/nl/lumc/sasc/biopet/pipelines/basty/Basty.scala @@ -33,7 +33,7 @@ class Basty(val root: Configurable) extends QScript with MultiSampleQScript { def makeSample(id: String) = new Sample(id) class Sample(sampleId: String) extends AbstractSample(sampleId) { def makeLibrary(id: String) = new Library(id) - class Library(libraryId: String) extends AbstractLibrary(libraryId) { + class Library(libId: String) extends AbstractLibrary(libId) { protected def addJobs(): Unit = {} } diff --git a/protected/biopet-gatk-pipelines/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala b/protected/biopet-gatk-pipelines/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala index 88600e576..ee814bd31 100644 --- a/protected/biopet-gatk-pipelines/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala +++ b/protected/biopet-gatk-pipelines/src/main/scala/nl/lumc/sasc/biopet/pipelines/gatk/GatkPipeline.scala @@ -43,10 +43,10 @@ class GatkPipeline(val root: Configurable) extends QScript with MultiSampleQScri def makeSample(id: String) = new Sample(id) class Sample(sampleId: String) extends AbstractSample(sampleId) { def makeLibrary(id: String) = new Library(id) - class Library(libraryId: String) extends AbstractLibrary(libraryId) { + class Library(libId: String) extends AbstractLibrary(libId) { val mapping = new Mapping(qscript) mapping.sampleId = sampleId - mapping.libraryId = libraryId + mapping.libId = libId mapping.outputDir = libDir + "/variantcalling/" /** Library variantcalling */ @@ -67,8 +67,8 @@ class GatkPipeline(val root: Configurable) extends QScript with MultiSampleQScri if (!bamFile.exists) throw new IllegalStateException("Bam in config does not exist, file: " + bamFile) if (config("bam_to_fastq", default = false).asBoolean) { - val samToFastq = SamToFastq(qscript, bamFile, libDir + sampleId + "-" + libraryId + ".R1.fastq", - libDir + sampleId + "-" + libraryId + ".R2.fastq") + val samToFastq = SamToFastq(qscript, bamFile, libDir + sampleId + "-" + libId + ".R1.fastq", + libDir + sampleId + "-" + libId + ".R2.fastq") samToFastq.isIntermediate = true qscript.add(samToFastq) mapping.input_R1 = samToFastq.fastqR1 @@ -83,17 +83,17 @@ class GatkPipeline(val root: Configurable) extends QScript with MultiSampleQScri val header = inputSam.getFileHeader.getReadGroups for (readGroup <- inputSam.getFileHeader.getReadGroups) { if (readGroup.getSample != sampleId) logger.warn("Sample ID readgroup in bam file is not the same") - if (readGroup.getLibrary != libraryId) logger.warn("Library ID readgroup in bam file is not the same") - if (readGroup.getSample != sampleId || readGroup.getLibrary != libraryId) readGroupOke = false + if (readGroup.getLibrary != libId) logger.warn("Library ID readgroup in bam file is not the same") + if (readGroup.getSample != sampleId || readGroup.getLibrary != libId) readGroupOke = false } inputSam.close if (!readGroupOke) { if (config("correct_readgroups", default = false)) { logger.info("Correcting readgroups, file:" + bamFile) - val aorrg = AddOrReplaceReadGroups(qscript, bamFile, new File(libDir + sampleId + "-" + libraryId + ".bam")) - aorrg.RGID = sampleId + "-" + libraryId - aorrg.RGLB = libraryId + val aorrg = AddOrReplaceReadGroups(qscript, bamFile, new File(libDir + sampleId + "-" + libId + ".bam")) + aorrg.RGID = sampleId + "-" + libId + aorrg.RGLB = libId aorrg.RGSM = sampleId aorrg.isIntermediate = true qscript.add(aorrg) @@ -106,7 +106,7 @@ class GatkPipeline(val root: Configurable) extends QScript with MultiSampleQScri Some(bamFile) } } else { - logger.error("Sample: " + sampleId + ": No R1 found for run: " + libraryId) + logger.error("Sample: " + sampleId + ": No R1 found for run: " + libId) None } diff --git a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/core/MultiSampleQScript.scala b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/core/MultiSampleQScript.scala index 0217615d1..eb76502c8 100644 --- a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/core/MultiSampleQScript.scala +++ b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/core/MultiSampleQScript.scala @@ -40,26 +40,26 @@ trait MultiSampleQScript extends BiopetQScript { /** * Library class with basic functions build in - * @param libraryId + * @param libId */ - abstract class AbstractLibrary(val libraryId: String) { + abstract class AbstractLibrary(val libId: String) { /** Overrules config of qscript with default sample and default library */ - val config = new ConfigFunctions(defaultSample = sampleId, defaultLibrary = libraryId) + val config = new ConfigFunctions(defaultSample = sampleId, defaultLibrary = libId) /** Adds the library jobs */ final def addAndTrackJobs(): Unit = { currentSample = Some(sampleId) - currentLib = Some(libraryId) + currentLib = Some(libId) addJobs() currentLib = None currentSample = None } /** Creates a library file with given suffix */ - def createFile(suffix: String): File = new File(libDir, sampleId + "-" + libraryId + suffix) + def createFile(suffix: String): File = new File(libDir, sampleId + "-" + libId + suffix) /** Returns library directory */ - def libDir = sampleDir + "lib_" + libraryId + File.separator + def libDir = sampleDir + "lib_" + libId + File.separator /** Function that add library jobs */ protected def addJobs() @@ -95,7 +95,7 @@ trait MultiSampleQScript extends BiopetQScript { /** function add all libraries in one call */ protected final def addPerLibJobs(): Unit = { - for ((libraryId, library) <- libraries) { + for ((libId, library) <- libraries) { library.addAndTrackJobs() } } diff --git a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/MultisamplePipelineTemplate.scala b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/MultisamplePipelineTemplate.scala index 271c06085..7d46c4b8f 100644 --- a/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/MultisamplePipelineTemplate.scala +++ b/public/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/pipelines/MultisamplePipelineTemplate.scala @@ -26,7 +26,7 @@ class MultisamplePipelineTemplate(val root: Configurable) extends QScript with M class Sample(sampleId: String) extends AbstractSample(sampleId) { def makeLibrary(id: String) = new Library(id) - class Library(libraryId: String) extends AbstractLibrary(libraryId) { + class Library(libId: String) extends AbstractLibrary(libId) { protected def addJobs(): Unit = { // Library jobs } 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 b06c01f57..007636a63 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 @@ -44,14 +44,14 @@ class Carp(val root: Configurable) extends QScript with MultiSampleQScript { def makeSample(id: String) = new Sample(id) class Sample(sampleId: String) extends AbstractSample(sampleId) { def makeLibrary(id: String) = new Library(id) - class Library(libraryId: String) extends AbstractLibrary(libraryId) { + class Library(libId: String) extends AbstractLibrary(libId) { val mapping = new Mapping(qscript) def addJobs(): Unit = { if (config.contains("R1")) { mapping.input_R1 = config("R1") if (config.contains("R2")) mapping.input_R2 = config("R2") - mapping.libraryId = libraryId + mapping.libId = libId mapping.sampleId = sampleId mapping.outputDir = libDir @@ -59,7 +59,7 @@ class Carp(val root: Configurable) extends QScript with MultiSampleQScript { mapping.biopetScript addAll(mapping.functions) - } else logger.error("Sample: " + sampleId + ": No R1 found for library: " + libraryId) + } else logger.error("Sample: " + sampleId + ": No R1 found for library: " + libId) } } diff --git a/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/Flexiprep.scala b/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/Flexiprep.scala index c8cb6f6a6..b6bf15303 100644 --- a/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/Flexiprep.scala +++ b/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/Flexiprep.scala @@ -38,7 +38,7 @@ class Flexiprep(val root: Configurable) extends QScript with BiopetQScript { /** Skip Clip fastq files */ var skipClip: Boolean = config("skip_clip", default = false) - // TODO: hide sampleId and libraryId from the command line so they do not interfere with our config values + // TODO: hide sampleId and libId from the command line so they do not interfere with our config values /** Sample name */ @Argument(doc = "Sample ID", shortName = "sample", required = true) @@ -46,7 +46,7 @@ class Flexiprep(val root: Configurable) extends QScript with BiopetQScript { /** Library name */ @Argument(doc = "Library ID", shortName = "library", required = true) - var libraryId: String = _ + var libId: String = _ var paired: Boolean = (input_R2 != null) var R1_ext: String = _ @@ -65,7 +65,7 @@ class Flexiprep(val root: Configurable) extends QScript with BiopetQScript { if (input_R1 == null) throw new IllegalStateException("Missing R1 on flexiprep module") if (outputDir == null) throw new IllegalStateException("Missing Output directory on flexiprep module") if (sampleId == null) throw new IllegalStateException("Missing Sample name on flexiprep module") - if (libraryId == null) throw new IllegalStateException("Missing Library name on flexiprep module") + if (libId == null) throw new IllegalStateException("Missing Library name on flexiprep module") else if (!outputDir.endsWith("/")) outputDir += "/" paired = (input_R2 != null) @@ -83,7 +83,7 @@ class Flexiprep(val root: Configurable) extends QScript with BiopetQScript { R2_name = R2_name.substring(0, R2_name.lastIndexOf(R2_ext)) } - summary.out = outputDir + sampleId + "-" + libraryId + ".qc.summary.json" + summary.out = outputDir + sampleId + "-" + libId + ".qc.summary.json" } def biopetScript() { diff --git a/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/FlexiprepSummary.scala b/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/FlexiprepSummary.scala index 84acd1820..353fb1acb 100644 --- a/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/FlexiprepSummary.scala +++ b/public/flexiprep/src/main/scala/nl/lumc/sasc/biopet/pipelines/flexiprep/FlexiprepSummary.scala @@ -121,7 +121,7 @@ class FlexiprepSummary(val root: Configurable) extends InProcessFunction with Co md5Summary() val summary = ("samples" := ( flexiprep.sampleId := - ("libraries" := ( flexiprep.libraryId := ( + ("libraries" := ( flexiprep.libId := ( ("flexiprep" := ( ("clipping" := !flexiprep.skipClip) ->: ("trimming" := !flexiprep.skipTrim) ->: diff --git a/public/kopisu/src/main/scala/nl/lumc/sasc/biopet/pipelines/kopisu/Kopisu.scala b/public/kopisu/src/main/scala/nl/lumc/sasc/biopet/pipelines/kopisu/Kopisu.scala index 5360b4ac7..94b3bcbe6 100644 --- a/public/kopisu/src/main/scala/nl/lumc/sasc/biopet/pipelines/kopisu/Kopisu.scala +++ b/public/kopisu/src/main/scala/nl/lumc/sasc/biopet/pipelines/kopisu/Kopisu.scala @@ -36,7 +36,7 @@ class Kopisu(val root: Configurable) extends QScript with MultiSampleQScript { def makeSample(id: String) = new Sample(id) class Sample(sampleId: String) extends AbstractSample(sampleId) { def makeLibrary(id: String) = new Library(id) - class Library(libraryId: String) extends AbstractLibrary(libraryId) { + class Library(libId: String) extends AbstractLibrary(libId) { def addJobs(): Unit = { } diff --git a/public/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala b/public/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala index 7ed86a669..c7e94ded3 100644 --- a/public/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala +++ b/public/mapping/src/main/scala/nl/lumc/sasc/biopet/pipelines/mapping/Mapping.scala @@ -67,7 +67,7 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { protected var readgroupId: String = _ /** Readgroup Library */ - var libraryId: String = _ + var libId: String = _ /** Readgroup Platform */ protected var platform: String = config("platform", default = "illumina") @@ -100,9 +100,9 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { if (input_R1 == null) throw new IllegalStateException("Missing FastQ R1 on mapping module") paired = (input_R2 != null) - if (libraryId == null) libraryId = config("library_id") + if (libId == null) libId = config("library_id") if (sampleId == null) sampleId = config("sample_id") - if (readgroupId == null && sampleId != null && libraryId != null) readgroupId = sampleId + "-" + libraryId + if (readgroupId == null && sampleId != null && libId != null) readgroupId = sampleId + "-" + libId else if (readgroupId == null) readgroupId = config("readgroup_id") if (outputName == null) outputName = readgroupId @@ -127,7 +127,7 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { flexiprep.input_R1 = input_R1 if (paired) flexiprep.input_R2 = input_R2 flexiprep.sampleId = this.sampleId - flexiprep.libraryId = this.libraryId + flexiprep.libId = this.libId flexiprep.init flexiprep.runInitialJobs } @@ -279,7 +279,7 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { var RG: String = "ID:" + readgroupId + "," RG += "SM:" + sampleId + "," - RG += "LB:" + libraryId + "," + RG += "LB:" + libId + "," if (readgroupDescription != null) RG += "DS" + readgroupDescription + "," RG += "PU:" + platformUnit + "," if (predictedInsertsize.getOrElse(0) > 0) RG += "PI:" + predictedInsertsize.get + "," @@ -330,7 +330,7 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { addOrReplaceReadGroups.createIndex = true addOrReplaceReadGroups.RGID = readgroupId - addOrReplaceReadGroups.RGLB = libraryId + addOrReplaceReadGroups.RGLB = libId addOrReplaceReadGroups.RGPL = platform addOrReplaceReadGroups.RGPU = platformUnit addOrReplaceReadGroups.RGSM = sampleId @@ -344,7 +344,7 @@ class Mapping(val root: Configurable) extends QScript with BiopetQScript { def getReadGroup(): String = { var RG: String = "@RG\\t" + "ID:" + readgroupId + "\\t" - RG += "LB:" + libraryId + "\\t" + RG += "LB:" + libId + "\\t" RG += "PL:" + platform + "\\t" RG += "PU:" + platformUnit + "\\t" RG += "SM:" + sampleId + "\\t" diff --git a/public/sage/src/main/scala/nl/lumc/sasc/biopet/pipelines/sage/Sage.scala b/public/sage/src/main/scala/nl/lumc/sasc/biopet/pipelines/sage/Sage.scala index 03be2ded6..1bc0837a3 100644 --- a/public/sage/src/main/scala/nl/lumc/sasc/biopet/pipelines/sage/Sage.scala +++ b/public/sage/src/main/scala/nl/lumc/sasc/biopet/pipelines/sage/Sage.scala @@ -62,16 +62,16 @@ class Sage(val root: Configurable) extends QScript with MultiSampleQScript { def makeSample(id: String) = new Sample(id) class Sample(sampleId: String) extends AbstractSample(sampleId) { def makeLibrary(id: String) = new Library(id) - class Library(libraryId: String) extends AbstractLibrary(libraryId) { + class Library(libId: String) extends AbstractLibrary(libId) { val inputFastq: File = config("R1", required = true) val prefixFastq: File = createFile(".prefix.fastq") val flexiprep = new Flexiprep(qscript) flexiprep.sampleId = sampleId - flexiprep.libraryId = libraryId + flexiprep.libId = libId val mapping = new Mapping(qscript) - mapping.libraryId = libraryId + mapping.libId = libId mapping.sampleId = sampleId protected def addJobs(): Unit = { @@ -96,8 +96,8 @@ class Sage(val root: Configurable) extends QScript with MultiSampleQScript { qscript.addAll(mapping.functions) if (config("library_counts", default = false).asBoolean) { - addBedtoolsCounts(mapping.finalBamFile, sampleId + "-" + libraryId, libDir) - addTablibCounts(pf.outputFastq, sampleId + "-" + libraryId, libDir) + addBedtoolsCounts(mapping.finalBamFile, sampleId + "-" + libId, libDir) + addTablibCounts(pf.outputFastq, sampleId + "-" + libId, libDir) } } } diff --git a/public/yamsvp/src/main/scala/nl/lumc/sasc/biopet/pipelines/yamsvp/Yamsvp.scala b/public/yamsvp/src/main/scala/nl/lumc/sasc/biopet/pipelines/yamsvp/Yamsvp.scala index 8e3ef0af5..206c92a69 100644 --- a/public/yamsvp/src/main/scala/nl/lumc/sasc/biopet/pipelines/yamsvp/Yamsvp.scala +++ b/public/yamsvp/src/main/scala/nl/lumc/sasc/biopet/pipelines/yamsvp/Yamsvp.scala @@ -124,11 +124,11 @@ class Yamsvp(val root: Configurable) extends QScript with BiopetQScript { //with // Called for each run from a sample - def runSingleLibraryJobs(libraryId: String, sampleID: String): LibraryOutput = { + def runSingleLibraryJobs(libId: String, sampleID: String): LibraryOutput = { val libraryOutput = new LibraryOutput val alignmentDir: String = outputDir + sampleID + "/alignment/" - val runDir: String = alignmentDir + "run_" + libraryId + "/" + val runDir: String = alignmentDir + "run_" + libId + "/" if (config.contains("R1")) { val mapping = new Mapping(this) @@ -140,7 +140,7 @@ class Yamsvp(val root: Configurable) extends QScript with BiopetQScript { //with mapping.input_R1 = config("R1") mapping.input_R2 = config("R2") mapping.paired = (mapping.input_R2 != null) - mapping.RGLB = libraryId + mapping.RGLB = libId mapping.RGSM = sampleID mapping.RGPL = config("PL") mapping.RGPU = config("PU") @@ -154,7 +154,7 @@ class Yamsvp(val root: Configurable) extends QScript with BiopetQScript { //with // start sambamba dedup libraryOutput.mappedBamFile = mapping.outputFiles("finalBamFile") - } else this.logger.error("Sample: " + sampleID + ": No R1 found for library: " + libraryId) + } else this.logger.error("Sample: " + sampleID + ": No R1 found for library: " + libId) return libraryOutput // logger.debug(outputFiles) // return outputFiles -- GitLab