Question: Need help with "Unified Genotyper" tool
0
gravatar for jason.becker
3.9 years ago by
United States
jason.becker0 wrote:

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.1k views
ADD COMMENTlink modified 3.9 years ago by Jennifer Hillman Jackson25k • written 3.9 years ago by jason.becker0
1
gravatar for Jennifer Hillman Jackson
3.9 years ago by
United States
Jennifer Hillman Jackson25k 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 COMMENTlink written 3.9 years ago by Jennifer Hillman Jackson25k

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

ADD REPLYlink written 3.9 years ago by jason.becker0

Very glad all sorted out! Jen

ADD REPLYlink written 3.9 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 90 users visited in the last hour