Question: How to configure genome in "Map with BWA" (existing tool) after adding it to my local section
0
gravatar for 18Kbeyond
7 months ago by
18Kbeyond20
New York
18Kbeyond20 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 COMMENTlink modified 7 months ago by Jennifer Hillman Jackson24k • written 7 months ago by 18Kbeyond20
0
gravatar for Jennifer Hillman Jackson
7 months ago by
United States
Jennifer Hillman Jackson24k 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 COMMENTlink written 7 months ago by Jennifer Hillman Jackson24k

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 REPLYlink modified 6 months ago • written 7 months ago by 18Kbeyond20

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 REPLYlink modified 6 months ago • written 6 months ago by 18Kbeyond20
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: 92 users visited in the last hour