UnifiedGenotyperAllele.scala 827 Bytes
Newer Older
Peter van 't Hof's avatar
Peter van 't Hof committed
1
2
3
4
5
/**
 * Due to the license issue with GATK, this part of Biopet can only be used inside the
 * LUMC. Please refer to https://git.lumc.nl/biopet/biopet/wikis/home for instructions
 * on how to use this protected part of biopet or contact us at sasc@lumc.nl
 */
6
package nl.lumc.sasc.biopet.pipelines.shiva.variantcallers
7

Peter van 't Hof's avatar
Peter van 't Hof committed
8
import nl.lumc.sasc.biopet.extensions.gatk.broad
9
import nl.lumc.sasc.biopet.utils.config.Configurable
10
11
12
13
14
15
16

/** 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
17
    val ug = broad.UnifiedGenotyper(this, inputBams.values.toList, outputFile)
18
    ug.alleles = config("input_alleles")
Peter van 't Hof's avatar
Peter van 't Hof committed
19
    ug.genotyping_mode = Some("GENOTYPE_GIVEN_ALLELES")
20
21
22
    add(ug)
  }
}