Mercurial > repos > scottx611x > qualimap2_bamqc
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' )