Skip to content
Snippets Groups Projects
Commit 48b7822f authored by Peter van 't Hof's avatar Peter van 't Hof
Browse files

Added wigle and tdf files to summary

parent 18184194
No related branches found
No related tags found
No related merge requests found
...@@ -38,6 +38,10 @@ class Bam2Wig(val root: Configurable) extends QScript with BiopetQScript { ...@@ -38,6 +38,10 @@ class Bam2Wig(val root: Configurable) extends QScript with BiopetQScript {
inputFiles :+= new InputFile(bamFile) inputFiles :+= new InputFile(bamFile)
} }
def outputWigleFile = new File(outputDir, bamFile.getName + ".wig")
def outputTdfFile = new File(outputDir, bamFile.getName + ".tdf")
def outputBwFile = new File(outputDir, bamFile.getName + ".bw")
def biopetScript(): Unit = { def biopetScript(): Unit = {
val bs = new BamToChromSizes(this) val bs = new BamToChromSizes(this)
bs.bamFile = bamFile bs.bamFile = bamFile
...@@ -48,14 +52,14 @@ class Bam2Wig(val root: Configurable) extends QScript with BiopetQScript { ...@@ -48,14 +52,14 @@ class Bam2Wig(val root: Configurable) extends QScript with BiopetQScript {
val igvCount = new IGVToolsCount(this) val igvCount = new IGVToolsCount(this)
igvCount.input = bamFile igvCount.input = bamFile
igvCount.genomeChromSizes = bs.chromSizesFile igvCount.genomeChromSizes = bs.chromSizesFile
igvCount.wig = Some(new File(outputDir, bamFile.getName + ".wig")) igvCount.wig = Some(outputWigleFile)
igvCount.tdf = Some(new File(outputDir, bamFile.getName + ".tdf")) igvCount.tdf = Some(outputTdfFile)
add(igvCount) add(igvCount)
val wigToBigWig = new WigToBigWig(this) val wigToBigWig = new WigToBigWig(this)
wigToBigWig.inputWigFile = igvCount.wig.get wigToBigWig.inputWigFile = igvCount.wig.get
wigToBigWig.inputChromSizesFile = bs.chromSizesFile wigToBigWig.inputChromSizesFile = bs.chromSizesFile
wigToBigWig.outputBigWig = new File(outputDir, bamFile.getName + ".bw") wigToBigWig.outputBigWig = outputBwFile
add(wigToBigWig) add(wigToBigWig)
} }
} }
......
...@@ -109,7 +109,14 @@ class Mapping(val root: Configurable) extends QScript with SummaryQScript with S ...@@ -109,7 +109,14 @@ class Mapping(val root: Configurable) extends QScript with SummaryQScript with S
/** File to add to the summary */ /** File to add to the summary */
def summaryFiles: Map[String, File] = Map("output_bam" -> finalBamFile, "input_R1" -> inputR1, def summaryFiles: Map[String, File] = Map("output_bam" -> finalBamFile, "input_R1" -> inputR1,
"reference" -> referenceFasta()) ++ "reference" -> referenceFasta()) ++
(if (inputR2.isDefined) Map("input_R2" -> inputR2.get) else Map()) (if (inputR2.isDefined) Map("input_R2" -> inputR2.get) else Map()) ++
(bam2wig match {
case Some(b) => Map(
"output_wigle" -> b.outputWigleFile,
"output_tdf" -> b.outputTdfFile,
"output_bigwig" -> b.outputBwFile)
case _ => Map()
})
/** Settings to add to summary */ /** Settings to add to summary */
def summarySettings = Map( def summarySettings = Map(
...@@ -270,12 +277,15 @@ class Mapping(val root: Configurable) extends QScript with SummaryQScript with S ...@@ -270,12 +277,15 @@ class Mapping(val root: Configurable) extends QScript with SummaryQScript with S
add(gears) add(gears)
} }
if (config("generate_wig", default = false).asBoolean) bam2wig.foreach(add(_))
addAll(Bam2Wig(this, finalBamFile).functions)
addSummaryJobs() addSummaryJobs()
} }
protected lazy val bam2wig = if (config("generate_wig", default = false)) {
Some(Bam2Wig(this, finalBamFile))
} else None
/** Add bwa aln jobs */ /** Add bwa aln jobs */
def addBwaAln(R1: File, R2: Option[File], output: File): File = { def addBwaAln(R1: File, R2: Option[File], output: File): File = {
val bwaAlnR1 = new BwaAln(this) val bwaAlnR1 = new BwaAln(this)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment