Sequences Are Not Currently Available For The Specified Build

Question: Sequences Are Not Currently Available For The Specified Build

7.1 years ago by

Mark Silby • 10 wrote:

Hi all, We are just beginning to analyze RNA-seq data and hope to make use of the Galaxy tools. We are not using a commonly studied organism such as human. Rather, our data are from the bacterium Pseudomonas fluorescens Pf0-1. The genome for P. fluorescens Pf0-1 appears in the "database/build" menu, and we have selected this in the .sam file attributes. When we run the sam to bam conversion, we get the error "Sequences are not currently available for the specified build." Is this something that can be corrected easily? Many thanks for your help. Mark Silby

galaxy • 1.4k views

ADD COMMENT • link •

modified 6.5 years ago by Jennifer Hillman Jackson ♦ 25k • written 7.1 years ago by Mark Silby • 10

7.1 years ago by

Jennifer Hillman Jackson ♦ 25k

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hello Mark, Not all genomes in the database assignment list are pre-index for all tools. That is the case for this genome. We can consider add this in, but it would actually be quicker for you to simply load a fasta version of the genome into your history and use it as a custom reference genome from there. Most tools have some variation of the choice on the SAM->BAM tool form's option: "Choose the source for the reference list:" set as "History". Load the genome using FTP, as explained here: http://wiki.g2.bx.psu.edu/Learn/Upload%20via%20FTP Best wishes for your project, Jen Galaxy team -- Jennifer Jackson http://usegalaxy.org http://galaxyproject.org/wiki/Support

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

Hi Jen,

I encountered the same issue and the genome I'm using is mm9. Can this be added to Galaxy? What is a fasta version of the genome mm9? Isn't the reference genome is the one we used for aligning?

ADD REPLY • link written 3.6 years ago by minh.nguyen • 20

6.5 years ago by

Richard Linchangco • 20

Richard Linchangco • 20 wrote:

Dear Galaxy, I've been searching the lists for this type of issue and only found one solution thus far which is the use of a custom reference. It doesn't make sense in my situation because the reference I used was from Galaxy itself when I mapped my data. I'm now trying to use GATK to find SNPs but no matter what I've tried I can't get past this issue. I'm trying to use the Count Covariates and the Unified Genotyper but to no avail. The only issue appears to be that "Sequences are not currently available for the specified build." Any help would be much appreciated. Thanks Sincerely, Rich -- Richard Linchangco PSM in Bioinformatics College of Computing and Informatics University of North Carolina at Charlotte Tel: (630)440-7010 rlinchan@uncc.edu, rlinch2@gmail.com

ADD COMMENT • link written 6.5 years ago by Richard Linchangco • 20

Hi Rich, Additional genomes will be specially sorted, indexed, and added to the GATK tool suite as it moves out of beta status. Hg19 is short-listed for addition near term. We do take requests to have genome added to tools and consider these when ranking our prioritization lists. Which genome did you want to use? One small warning when using a custom reference genome with this particular tool set - be sure to visit the GATK web site links directly to understand the sorting criteria for genomes. It can be different than how Galaxy, UCSC, and many of the existing tools already sort or instruct users to sort genomes or data. In short, the genome must be sorted in the exact order that it was originally released, but even this can be slightly confusing, especially if working with a non-human genome as there are few examples. Still, the documentation can help and tools are easily tested (if the sorting is wrong, the tool will fail and let you know). If others have requests for GATK native genomes, they are also welcome to reply. In general, key model organisms would be ranked highest in priority. We also try to get the largest genomes loaded natively first (for purely practical reasons). Good question, thanks! Jen Galaxy team -- Jennifer Jackson http://galaxyproject.org

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

6.5 years ago by

Jennifer Hillman Jackson ♦ 25k

United States

Jennifer Hillman Jackson ♦ 25k wrote:

Hi Rich, You could do this or you could switch to this version of the human genome (is in all of Galaxy's tools, but will not be useable at UCSC): Human (Homo sapiens) (b37): hg_g1k_v37 Thanks, Jen Galaxy team -- Jennifer Jackson http://galaxyproject.org

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

Please log in to add an answer.

Similar posts • Search »