annotate jammwrapper.py @ 0:d42f4d78c85e draft

Uploaded
author messersc
date Wed, 17 Dec 2014 10:40:23 -0500
parents
children 243f75d0ed6e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
1 # This is a wrapper for JAMM to use it in Galaxy
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
2 # It will be accompanied with a jamm.xml file, which specifies the interface that Galaxy is showing to the user
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
3 # as well as the JAMM software including JAMM.sh, the bash script that is actually called.
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
4
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
5 # This wrapper does the following things:
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
6 # map the files (from Galaxy history) to a directory
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
7 # pass the parameters from the GUI to JAMM.sh
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
8 # call JAMM.sh
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
9 # maybe map the resulting tabular files back to history items, the XML could take care of that if output is fixed
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
10
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
11 #import optparse
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
12 import argparse, os, shutil, subprocess, sys, tempfile
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
13 import shlex
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
14 # importing some of the modules used, especially for symlinking, dir manipulation and tool calling.
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
15 # since python2.7, argparse is preferred over optparse.
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
16
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
17 # for reference wrappers, see
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
18 # https://bitbucket.org/fubar/rossdev/src/4a91f99b5e1953270c9b35d2ca70c325a10fcfcb/rgdev/bwa_wrapper.py?at=default
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
19 # https://www.e-biogenouest.org/wiki/WrapperpythonforGalaxy:Asyntaxexample
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
20
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
21 def main():
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
22
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
23 #Read command line arguments
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
24 parser = argparse.ArgumentParser()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
25 parser.add_argument('-i', dest = 'input', nargs='+')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
26 parser.add_argument('-g', dest = 'gsize')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
27 parser.add_argument('-z', dest = 'peakfile')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
28
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
29 parser.add_argument('-m', dest = 'mode')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
30 parser.add_argument('-r', dest = 'resolution')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
31 parser.add_argument('-p', dest = 'processes')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
32 #parser.add_argument('-b', dest = 'binSize')
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
33
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
34 args = parser.parse_args()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
35
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
36 print "\n"
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
37 print "Wrapper debugging"
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
38 print "Files to be used:"
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
39 for j in args.input:
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
40 print j
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
41
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
42 print "output file:"
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
43 print args.peakfile
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
44
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
45 # depracted, can still be found in may of the example wrappers
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
46 #parser = optparse.OptionParser()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
47 #parser.add_option( '-i', dest='input', help='input bed files' )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
48 #parser.add_option( '-c', dest='csize', help='chr sizes' )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
49 #(options, args) = parser.parse_args()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
50
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
51 # create temp dir
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
52 tmp_dir = tempfile.mkdtemp()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
53 os.chdir(tmp_dir)
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
54 # symlink creation
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
55 for file in args.input:
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
56 filen = tmp_dir + "/" + os.path.basename(os.path.splitext(file)[0])+".bed"
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
57 os.symlink(file, filen)
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
58
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
59 # in case temp files should have random names
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
60 # ref_file = tempfile.NamedTemporaryFile( dir=tmp_dir )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
61 # ref_file_name = ref_file.name
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
62 # ref_file.close()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
63 # os.symlink( file, ref_file_name )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
64
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
65 command = ( "/home/cmesser/galaxy/tools/jamm/JAMM.sh -s %s -g %s -o results -m %s -r %s -p %s"
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
66 % ( tmp_dir, args.gsize, args.mode, args.resolution, args.processes ) )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
67 print command
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
68 # depending on how your programm is called, it may be necessary to use shlex.split on the command string before
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
69 # in this case, this was actually harmful. idk why
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
70 # command = shlex.split(command)
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
71 # Please read the docs for subprocess.Popen. You can't just pass as string as variables need to be extended
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
72 # Note that shell = True can be a security hazard.
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
73 # proc = subprocess.Popen( command, executable = '/bin/bash', shell=True )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
74 proc = subprocess.Popen( command, shell=True)
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
75 returncode = proc.wait()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
76
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
77 #mv files to a place where galaxy wrapper can find them
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
78 mvcommand = "mv %s/results/peaks/all.peaks.narrowPeak %s" % ( tmp_dir, args.peakfile )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
79 os.system(mvcommand)
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
80
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
81 # clean up temp dir
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
82 if os.path.exists( tmp_dir ):
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
83 shutil.rmtree( tmp_dir )
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
84
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
85
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
86 if __name__ == "__main__":
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
87 main()
d42f4d78c85e Uploaded
messersc
parents:
diff changeset
88