UnifiedGenotyperAllele.scala 700 Bytes
Newer Older
1
2
3
4
package nl.lumc.sasc.biopet.pipelines.gatk.variantcallers

import nl.lumc.sasc.biopet.pipelines.shiva.variantcallers.Variantcaller
import nl.lumc.sasc.biopet.utils.config.Configurable
Peter van 't Hof's avatar
Peter van 't Hof committed
5
import nl.lumc.sasc.biopet.extensions.gatk.broad
6
7
8
9
10
11
12

/** Allele mode for GenotyperAllele */
class UnifiedGenotyperAllele(val root: Configurable) extends Variantcaller {
  val name = "unifiedgenotyper_allele"
  protected def defaultPrio = 9

  def biopetScript() {
Peter van 't Hof's avatar
Peter van 't Hof committed
13
    val ug = broad.UnifiedGenotyper(this, inputBams.values.toList, outputFile)
14
15
16
17
18
    ug.alleles = config("input_alleles")
    ug.genotyping_mode = org.broadinstitute.gatk.tools.walkers.genotyper.GenotypingOutputMode.GENOTYPE_GIVEN_ALLELES
    add(ug)
  }
}