# HG changeset patch # User greg # Date 1495462310 14400 # Node ID 9a105e8f08e55d100ab32a690255d81e712cd7bf # Parent 0f65703d5bb8c94e895c179b0d80f87f3464245b Uploaded diff -r 0f65703d5bb8 -r 9a105e8f08e5 assembly_post_processor.py --- a/assembly_post_processor.py Mon May 15 08:26:03 2017 -0400 +++ b/assembly_post_processor.py Mon May 22 10:11:50 2017 -0400 @@ -1,8 +1,11 @@ #!/usr/bin/env python import argparse +import os import utils +OUTPUT_DIR = 'assemblyPostProcessing_dir' + parser = argparse.ArgumentParser() parser.add_argument('--dereplicate', dest='dereplicate', default=None, help='Remove duplicate sequences') parser.add_argument('--gap_trimming', dest='gap_trimming', type=float, default=0, help='Trim alignments') @@ -10,6 +13,8 @@ parser.add_argument('--method', dest='method', default=None, help='Protein clustering method') parser.add_argument('--min_length', dest='min_length', type=int, default=0, help='Minimum sequence length') parser.add_argument('--num_threads', dest='num_threads', type=int, help='Number of processors') +parser.add_argument('--output_pttgf', dest='output_pttgf', default=None, help='Primary targeted gene families dataset') +parser.add_argument('--output_pttgf_dir', dest='output_pttgf_dir', default=None, help='Directory hierarchy of targeted gene family datasets') parser.add_argument('--prediction_method', dest='prediction_method', help='Coding regions prediction method') parser.add_argument('--scaffold', dest='scaffold', default=None, help='Gene family scaffold') parser.add_argument('--score_matrices', dest='score_matrices', default=None, help='Scores matrices') @@ -41,3 +46,8 @@ cmd += ' --transcripts %s' % args.transcripts # Run the command. utils.run_command(cmd) + +# Handle outputs. +if args.output_pttgf is not None and args.output_pttgf_dir is not None: + src_output_dir = os.path.join(OUTPUT_DIR, 'targeted_gene_families') + utils.move_directory_files(src_output_dir, args.output_pttgf_dir)