|
0
|
1 #!/usr/bin/env python
|
|
|
2 import argparse
|
|
|
3 import subprocess
|
|
|
4
|
|
1
|
5 import utils
|
|
0
|
6
|
|
|
7 OUTPUT_DIR = 'integratedGeneFamilies_dir'
|
|
|
8
|
|
|
9 parser = argparse.ArgumentParser()
|
|
|
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')
|
|
|
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")
|
|
|
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'
|
|
|
26 # Run the command.
|
|
|
27 proc = subprocess.Popen(args=cmd, stderr=subprocess.PIPE, stdout=subprocess.PIPE, shell=True)
|
|
|
28 rc = proc.wait()
|
|
|
29 utils.check_execution_errors(rc, proc.stderr)
|
|
6
|
30 utils.move_directory_files(os.path.join(OUTPUT_DIR, 'orthogroups_fasta'), args.output_dir)
|
|
0
|
31 utils.write_html_output(args.output, 'Integrated gene family sequences', args.output_dir)
|
