The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Register now for the upcoming GATK Best Practices workshop, Nov 7-8 at the Broad in Cambridge, MA. Open to all comers! More info and signup at

VCF from GenotypeGVCFs missing fields

dklevebringdklevebring Posts: 79Member

So I've just run GenotypeGVCFs and get some odd looking things. The FORMAT field is GT:AD:DP:GQ:PL and on most occasions, the values are 1/1:0,1,0:1:3:10,3,0, i.e. 5 fields. However, sometimes I get ./.:.:317, only 3 fields.

I use queue to run this, so my case class is this:

case class genotypeGVCF(gvcfs:List[File], outvcf:File) extends GenotypeGVCFs with ExternalCommonArgs with FourteenCoreJob with NineDayJob {
    @Input(doc="list of gvcfs to merge") val inGVCFs = gvcfs
    @Output(doc="SNP calls in VCF format") val outVCF = outvcf
    this.variant = inGVCFs
    this.out = outVCF
    this.R = refGenome
    this.dbsnp = dbSNP
    this.jobName = "genotypeGVCF " + outvcf.getName()
    this.analysisName = this.jobName
    this.isIntermediate = false
    this.nt = 14
    this.intervals = Seq(targetIntervals)
    this.memoryLimit = 112

Best Answer


Sign In or Register to comment.