Question: How to configure genome in "Map with BWA" (existing tool) after adding it to my local section
0
gravatar for 18Kbeyond
14 months ago by
18Kbeyond • 20
New York
18Kbeyond • 20 wrote:

Hi. I have added the lines below to the section I made in config/tool_conf.xml file and "Map with BWA" is showing up on galaxy interface when I execute run.sh. However, I can't pick a genome in the "Using reference genome" box.

What modification should I make in shed_tools/toolshed.g2.bx.psu.edu/repos/devteam/bwa/051eba708f43/bwa/bwa.xml (or is there something else I need to do ) so as to add reference genome path? Thank you.

<tool file="../shed_tools/toolshed.g2.bx.psu.edu/repos/devteam/bwa/051eba708f43/bwa/bwa.xml" guid="toolshed.g2.bx.psu.edu/repos/devteam/bwa/bwa/0.7.15.2">
    <tool_shed>toolshed.g2.bx.psu.edu</tool_shed>
            <repository_name>bwa</repository_name>
            <repository_owner>devteam</repository_owner>
            <installed_changeset_revision>051eba708f43</installed_changeset_revision>
            <id>toolshed.g2.bx.psu.edu/repos/devteam/bwa/bwa/0.7.15.2</id>
            <version>0.7.15.2</version>
</tool>
ADD COMMENT • link • modified 14 months ago by Jennifer Hillman Jackson ♦ 25k • written 14 months ago by 18Kbeyond • 20
0
gravatar for Jennifer Hillman Jackson
14 months ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

Did you load any genomes and create indexes yet?

If not, install Data Managers from the Tool Shed and run them from the Admin area (local data). The usual order is:

  • Fetch the genome and create a dbkey
  • Samtools index
  • Picard index
  • TwoBit index
  • Followed by other tool indexes, including BWA

BWA for Illumina has been deprecated - upgrading to BWA/BWA-MEM for both the mapping tool and the Data Manager indexes is strongly recommended. Indexes for any tool are usually tool and tool version specific. Some exceptions are those that make use of the twoBit index (a few different tools use these same data).

Thanks, Jen, Galaxy team

ADD COMMENT • link written 14 months ago by Jennifer Hillman Jackson ♦ 25k

Thank you. Instead of installing data manager from the tool shed, I made my own genome and index and make changes in tool-data/bwa_mem_index.loc

hg19bwa hg19 hg19 Canonical /home/galaxy/galaxy_data/genomes/hg19/bwa/hg19.fa

mm10bwa mm10 mm10 Canonical /home/galaxy/galaxy_data/genomes/mm10/bwa/mm10.fa

ADD REPLY • link modified 14 months ago • written 14 months ago by 18Kbeyond • 20

However, when I execute the job, it generates the following error:

Traceback (most recent call last):

File "/data/galaxy/lib/galaxy/jobs/runners/local.py", line 130, in queue_job

job_wrapper.finish( stdout, stderr, exit_code )

File "/data/galaxy/lib/galaxy/jobs/__init__.py", line 1266, in finish

dataset.datatype.set_meta( dataset, overwrite=False )

File "/data/galaxy/lib/galaxy/datatypes/binary.py", line 394, in set_meta

exit_code = subprocess.call( args=command, stderr=open( stderr_name, 'wb' ) )

File "/usr/local/lib/python2.7/subprocess.py", line 522, in call

return Popen(popenargs, *kwargs).wait()

File "/usr/local/lib/python2.7/subprocess.py", line 710, in __init__ errread, errwrite)

File "/usr/local/lib/python2.7/subprocess.py", line 1335, in _execute_child raise child_exception

OSError: [Errno 2] No such file or directory

ADD REPLY • link modified 14 months ago • written 14 months ago by 18Kbeyond • 20
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: 172 users visited in the last hour