Forum Login Issue:
Currently the "Log in with Google" button redirects you to a "Page not found." This is an issue that our forum vendors are working on fixing. In the meantime, while on the "Page not found" you can edit the URL to delete the second gatk, firecloud, or wdl (depending on what subforum you are acessing).
ex: https://gatkforums.broadinstitute.org/gatk/gatk/entry/...

cannot find reference for member access

Hi,

I'm trying to build workflow to scatter genotyping over regions, but I keep getting the following error. I'm not sure how to fix it. Neither my editor nor womtool can give me more info.

Womtool validate Error:

ERROR: Cannot find reference to 'jGT' for member access 'jGT.GenomicsDBImport' (line 24, col 37):

             InputGDB         = jGT.GenomicsDBImport.genomicsDB,
                                    ^

Workflow:

import "jointGenoType.wdl" as jGT

workflow EOG {

  File RefFasta
  File RegionsFile
  File InputSamplesFile
  String MergedVCFname
  File? dbSNP

  scatter (region in read_lines("${RegionsFile}")){
    call jGT.ParseBedRegion{
      input: BedRegion = region
    }

    call jGT.GenomicsDBImport{
      input: RefFasta         = RefFasta,
             InputSamplesFile = InputSamplesFile,
             Region           = jGT.ParseBedRegion.region
    }

    call jGT.GenotypeFromGenomicsDB{
      input: RefFasta         = RefFasta,
             InputGDB         = jGT.GenomicsDBImport.genomicsDB,
             dbSNP            = dbSNP
    }
  }

  call jGT.GatherVCF{
    input: RefFasta      = RefFasta,
           OutputVCFname = MergedVCFname,
           InputVCFs     = jGT.GenotypeFromGenomicsDB.OutputVCF
  }

  call decompose_normalize{
    input: RefFasta      = RefFasta,
           InputVCF      = jGT.GatherVCF.OutputVCF
  }

  call jGT.MakeSitesOnlyVCF{
    input: RefFasta = RefFasta,
           InputVCF = decompose_normalize.OutputVCF
  }
}


task decompose_normalize{
  File InputVCF
  File RefFasta
  String OutputVCFname = sub(InputVCF, "\\.*", "_decomposed.vcf.gz")

  command{
    vt decompose \
      -s ${InputVCF} \
    | vt normalize \
    -r ${RefFasta}  - \
    | bgzip -c ${OutputVCFname}
   }

   runtime {
    docker: "quay.io/biocontainers/vt:2015.11.10--2"
  }

  output {
    File OutputVCF = OutputVCFname
  }
}

Also, since this is my first workflow, all protips are very welcome.

Thanks
M

Sign In or Register to comment.