Question: import of refrence genome
gravatar for k_s_87
10 months ago by
k_s_870 wrote:

How can I add reference human genome to local galaxy instance? I've got fasta file and bwa index, but how can I upload it to the data mannager?

thanks kate

rna-seq • 289 views
ADD COMMENTlink modified 9 months ago • written 10 months ago by k_s_870
gravatar for Jennifer Hillman Jackson
9 months ago by
United States
Jennifer Hillman Jackson25k wrote:


Even if you already have the human genome, you might want to let the Data manager fetch the genome from the source (if it is a source listed on the DM). Or you can put your fasta version into the same history that you run DMs from, and input the fasta on the DM form as the "source from the history".

Rerunning the DM for BWA would be an easier way of getting the index into Galaxy (instead of manually adding in an existing index). But this is your decision.

Both choices are explained here, plus there is a DM tutorial:

Thanks! Jen, Galaxy team

ADD COMMENTlink written 9 months ago by Jennifer Hillman Jackson25k
gravatar for k_s_87
9 months ago by
k_s_870 wrote:

Hi, Thanks for your answer. I did as you mentioned and there is both reference genome (.fa) and index in data library. However, when I want to use FreeBayes tool, there is no option to use my local reference genome. How can I fix it?

Thanks, Kate

ADD COMMENTlink written 9 months ago by k_s_870

Any genome added to fasta_indexes.loc will be available to Freebayes. The tool does not use special tool-specific indexes.

You can also select a genome/transcriptome/exome (really any fasta target) "from the history" with this tool to create an index during runtime. As is with most Galaxy tools. How to:

ADD REPLYlink modified 9 months ago • written 9 months ago by Jennifer Hillman Jackson25k

Hi - Dan from our team reminded me that running Samtools indexes on an already fetched genome also populates the table that Freebayes uses for indexes. So try doing that.

I tend to load all genomes with these DMs, in this order:

1 - Fetch 
2 - Samtools 
3 - Picard
4 - 2bit
5 - Other tool-specific indexes (BWA, etc), any order

This helps to make sure most tools/functions get baseline data index files/tables installed and usable. Those that 1) come default in Galaxy or 2) are imported from the tool shed but do not have a specific data manager or index instructions the author added to the README. You can always look at the tip files for a repo for that info to check to see where it is getting index data by table name. This one is a bit confusing but most others match up well.

I keep track making indexes for in a ticket which might help. You can also see the open minor issues that have come up (linked tickets).

ADD REPLYlink modified 9 months ago • written 9 months ago by Jennifer Hillman Jackson25k
Please log in to add an answer.


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