Mercurial > repos > scottx611x > qualimap2_bamqc
comparison qualimap_bamqc.py @ 41:e856fda6f375 draft
planemo upload for repository https://github.com/scottx611x/qualimap2 commit dc78b7c4b1780b316ca4aba2be247969ac1100ec-dirty
author | scottx611x |
---|---|
date | Fri, 27 Jul 2018 11:53:29 -0400 |
parents | 98944e914637 |
children | c53375e61b33 |
comparison
equal
deleted
inserted
replaced
40:98944e914637 | 41:e856fda6f375 |
---|---|
27 description="Generate Bam Quality Statistics" | 27 description="Generate Bam Quality Statistics" |
28 ) | 28 ) |
29 parser.add_argument('--input_file') | 29 parser.add_argument('--input_file') |
30 parser.add_argument('--out_genome_file') | 30 parser.add_argument('--out_genome_file') |
31 parser.add_argument('--out_dir') | 31 parser.add_argument('--out_dir') |
32 parser.add_argument('--out_zip') | |
33 parser.add_argument('--out_html') | |
34 parser.add_argument('--out_results') | |
32 parser.add_argument('--java_mem_size') | 35 parser.add_argument('--java_mem_size') |
33 | 36 |
34 args = parser.parse_args() | 37 args = parser.parse_args() |
35 print(args) | 38 print(args) |
36 | |
37 output_dir = os.path.dirname(args.out_genome_file) | |
38 | 39 |
39 qualimap_bamqc( | 40 qualimap_bamqc( |
40 args.input_file, | 41 args.input_file, |
41 args.out_genome_file, | 42 args.out_genome_file, |
42 args.out_dir, | 43 args.out_dir, |
44 ) | 45 ) |
45 | 46 |
46 shutil.make_archive( | 47 shutil.make_archive( |
47 'raw_data_qualimapReport', | 48 'raw_data_qualimapReport', |
48 'zip', | 49 'zip', |
49 args.out_dir + '/raw_data_qualimapReport' | 50 os.path.join(args.out_dir,'raw_data_qualimapReport') |
50 ) | 51 ) |
51 | 52 |
53 shutil.move("raw_data_qualimapReport.zip", args.out_zip) | |
52 shutil.move( | 54 shutil.move( |
53 "raw_data_qualimapReport.zip", | 55 os.path.join(args.out_dir, "genome_results.txt"), |
54 os.path.join(output_dir, "raw_data_qualimapReport.zip") | 56 args.out_results |
55 ) | 57 ) |
56 | 58 shutil.move( |
57 [ | 59 os.path.join(args.out_dir, "qualimapReport.html"), |
58 shutil.move( | 60 args.out_html |
59 os.path.join(args.out_dir, file_name), | 61 ) |
60 os.path.join(output_dir, file_name) | |
61 ) | |
62 for file_name in ["genome_results.txt", "qualimapReport.html"] | |
63 ] | |
64 | 62 |
65 if __name__ == "__main__": | 63 if __name__ == "__main__": |
66 main() | 64 main() |