LSHTM_analysis/mcsm_na/run_submit_mcsm_na.py

49 lines
1.5 KiB
Python
Executable file

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Fri Feb 12 12:15:26 2021
@author: tanu
"""
#%% load packages
import os
homedir = os.path.expanduser('~')
os.chdir (homedir + '/git/LSHTM_analysis/mcsm_na')
from submit_mcsm_na import *
########################################################################
# variables
my_host = 'http://biosig.unimelb.edu.au'
my_prediction_url = f"{my_host}/mcsm_na/run_prediction_list"
print(my_prediction_url)
# TODO: add cmd line args
#gene = 'gid'
drug = ''
datadir = homedir + '/git/Data/'
indir = datadir + drug + 'input/'
outdir = datadir + drug + 'output/'
outdir_mcsm_na = outdir + 'mcsm_na_results/'
my_nuc_type = 'RNA'
my_pdb_file = indir + gene.lower() + '_complex.pdb'
#=============================================================================
# batch 26: 25.txt # RAN: 16 Feb:
# batch 27: 26.txt # RAN: 6 Aug:
# off by one
my_mutation_list = outdir + '/snp_batches/20/snp_batch_26.txt'
my_suffix = 'gid_b27'
#==============================================================================
#==========================
# CALL: submit_mcsm_na()
# Data: gid+streptomycin
#==========================
submit_mcsm_na(host_url = my_host
, pdb_file = my_pdb_file
, mutation_list = my_mutation_list
, nuc_type = my_nuc_type
, prediction_url = my_prediction_url
, output_dir = outdir_mcsm_na
, outfile_suffix = my_suffix)
#%%#####################################################################