Hi,
I think I may have spotted the problem. It looks like there are two tabs between "hg18" and "/share/apps/genome/human/bowtie/hg18/ hg18.fa", and I think this might cause the loc parsing to behave badly.
Otherwise, your loc file looks fine to me (assuming you mean to have it in a directory called bowtie). And hg18 is the correct way to refer to Human Mar. 2006, so everything should be matching up.
Let me know if that works or not.
Regards, Kelly
On Nov 9, 2009, at 2:31 PM, Juan Perin wrote:
Unfortunately that's the first thing I tried. The initial errors seemed to suggest that from the web interface. However, I made the addition of the reference and the indexed file (.fai) in my .loc file. I even went back to ensure I started the alignment from the beginning, just to be sure that the alignment was associated with the right reference file. I definitely have the two files (.fa and .fai) in the specified location and they are both fully readable by all users. I've tried this with two different reference sets as well.
Is there anywhere else I might find error oupput? Its difficult to traverse through the running output from run.sh, as it contains all the SQL etc... and the output that does exist seems to be sparce. The python error that I sent, which Nate, pointed out would fail on the command line without the proper env. variable was just a bad lead. I haven't gotten that to work on the command line yet, however, to see if it may give me more info.
My only guess, at the moment, is that maybe the reference that my alignment is associated to (Human Mar. 2006 (hg18)) does not match up with the hg18 that I have listed in my .loc files? My loc file entry looks like this:
index hg18 /share/apps/genome/human/bowtie/hg18/hg18.fa index Human_male_hg18 /share/apps/genome/human/bwa/hg18_male/ human_b36_male.fa index ZebraFish_v8 /share/apps/genome/zebrafish/v8/ZebrafishFull.fa
The dropdown menu that tells Galaxy which build the alignment refers to is based on the original hard coded hg18 listed above, so I'm thinking either the way I've named it is wrong and Galaxy doesn't know that my hg18 is the same as the hg18 listed, or I'm doing something else totally wrong. Thanks again.
Juan
On Nov 9, 2009, at 1:18 PM, Kelly Vincent wrote:
Hi Juan,
I believe that the problem is actually not the error that is printing below. Do you have the Samtools indices in place (specifically for hg18) and the loc file sam_fa_indices.loc created in the tool-data directory? Specifically, from my experience, it appears that you have to have both the index (hg18.fa.fai) and the original fasta file (hg18.fa) in the directory you specify in the loc file.
Let us know if that doesn't help, or if you need more information about setting up the loc file and indices.
Regards, Kelly
On Nov 9, 2009, at 10:59 AM, juan perin wrote:
Hi again,
I'm now stuck trying to get the samtools app's to work. I've finally gotten the bwa and bowtie tools to run an alignment, but my attempts to convert the resulting sam file to bam, using sam_to_bam fails. The error I get in the browser is:
An error occurred running this job: No sequences are available for 'hg18', request them by reporting this error.
The error, (after searching through the output for the command), seemed to suggest that I was missing an 'egg': This was the output:
python /home/perin/galaxy-dist/tools/samtools/sam_to_bam.py -- input1=/home/perin/galaxy-dist/database/files/000/dataset_70.dat -- dbkey=hg18 --ref_file="None" --output1=/home/perin/galaxy-dist/ database/files/000/dataset_72.dat --index_dir=/home/perin/galaxy- dist/tool-data Traceback (most recent call last): File "/home/perin/galaxy-dist/tools/samtools/sam_to_bam.py", line 17, in ? from galaxy import eggs ImportError: No module named galaxy
I'm guessing I'm either missing a config parameter that tells samtools how to work, or that I need to scramble another egg. Not sure, however, and was hoping for some help. I can't seem to find anything on the web. Thanks in advance.
Juan Perin
galaxy-dev mailing list galaxy-dev@bx.psu.edu http://mail.bx.psu.edu/cgi-bin/mailman/listinfo/galaxy-dev