The current GATK version is 3.7-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!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

#### ☞ Got a problem?

1. Search using the upper-right search box, e.g. using the error message.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

#### ☞ Formatting tip!

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks (  ) each to make a code block as demonstrated here.

**GATK4-BETA.2** is here. That's TWO, as in the second beta release. Be sure to read about the known issues before test driving. See Article#9881 to start and https://github.com/broadinstitute/gatk/blob/master/README.md for details.

# GATK to BEAGLE interface

University of PennsylvaniaMember

Greetings!

I'm trying to phase GATK genotype and to impute some SNP calls. Before I could do that, I must convert GATK results to an acceptable BEAGLE input format. What's the difference between VariantsToBeagleUnphased and ProduceBeagleInput? I know the latter outputs a file with genotype likelihoods. Incidentally, using that file didn't work in BEAGLE and produced the following log and error files. Can anyone give any pointers? Thanks in advance!

[stechen@node24 ~]$more beagle_run_410.o720239 Beagle version 3.3.2 (31 Oct 2011) Enter "java -jar beagle.jar" for summary of command line arguments. Start time: 11:59 AM EDT on 07 Aug 2013 Command line: java -Xmx7281m -jar beagle.jar like=beagle_input_410_impute phased=~/stechen/phase_ref/ALL.chr1.phase1_release_v3.20101123.filt.bgl markers=~stechen/phase_ref/ALL.chr1.phase1_release_v3.20101123.filt.markers missing=? out=beagle_output_410_chr1 [stechen@node24 ~]$ more beagle_run_410.e720239
bash: module: line 1: syntax error: unexpected end of file
bash: error importing function definition for module'
Picked up _JAVA_OPTIONS: -Djava.io.tmpdir=~stechen/tmp -Xms256m -Xmx8G
at phaser.y.a(Unknown Source)
at phaser.y.(Unknown Source)
at phaser.H.a(Unknown Source)
at phaser.v.(Unknown Source)
at phaser.PhaseMain.(Unknown Source)
at phaser.PhaseMain.main(Unknown Source)

Tagged:

Hi there,

For usage of the Beagle-related tools, we have a presentation that explains it in detail. Have a look at our Guide section, category "Events", under Best Practices Workshop. You'll find several presentations including one on genotype refinement. Let me know if you have trouble finding it. I'd give you a direct link but I'm doing this on a smartphone (which is acting not so smart right now).

For help with running Beagle software itself, please ask the Beagle maintainers for support.

• University of PennsylvaniaMember

Thanks for the info, Geraldine! May I ask how do I reach the Beagle maintainers?

• University of PennsylvaniaMember
edited August 2013

Also, I have indeed seen the workshop slides. Still, I'm wondering what's the difference between VariantsToBeagleUnphased and ProduceBeagleInput?

Also, in the GATK page, the command is given as "java -Xmx4000m -jar path_to_beagle/beagle.jar like=path_to_beagle_output/beagle_output out=myrun" What if I wanted to use a phased file as the template? It also seems to not work without a marker argument, is that correct?

Thanks again!