Skip to content
Snippets Groups Projects
Commit 0de6794f authored by wyleung's avatar wyleung
Browse files

Converting GATK VariantEval to Queue, wraps around to respect the

config, issue #82
parent 14790fa3
No related branches found
No related tags found
No related merge requests found
......@@ -25,6 +25,5 @@ Any value defined here will override the pom.xml file value but is only applicab
<com-junichi11-netbeans-changelf.use-global>false</com-junichi11-netbeans-changelf.use-global>
<com-junichi11-netbeans-changelf.show-dialog>true</com-junichi11-netbeans-changelf.show-dialog>
<org-netbeans-modules-javascript2-requirejs.enabled>true</org-netbeans-modules-javascript2-requirejs.enabled>
<netbeans.hint.jdkPlatform>JDK_1.8</netbeans.hint.jdkPlatform>
</properties>
</project-shared-configuration>
package nl.lumc.sasc.biopet.extensions.gatk
import java.io.File
import nl.lumc.sasc.biopet.core.config.Configurable
class VariantEval(val root: Configurable) extends org.broadinstitute.gatk.queue.extensions.gatk.VariantEval with GatkGeneral {
override def afterGraph {
super.afterGraph
}
}
object VariantEval {
def apply(root: Configurable, sample: File, compareWith: File,
output: File): VariantEval = {
val vareval = new VariantEval(root)
vareval.eval = Seq(sample)
vareval.comp = Seq(compareWith)
vareval.out = output
vareval.afterGraph
return vareval
}
def apply(root: Configurable, sample: File, compareWith: File,
output: File, ST: Seq[String], EV: Seq[String]): VariantEval = {
val vareval = new VariantEval(root)
vareval.eval = Seq(sample)
vareval.comp = Seq(compareWith)
vareval.out = output
vareval.noST = true
vareval.ST = ST
vareval.noEV = true
vareval.EV = EV
vareval.afterGraph
return vareval
}
}
\ No newline at end of file
......@@ -5,8 +5,9 @@ import java.io.File
import nl.lumc.sasc.biopet.core.config.Configurable
import nl.lumc.sasc.biopet.extensions.gatk.CombineVariants
import nl.lumc.sasc.biopet.extensions.gatk.SelectVariants
import nl.lumc.sasc.biopet.extensions.gatk.VariantEval
import org.broadinstitute.gatk.queue.QScript
import org.broadinstitute.gatk.queue.extensions.gatk.{ CommandLineGATK, VariantEval }
import org.broadinstitute.gatk.queue.extensions.gatk.{ CommandLineGATK }
import org.broadinstitute.gatk.utils.commandline.{ Input, Argument }
class GatkVcfSampleCompare(val root: Configurable) extends QScript with BiopetQScript {
......@@ -30,8 +31,6 @@ class GatkVcfSampleCompare(val root: Configurable) extends QScript with BiopetQS
trait gatkArguments extends CommandLineGATK {
this.reference_sequence = reference
this.memoryLimit = 2
this.jobResourceRequests :+= "h_vmem=4G"
}
def init() {
......@@ -70,14 +69,13 @@ class GatkVcfSampleCompare(val root: Configurable) extends QScript with BiopetQS
for ((sample, sampleVcf) <- sampleVcfs) {
val sampleDir = generalSampleDir + sample + File.separator
for ((compareSample, compareSampleVcf) <- sampleVcfs) {
val variantEval = new VariantEval with gatkArguments
variantEval.eval = Seq(sampleVcf)
variantEval.comp = Seq(compareSampleVcf)
variantEval.out = new File(sampleDir + sample + "-" + compareSample + ".eval.txt")
variantEval.noST = true
variantEval.ST = Seq("VariantType", "CompRod")
variantEval.noEV = true
variantEval.EV = Seq("CompOverlap")
val variantEval = VariantEval(this,
sampleVcf,
compareSampleVcf,
new File(sampleDir + sample + "-" + compareSample + ".eval.txt"),
Seq("VariantType", "CompRod"),
Seq("CompOverlap")
)
if (targetBed != null) variantEval.L = targetBed
add(variantEval)
sampleCompareMetrics.deps ::= variantEval.out
......
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