in reply to Re: Command not executing
in thread Command not executing

Thank you for the replies. I forgot to check the output in @out. Turns out the program rsem cannot find the trans.fa that I'm using as input (even though it's there). Oddly enough, when I run the code again using the konsole and not via Galaxy, it works fine.

Replies are listed 'Best First'.
Re^3: Command not executing
by BrowserUk (Patriarch) on Jun 22, 2012 at 04:44 UTC
    program rsem cannot find the trans.fa that I'm using as input (even though it's there). Oddly enough, when I run the code again using the konsole and not via Galaxy, it works fine.

    That sounds like a classic "permissions issue".

    One way to track such issues is to add a couple of lines to the top of the script (rsem.pl in this case), that display the permissions available for each component of the path(s). Eg.

    sub checkPathPermissions { my( $path, @bits ) = split '[/\\]', shift; $path = '/' . $path; while( @bits ) { print STDERR "$path: ", join ' : ', stat( $path ); $path .= '/' . shift @bits; } }

    Note: That untested code; to illustrate the suggestion only.

    If you add that to the top of rsem and call it on each of the input paths; then run it from your console and via galaxy, the difference between the two runs will probably go a long way to diagnosing the issue.


    With the rise and rise of 'Social' network sites: 'Computers are making people easier to use everyday'
    Examine what is said, not who speaks -- Silence betokens consent -- Love the truth but pardon error.
    "Science is about questioning the status quo. Questioning authority".
    In the absence of evidence, opinion is indistinguishable from prejudice.

    The start of some sanity?