Mercurial > repos > dfornika > mentalist
view data_managers/data_manager_mentalist_build_db/data_manager/mentalist_build_db.py @ 0:1d9e3950ce61 draft
planemo upload for repository https://github.com/WGS-TB/MentaLiST/tree/master/galaxy commit 77d12d12820477cdb046be90175f5d2dab52e620
author | dfornika |
---|---|
date | Wed, 20 Dec 2017 17:36:03 -0500 |
parents | |
children |
line wrap: on
line source
#!/usr/bin/env python from __future__ import print_function import argparse import os import subprocess import sys import errno from json import dumps, loads DEFAULT_DATA_TABLE_NAMES = ["mentalist_databases"] def mentalist_build_db( data_manager_dict, database_name, kmer_size, profile, fasta_files, params, target_directory, data_table_names=DEFAULT_DATA_TABLE_NAMES ): args = [ 'mentalist', 'build_db', '--db', database_name, '-k', str(kmer_size)] if profile: args += ['--profile', profile] print(args) args += ['--fasta_files'] + fasta_files proc = subprocess.Popen( args=args, shell=False, cwd=target_directory ) return_code = proc.wait() if return_code: print("Error building database.", file=sys.stderr) sys.exit( return_code ) data_table_entry = dict( value=database_name, dbkey=database_name, name=database_name, path=database_name ) for data_table_name in data_table_names: _add_data_table_entry( data_manager_dict, data_table_name, data_table_entry ) def _add_data_table_entry( data_manager_dict, data_table_name, data_table_entry ): data_manager_dict['data_tables'] = data_manager_dict.get( 'data_tables', {} ) data_manager_dict['data_tables'][ data_table_name ] = data_manager_dict['data_tables'].get( data_table_name, [] ) data_manager_dict['data_tables'][ data_table_name ].append( data_table_entry ) return data_manager_dict def main(): parser = argparse.ArgumentParser() parser.add_argument('params') parser.add_argument( '-d', '--db', dest='database_name', default=None, help='Database Name' ) parser.add_argument( '-f', '--fasta_files', dest='fasta_files', nargs='+', default=None, help='FASTA Filenames' ) parser.add_argument( '-k', '--kmer_size', dest='kmer_size', type=int, default=None, help='kmer Size' ) parser.add_argument( '-p', '--profile', dest='profile', type=int, default=None, help='Profile' ) args = parser.parse_args() params = loads( open( args.params ).read() ) target_directory = params[ 'output_data' ][0]['extra_files_path'] try: os.mkdir( target_directory ) except OSError as exc: if exc.errno == errno.EEXIST and os.path.isdir( target_directory ): pass else: raise data_manager_dict = {} # build the index mentalist_build_db( data_manager_dict, args.database_name, args.kmer_size, args.profile, args.fasta_files, params, target_directory, DEFAULT_DATA_TABLE_NAMES ) # save info to json file open( args.params, 'wb' ).write( dumps( data_manager_dict ) ) if __name__ == "__main__": main()