Question: import of refrence genome
0
gravatar for k_s_87
7 weeks ago by
k_s_870
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 • 107 views
ADD COMMENTlink modified 6 weeks ago • written 7 weeks ago by k_s_870
0
gravatar for Jennifer Hillman Jackson
7 weeks ago by
United States
Jennifer Hillman Jackson23k wrote:

Hello,

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: https://galaxyproject.org/admin/#data-resources

Thanks! Jen, Galaxy team

ADD COMMENTlink written 7 weeks ago by Jennifer Hillman Jackson23k
0
gravatar for k_s_87
6 weeks ago by
k_s_870
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 6 weeks 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: https://galaxyproject.org/learn/custom-genomes/

ADD REPLYlink modified 6 weeks ago • written 6 weeks ago by Jennifer Hillman Jackson23k

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 http://usegalaxy.org in a ticket which might help. You can also see the open minor issues that have come up (linked tickets). https://github.com/galaxyproject/galaxy/issues/1470

ADD REPLYlink modified 6 weeks ago • written 6 weeks ago by Jennifer Hillman Jackson23k
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: 104 users visited in the last hour