diff zinbra.py @ 6:7d532e820976

Support comparison in zinbra. 3fba761
author oleg@oleg-desktop.times.labs.intellij.net
date Wed, 28 Oct 2015 20:08:20 +0300
parents 0eb50728861a
children
line wrap: on
line diff
--- a/zinbra.py	Wed Oct 28 19:52:30 2015 +0300
+++ b/zinbra.py	Wed Oct 28 20:08:20 2015 +0300
@@ -1,7 +1,6 @@
 #!/usr/bin/env python
 """
 Python wrapper for zinbra.xml
-Usage: zinbra.py "${genome}" "${bed}" "${bin}" "${fdr}"
 """
 
 import os
@@ -10,7 +9,14 @@
 
 argv = sys.argv[1:]
 print 'Arguments {0}'.format(argv)
-genome, bed, bin, fdr = argv
+
+# Check command
+if len(argv) == 5:
+    genome, bin, fdr, action, bed = argv
+    assert action == "analyze"
+else:
+    genome, bin, fdr, action, bed1, bed2 = argv
+    assert action == "compare"
 
 # Configure main jar path
 jar = os.environ.get("INTEGRATION_JAR")
@@ -21,10 +27,11 @@
 subprocess.check_call(cmd, cwd=None, shell=True)
 
 # See https://github.com/JetBrains-Research/zinbra for command line options
-# cla.argument_string_list() is configured at ZinbraApplications#rebuildArgumentStringsFromVisualOptions
-cmd = 'java -cp {0} org.jetbrains.bio.zinbra.ZinbraCLA ' \
-      'analyze -i {1} -bed result.bed -r reference.2bit -b {3} -fdr {4}'.format(jar,
-                                                                                bed,
-                                                                                bin, fdr)
+if action == "analyze":
+    cmd = 'java -cp {0} org.jetbrains.bio.zinbra.ZinbraCLA ' \
+          'analyze -i {1} -bed result.bed -r reference.2bit -b {3} -fdr {4}'.format(jar, bed, bin, fdr)
+else:
+    cmd = 'java -cp {0} org.jetbrains.bio.zinbra.ZinbraCLA ' \
+          'compare -1 {1} -2 {2} -bed result.bed -r reference.2bit -b {3} -fdr {4}'.format(jar, bed1, bed2, bin, fdr)
 print 'Launching zinbra: {0}'.format(cmd)
 subprocess.check_call(cmd, cwd=None, shell=True)