diff qualimap_bamqc.py @ 34:7e359370b5a2 draft

planemo upload for repository https://github.com/scottx611x/qualimap2 commit dc78b7c4b1780b316ca4aba2be247969ac1100ec-dirty
author scottx611x
date Thu, 26 Jul 2018 16:52:28 -0400
parents 0ea51e29a0ce
children a3e9c80894dd
line wrap: on
line diff
--- a/qualimap_bamqc.py	Thu Jul 26 16:47:12 2018 -0400
+++ b/qualimap_bamqc.py	Thu Jul 26 16:52:28 2018 -0400
@@ -9,12 +9,12 @@
 OUTPUT_DIR = "qualimap_results"
 
 
-def qualimap_bamqc(bam_filename, genomecov_file, jv_mem_size):
+def qualimap_bamqc(bam_filename, genomecov_file, out_dir, jv_mem_size):
     qualimap_command = [
         "qualimap", "bamqc",
         "-bam " + bam_filename,
         "-oc " + genomecov_file,
-        "-outdir " + OUTPUT_DIR,
+        "-outdir " + out_dir,
         "--java-mem-size=" + jv_mem_size
     ]
 
@@ -30,24 +30,26 @@
     )
     parser.add_argument('--input_file')
     parser.add_argument('--out_genome_file')
+    parser.add_argument('--out_dir')
     parser.add_argument('--java_mem_size')
 
     args = parser.parse_args()
 
-    print(os.listdir(os.getcwd()))
+    print(os.listdir(args.out_dir))
 
     qualimap_bamqc(
         args.input_file,
         args.out_genome_file,
+        args.out_dir,
         args.java_mem_size
     )
 
-    shutil.move(OUTPUT_DIR + "/genome_results.txt", "./genome_results.txt")
-    shutil.move(OUTPUT_DIR + "/qualimapReport.html", "./qualimapReport.html")
+    shutil.move(args.out_dir + "/genome_results.txt", "./genome_results.txt")
+    shutil.move(args.out_dir + "/qualimapReport.html", "./qualimapReport.html")
     shutil.make_archive(
         './raw_data_qualimapReport',
         'zip',
-        OUTPUT_DIR + '/raw_data_qualimapReport'
+        args.out_dir + '/raw_data_qualimapReport'
     )