Need help with "Unified Genotyper" tool

Question: Need help with "Unified Genotyper" tool

4.3 years ago by

United States

Every time I try to use the Unified Genotyper Tool under GATK, I receive the following error message: "Sequences are not currently available for the specified build." I have utilized the BAM files for other analysis and they work fine. In addition, I replaced the read groups using PICARD and the modified BAM file was still not accessible by the Unified Genotyper Tool.

Tool name: Unified Genotyper
Tool version: 0.0.6
Tool ID: toolshed.g2.bx.psu.edu/repos/devteam/unified_genotyper/gatk_unified_genotyper/0.0.6
ToolShed URL: https://toolshed.g2.bx.psu.edu/view/devteam/unified_genotyper

unified-genotyper • 1.3k views

ADD COMMENT • link •

modified 4.3 years ago by Jennifer Hillman Jackson ♦ 25k • written 4.3 years ago by jason.becker • 0

4.3 years ago by

Jennifer Hillman Jackson ♦ 25k

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

When using the public Main Galaxy instance at http://usegalaxy.org, this tool is natively indexed for the reference genome "hg_g1k_b37" when using the option "Choose the source for the reference list: Locally cashed". This is the 1000 genomes version of the GrCH37 Human build, sorted in a "GATK" specific manor.

If not using this specific build or if there is a mismatch, the problem you describe will result with this tool. Specifically, when using "Choose the source for the reference list:" as "locally cashed" with nothing selected for ""Using reference file" or if hg_g1k_b37 is a mismatch for the input. Or when using "Choose the source for the reference list:" as "history" without specifying a Custom reference genome (option "Using reference file: <fasta dataset in your history>").

More about Custom reference genomes:
http://wiki.galaxyproject.org/Support#Custom_reference_genome

Perhaps one of these is the issue? Please examine the form parameters you used (click on the "re-run" icon, the blue double arrow in the dataset box) to check, and adjust as needed.

These tools (all of GATK Beta - v1.4, or you can try GATK2), can be used on a cloud or production local Galaxy, where additional local reference genomes can be added (necessary for larger genomes - smaller Custom reference genomes tend to work best with this tool due to the resources it requires). Should you decide to go this route, here are some links to get started:
http://usegalaxy.org/cloud or http://getgalaxy.org (Galaxy)
http:///usegalaxy.org/cloud (Tools, these and many others)

Please let us know if you need more assistance, Jen, Galaxy team

ADD COMMENT • link written 4.3 years ago by Jennifer Hillman Jackson ♦ 25k

I realigned using hg_g1k_b37 and GATK worked great. Thanks for the help.

ADD REPLY • link written 4.3 years ago by jason.becker • 0

Very glad all sorted out! Jen

ADD REPLY • link written 4.3 years ago by Jennifer Hillman Jackson ♦ 25k

Similar posts • Search »