Mercurial > repos > peterjc > mira4_assembler
diff tools/mira4/mira4.py @ 21:4abe8d59a438 draft
Uploaded v0.0.4 preview 1; fix getting BAM without MAF
author | peterjc |
---|---|
date | Tue, 28 Oct 2014 08:29:59 -0400 |
parents | aeb3e35f8236 |
children |
line wrap: on
line diff
--- a/tools/mira4/mira4.py Tue Jun 10 10:11:58 2014 -0400 +++ b/tools/mira4/mira4.py Tue Oct 28 08:29:59 2014 -0400 @@ -58,6 +58,9 @@ parser.add_option("--log", dest="log", default="-", metavar="FILE", help="MIRA logging output filename") +parser.add_option("-v", "--version", dest="version", + default=False, action="store_true", + help="Show version and quit") options, args = parser.parse_args() manifest = options.manifest out_maf = options.maf @@ -84,7 +87,7 @@ mira_convert_ver = get_version(mira_convert) if not mira_convert_ver.strip().startswith("4.0"): stop_err("This wrapper is for MIRA V4.0, not:\n%s\n%s" % (mira_ver, mira_convert)) -if "-v" in sys.argv or "--version" in sys.argv: +if options.version: print "%s, MIRA wrapper version %s" % (mira_ver, WRAPPER_VER) if mira_ver != mira_convert_ver: print "WARNING: miraconvert %s" % mira_convert_ver @@ -192,7 +195,11 @@ #For mapping mode, probably most people would expect a BAM file #using the reference FASTA file... if out_bam and out_bam != "-": - msg = make_bam(mira_convert, out_maf, ref_fasta, out_bam, handle) + if out_maf and out_maf != "-": + msg = make_bam(mira_convert, out_maf, ref_fasta, out_bam, handle) + else: + #Not collecting the MAF file, use original location + msg = make_bam(mira_convert, old_maf, ref_fasta, out_bam, handle) if msg: stop_err(msg) @@ -216,7 +223,9 @@ assert os.path.isdir(temp) d = "%s_assembly" % name -assert not os.path.isdir(d), "Path %s already exists" % d +#This can fail on my development machine if stale folders exist +#under Galaxy's .../database/job_working_directory/ tree: +assert not os.path.isdir(d), "Path %r already exists:\n%s" % (d, os.path.abspath(d)) try: #Check path access os.mkdir(d)