comparison tools/mira3/mira.py @ 28:e28d6015fc56 draft default tip

planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira3/ commit 89578746a1c5b29c84a173d8b2709f086f69a7b6
author peterjc
date Mon, 03 Jun 2019 13:28:13 -0400
parents 85b86f31327c
children
comparison
equal deleted inserted replaced
27:85b86f31327c 28:e28d6015fc56
16 """Run MIRA to find its version number.""" 16 """Run MIRA to find its version number."""
17 # At the commend line I would use: mira -v | head -n 1 17 # At the commend line I would use: mira -v | head -n 1
18 # however there is some pipe error when doing that here. 18 # however there is some pipe error when doing that here.
19 cmd = ["mira", "-v"] 19 cmd = ["mira", "-v"]
20 try: 20 try:
21 child = subprocess.Popen(cmd, universal_newlines=True, 21 child = subprocess.Popen(
22 stdout=subprocess.PIPE, 22 cmd,
23 stderr=subprocess.STDOUT) 23 universal_newlines=True,
24 stdout=subprocess.PIPE,
25 stderr=subprocess.STDOUT,
26 )
24 except Exception as err: 27 except Exception as err:
25 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (" ".join(cmd), err)) 28 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (" ".join(cmd), err))
26 sys.exit(1) 29 sys.exit(1)
27 ver, tmp = child.communicate() 30 ver, tmp = child.communicate()
28 del child 31 del child
43 if not os.path.isdir(f): 46 if not os.path.isdir(f):
44 sys.exit("Missing output folder") 47 sys.exit("Missing output folder")
45 if not os.listdir(f): 48 if not os.listdir(f):
46 sys.exit("Empty output folder") 49 sys.exit("Empty output folder")
47 missing = [] 50 missing = []
48 for old, new in [("%s/%s_out.unpadded.fasta" % (f, name), out_fasta), 51 for old, new in [
49 ("%s/%s_out.unpadded.fasta.qual" % (f, name), out_qual), 52 ("%s/%s_out.unpadded.fasta" % (f, name), out_fasta),
50 ("%s/%s_out.wig" % (f, name), out_wig), 53 ("%s/%s_out.unpadded.fasta.qual" % (f, name), out_qual),
51 ("%s/%s_out.caf" % (f, name), out_caf), 54 ("%s/%s_out.wig" % (f, name), out_wig),
52 ("%s/%s_out.ace" % (f, name), out_ace)]: 55 ("%s/%s_out.caf" % (f, name), out_caf),
56 ("%s/%s_out.ace" % (f, name), out_ace),
57 ]:
53 if not os.path.isfile(old): 58 if not os.path.isfile(old):
54 missing.append(os.path.splitext(old)[-1]) 59 missing.append(os.path.splitext(old)[-1])
55 else: 60 else:
56 shutil.move(old, new) 61 shutil.move(old, new)
57 if missing: 62 if missing:
88 # print(cmd) 93 # print(cmd)
89 94
90 handle = open(out_log, "w") 95 handle = open(out_log, "w")
91 try: 96 try:
92 # Run MIRA 97 # Run MIRA
93 child = subprocess.Popen(cmd_list, universal_newlines=True, 98 child = subprocess.Popen(
94 stdout=handle, 99 cmd_list, universal_newlines=True, stdout=handle, stderr=subprocess.STDOUT
95 stderr=subprocess.STDOUT) 100 )
96 except Exception as err: 101 except Exception as err:
97 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err)) 102 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err))
98 # TODO - call clean up? 103 # TODO - call clean up?
99 handle.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err)) 104 handle.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err))
100 handle.close() 105 handle.close()