|
0
|
1 #!/usr/bin/env python
|
|
|
2 import argparse
|
|
7
|
3 import os
|
|
0
|
4
|
|
1
|
5 import utils
|
|
0
|
6
|
|
|
7 OUTPUT_DIR = 'integratedGeneFamilies_dir'
|
|
|
8
|
|
|
9 parser = argparse.ArgumentParser()
|
|
13
|
10 parser.add_argument('--orthogroup_faa', dest='orthogroup_faa', help='Directory of input fasta datasets')
|
|
6
|
11 parser.add_argument('--scaffold', dest='scaffold', help='Orthogroups or gene families proteins scaffold')
|
|
0
|
12 parser.add_argument('--method', dest='method', help='Protein clustering method')
|
|
|
13 parser.add_argument('--orthogroup_fna', dest='orthogroup_fna', default=None, help='Use correspong coding sequences')
|
|
13
|
14 parser.add_argument('--output', dest='output', help='Output dataset')
|
|
|
15 parser.add_argument('--output_dir', dest='output_dir', help='Output dataset file_path directory')
|
|
0
|
16
|
|
|
17 args = parser.parse_args()
|
|
|
18
|
|
|
19 # Build the command line.
|
|
|
20 cmd = 'GeneFamilyIntegrator'
|
|
3
|
21 cmd += ' --orthogroup_faa %s' % args.orthogroup_faa
|
|
0
|
22 cmd += ' --scaffold %s' % args.scaffold
|
|
|
23 cmd += ' --method %s' % args.method
|
|
|
24 if args.orthogroup_fna is not None:
|
|
|
25 cmd += ' --orthogroup_fna'
|
|
13
|
26
|
|
0
|
27 # Run the command.
|
|
13
|
28 utils.run_command(cmd)
|
|
|
29
|
|
|
30 # Handle outputs.
|
|
6
|
31 utils.move_directory_files(os.path.join(OUTPUT_DIR, 'orthogroups_fasta'), args.output_dir)
|
|
0
|
32 utils.write_html_output(args.output, 'Integrated gene family sequences', args.output_dir)
|
