Mercurial > repos > peterjc > mira4_assembler
annotate tools/mira4_0/mira4.py @ 35:259891fce7fd draft
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
| author | peterjc |
|---|---|
| date | Wed, 09 Aug 2017 13:32:11 -0400 |
| parents | 0785a6537f3e |
| children | cee8f9005e43 |
| rev | line source |
|---|---|
| 25 | 1 #!/usr/bin/env python |
| 2 """A simple wrapper script to call MIRA and collect its output. | |
| 3 """ | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
4 |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
5 from __future__ import print_function |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
6 |
| 25 | 7 import os |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
8 import shutil |
| 25 | 9 import subprocess |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
10 import sys |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
11 import tempfile |
| 25 | 12 import time |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
13 |
| 25 | 14 from optparse import OptionParser |
| 15 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
16 # Do we need any PYTHONPATH magic? |
| 25 | 17 from mira4_make_bam import make_bam |
| 18 | |
|
35
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
19 WRAPPER_VER = "0.0.11" # Keep in sync with the XML file |
| 25 | 20 |
| 21 | |
| 22 def get_version(mira_binary): | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
23 """Run MIRA to find its version number.""" |
| 25 | 24 # At the commend line I would use: mira -v | head -n 1 |
| 25 # however there is some pipe error when doing that here. | |
| 26 cmd = [mira_binary, "-v"] | |
| 27 try: | |
|
34
0785a6537f3e
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit 6405ba93fcec7ea93452bf54d559c7507ee7a57c
peterjc
parents:
32
diff
changeset
|
28 child = subprocess.Popen(cmd, universal_newlines=True, |
| 25 | 29 stdout=subprocess.PIPE, |
| 30 stderr=subprocess.STDOUT) | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
31 except Exception as err: |
| 25 | 32 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (" ".join(cmd), err)) |
| 33 sys.exit(1) | |
| 34 ver, tmp = child.communicate() | |
| 35 del child | |
| 36 return ver.split("\n", 1)[0].strip() | |
| 37 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
38 |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
39 # Parse Command Line |
| 25 | 40 usage = """Galaxy MIRA4 wrapper script v%s - use as follows: |
| 41 | |
| 42 $ python mira4.py ... | |
| 43 | |
| 44 This will run the MIRA binary and collect its output files as directed. | |
| 45 """ % WRAPPER_VER | |
| 46 parser = OptionParser(usage=usage) | |
| 47 parser.add_option("-m", "--manifest", dest="manifest", | |
| 48 default=None, metavar="FILE", | |
| 49 help="MIRA manifest filename") | |
| 50 parser.add_option("--maf", dest="maf", | |
| 51 default="-", metavar="FILE", | |
| 52 help="MIRA MAF output filename") | |
| 53 parser.add_option("--bam", dest="bam", | |
| 54 default="-", metavar="FILE", | |
| 55 help="Unpadded BAM output filename") | |
| 56 parser.add_option("--fasta", dest="fasta", | |
| 57 default="-", metavar="FILE", | |
| 58 help="Unpadded FASTA output filename") | |
| 59 parser.add_option("--log", dest="log", | |
| 60 default="-", metavar="FILE", | |
| 61 help="MIRA logging output filename") | |
| 62 parser.add_option("-v", "--version", dest="version", | |
| 63 default=False, action="store_true", | |
| 64 help="Show version and quit") | |
| 65 options, args = parser.parse_args() | |
| 66 manifest = options.manifest | |
| 67 out_maf = options.maf | |
| 68 out_bam = options.bam | |
| 69 out_fasta = options.fasta | |
| 70 out_log = options.log | |
| 71 | |
|
35
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
72 if "MIRA4" in os.environ: |
| 25 | 73 mira_path = os.environ["MIRA4"] |
|
35
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
74 mira_binary = os.path.join(mira_path, "mira") |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
75 if not os.path.isfile(mira_binary): |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
76 sys.exit("Missing mira under $MIRA4, %r\nFolder contained: %s" |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
77 % (mira_binary, ", ".join(os.listdir(mira_path)))) |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
78 mira_convert = os.path.join(mira_path, "miraconvert") |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
79 if not os.path.isfile(mira_convert): |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
80 sys.exit("Missing miraconvert under $MIRA4, %r\nFolder contained: %s" |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
81 % (mira_convert, ", ".join(os.listdir(mira_path)))) |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
82 else: |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
83 sys.stderr.write("DEBUG: Since $MIRA4 is not set, assuming mira binaries are on $PATH.\n") |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
84 mira_path = None |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
85 mira_binary = "mira" |
|
259891fce7fd
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit e4c56df75150c82d3e9c4ac487c4209108e52412
peterjc
parents:
34
diff
changeset
|
86 mira_convert = "miraconvert" |
| 25 | 87 |
| 88 mira_ver = get_version(mira_binary) | |
| 89 if not mira_ver.strip().startswith("4.0"): | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
90 sys.exit("This wrapper is for MIRA V4.0, not:\n%s\n%s" % (mira_ver, mira_binary)) |
| 25 | 91 mira_convert_ver = get_version(mira_convert) |
| 92 if not mira_convert_ver.strip().startswith("4.0"): | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
93 sys.exit("This wrapper is for MIRA V4.0, not:\n%s\n%s" % (mira_ver, mira_convert)) |
| 25 | 94 if options.version: |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
95 print("%s, MIRA wrapper version %s" % (mira_ver, WRAPPER_VER)) |
| 25 | 96 if mira_ver != mira_convert_ver: |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
97 print("WARNING: miraconvert %s" % mira_convert_ver) |
| 25 | 98 sys.exit(0) |
| 99 | |
| 100 if not manifest: | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
101 sys.exit("Manifest is required") |
| 25 | 102 elif not os.path.isfile(manifest): |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
103 sys.exit("Missing input MIRA manifest file: %r" % manifest) |
| 25 | 104 |
| 105 | |
| 106 try: | |
| 107 threads = int(os.environ.get("GALAXY_SLOTS", "1")) | |
| 108 except ValueError: | |
| 109 threads = 1 | |
| 110 assert 1 <= threads, threads | |
| 111 | |
| 112 | |
| 113 def override_temp(manifest): | |
| 114 """Override ``-DI:trt=/tmp`` in manifest with environment variable. | |
| 115 | |
| 116 Currently MIRA 4 does not allow envronment variables like ``$TMP`` | |
| 117 inside the manifest, which is a problem if you need to override | |
| 118 the default at run time. | |
| 119 | |
| 120 The tool XML will ``/tmp`` and we replace that here with | |
| 121 ``tempfile.gettempdir()`` which will respect $TMPDIR, $TEMP, $TMP | |
| 122 as explained in the Python standard library documentation: | |
| 123 http://docs.python.org/2/library/tempfile.html#tempfile.tempdir | |
| 124 | |
| 125 By default MIRA 4 would write its temporary files within the output | |
| 126 folder, which is a problem if that is a network drive. | |
| 127 """ | |
| 128 handle = open(manifest, "r") | |
| 129 text = handle.read() | |
| 130 handle.close() | |
| 131 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
132 # At time of writing, this is at the end of a file, |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
133 # but could be followed by a space in future... |
| 25 | 134 text = text.replace("-DI:trt=/tmp", "-DI:trt=" + tempfile.gettempdir()) |
| 135 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
136 # Want to try to ensure this gets written to disk before MIRA attempts |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
137 # to open it - any networked file system may impose a delay... |
| 25 | 138 handle = open(manifest, "w") |
| 139 handle.write(text) | |
| 140 handle.flush() | |
| 141 os.fsync(handle.fileno()) | |
| 142 handle.close() | |
| 143 | |
| 144 | |
| 145 def log_manifest(manifest): | |
| 146 """Write the manifest file to stderr.""" | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
147 sys.stderr.write("\n%s\nManifest file\n%s\n" % ("=" * 60, "=" * 60)) |
| 25 | 148 with open(manifest) as h: |
| 149 for line in h: | |
| 150 sys.stderr.write(line) | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
151 sys.stderr.write("\n%s\nEnd of manifest\n%s\n" % ("=" * 60, "=" * 60)) |
| 25 | 152 |
| 153 | |
| 154 def collect_output(temp, name, handle): | |
| 155 """Moves files to the output filenames (global variables).""" | |
| 156 f = "%s/%s_assembly/%s_d_results" % (temp, name, name) | |
| 157 if not os.path.isdir(f): | |
| 158 log_manifest(manifest) | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
159 sys.exit("Missing output folder") |
| 25 | 160 if not os.listdir(f): |
| 161 log_manifest(manifest) | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
162 sys.exit("Empty output folder") |
| 25 | 163 missing = [] |
| 164 | |
| 165 old_maf = "%s/%s_out.maf" % (f, name) | |
| 166 if not os.path.isfile(old_maf): | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
167 # Triggered extractLargeContigs.sh? |
| 25 | 168 old_maf = "%s/%s_LargeContigs_out.maf" % (f, name) |
| 169 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
170 # De novo or single strain mapping, |
| 25 | 171 old_fasta = "%s/%s_out.unpadded.fasta" % (f, name) |
| 172 ref_fasta = "%s/%s_out.padded.fasta" % (f, name) | |
| 173 if not os.path.isfile(old_fasta): | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
174 # Mapping (StrainX versus reference) or de novo |
| 25 | 175 old_fasta = "%s/%s_out_StrainX.unpadded.fasta" % (f, name) |
| 176 ref_fasta = "%s/%s_out_StrainX.padded.fasta" % (f, name) | |
| 177 if not os.path.isfile(old_fasta): | |
| 178 old_fasta = "%s/%s_out_ReferenceStrain.unpadded.fasta" % (f, name) | |
| 179 ref_fasta = "%s/%s_out_ReferenceStrain.padded.fasta" % (f, name) | |
| 180 | |
| 181 missing = False | |
| 182 for old, new in [(old_maf, out_maf), | |
| 183 (old_fasta, out_fasta)]: | |
| 184 if not os.path.isfile(old): | |
| 185 missing = True | |
| 186 elif not new or new == "-": | |
| 187 handle.write("Ignoring %s\n" % old) | |
| 188 else: | |
| 189 handle.write("Capturing %s\n" % old) | |
| 190 shutil.move(old, new) | |
| 191 if missing: | |
| 192 log_manifest(manifest) | |
| 193 sys.stderr.write("Contents of %r:\n" % f) | |
| 194 for filename in sorted(os.listdir(f)): | |
| 195 sys.stderr.write("%s\n" % filename) | |
| 196 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
197 # For mapping mode, probably most people would expect a BAM file |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
198 # using the reference FASTA file... |
| 25 | 199 if out_bam and out_bam != "-": |
| 200 if out_maf and out_maf != "-": | |
| 201 msg = make_bam(mira_convert, out_maf, ref_fasta, out_bam, handle) | |
| 202 else: | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
203 # Not collecting the MAF file, use original location |
| 25 | 204 msg = make_bam(mira_convert, old_maf, ref_fasta, out_bam, handle) |
| 205 if msg: | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
206 sys.exit(msg) |
| 25 | 207 |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
208 |
| 25 | 209 def clean_up(temp, name): |
| 210 folder = "%s/%s_assembly" % (temp, name) | |
| 211 if os.path.isdir(folder): | |
| 212 shutil.rmtree(folder) | |
| 213 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
214 |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
215 # TODO - Run MIRA in /tmp or a configurable directory? |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
216 # Currently Galaxy puts us somewhere safe like: |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
217 # /opt/galaxy-dist/database/job_working_directory/846/ |
| 25 | 218 temp = "." |
| 219 | |
| 220 name = "MIRA" | |
| 221 | |
| 222 override_temp(manifest) | |
| 223 | |
| 224 start_time = time.time() | |
| 225 cmd_list = [mira_binary, "-t", str(threads), manifest] | |
| 226 cmd = " ".join(cmd_list) | |
| 227 | |
| 228 assert os.path.isdir(temp) | |
| 229 d = "%s_assembly" % name | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
230 # This can fail on my development machine if stale folders exist |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
231 # under Galaxy's .../database/job_working_directory/ tree: |
| 25 | 232 assert not os.path.isdir(d), "Path %r already exists:\n%s" % (d, os.path.abspath(d)) |
| 233 try: | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
234 # Check path access |
| 25 | 235 os.mkdir(d) |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
236 except Exception as err: |
| 25 | 237 log_manifest(manifest) |
| 238 sys.stderr.write("Error making directory %s\n%s" % (d, err)) | |
| 239 sys.exit(1) | |
| 240 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
241 # print(os.path.abspath(".")) |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
242 # print(cmd) |
| 25 | 243 |
| 244 if out_log and out_log != "-": | |
| 245 handle = open(out_log, "w") | |
| 246 else: | |
| 247 handle = open(os.devnull, "w") | |
| 248 handle.write("======================== MIRA manifest (instructions) ========================\n") | |
| 249 m = open(manifest, "rU") | |
| 250 for line in m: | |
| 251 handle.write(line) | |
| 252 m.close() | |
| 253 del m | |
| 254 handle.write("\n") | |
| 255 handle.write("============================ Starting MIRA now ===============================\n") | |
| 256 handle.flush() | |
| 257 try: | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
258 # Run MIRA |
| 25 | 259 child = subprocess.Popen(cmd_list, |
| 260 stdout=handle, | |
| 261 stderr=subprocess.STDOUT) | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
262 except Exception as err: |
| 25 | 263 log_manifest(manifest) |
| 264 sys.stderr.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err)) | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
265 # TODO - call clean up? |
| 25 | 266 handle.write("Error invoking command:\n%s\n\n%s\n" % (cmd, err)) |
| 267 handle.close() | |
| 268 sys.exit(1) | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
269 # Use .communicate as can get deadlocks with .wait(), |
| 25 | 270 stdout, stderr = child.communicate() |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
271 assert not stdout and not stderr # Should be empty as sent to handle |
| 25 | 272 run_time = time.time() - start_time |
| 273 return_code = child.returncode | |
| 274 handle.write("\n") | |
| 275 handle.write("============================ MIRA has finished ===============================\n") | |
| 276 handle.write("MIRA took %0.2f hours\n" % (run_time / 3600.0)) | |
| 277 if return_code: | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
278 print("MIRA took %0.2f hours" % (run_time / 3600.0)) |
| 25 | 279 handle.write("Return error code %i from command:\n" % return_code) |
| 280 handle.write(cmd + "\n") | |
| 281 handle.close() | |
| 282 clean_up(temp, name) | |
| 283 log_manifest(manifest) | |
|
31
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
284 sys.stderr.write("Return error code %i from command:\n" % return_code) |
|
fd95aaef8818
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit bc3d484c5cd68ddcf456db2fff489d584aa2034c
peterjc
parents:
25
diff
changeset
|
285 sys.stderr.write(cmd + "\n") |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
286 sys.exit(return_code) |
| 25 | 287 handle.flush() |
| 288 | |
| 289 if os.path.isfile("MIRA_assembly/MIRA_d_results/ec.log"): | |
| 290 handle.write("\n") | |
| 291 handle.write("====================== Extract Large Contigs failed ==========================\n") | |
| 292 e = open("MIRA_assembly/MIRA_d_results/ec.log", "rU") | |
| 293 for line in e: | |
| 294 handle.write(line) | |
| 295 e.close() | |
| 296 handle.write("============================ (end of ec.log) =================================\n") | |
| 297 handle.flush() | |
| 298 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
299 # print("Collecting output...") |
| 25 | 300 start_time = time.time() |
| 301 collect_output(temp, name, handle) | |
| 302 collect_time = time.time() - start_time | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
303 handle.write("MIRA took %0.2f hours; collecting output %0.2f minutes\n" |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
304 % (run_time / 3600.0, collect_time / 60.0)) |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
305 print("MIRA took %0.2f hours; collecting output %0.2f minutes\n" |
|
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
306 % (run_time / 3600.0, collect_time / 60.0)) |
| 25 | 307 |
| 308 if os.path.isfile("MIRA_assembly/MIRA_d_results/ec.log"): | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
309 # Treat as an error, but doing this AFTER collect_output |
| 25 | 310 sys.stderr.write("Extract Large Contigs failed\n") |
| 311 handle.write("Extract Large Contigs failed\n") | |
| 312 handle.close() | |
| 313 sys.exit(1) | |
| 314 | |
|
32
56b421d59805
planemo upload for repository https://github.com/peterjc/galaxy_mira/tree/master/tools/mira4_0 commit fd979d17340cde155de176604744831d9597c6b6
peterjc
parents:
31
diff
changeset
|
315 # print "Cleaning up..." |
| 25 | 316 clean_up(temp, name) |
| 317 | |
| 318 handle.write("\nDone\n") | |
| 319 handle.close() | |
| 320 print("Done") |