Mercurial > repos > johnheap > vapper_galaxy
view Vap.py @ 10:320bdfa4d927 draft default tip
Deleted selected files
author | johnheap |
---|---|
date | Wed, 04 Jul 2018 11:37:20 -0400 |
parents | 08b35d28e826 |
children |
line wrap: on
line source
""" * Copyright 2018 University of Liverpool * Author: John Heap, Computational Biology Facility, UoL * Based on original scripts of Sara Silva Pereira, Institute of Infection and Global Health, UoL * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. * """ #import subprocess #import re import os import sys #import pandas as pd #import numpy as np #import seaborn as sns #import matplotlib.pyplot as plt #from matplotlib.mlab import PCA import Tryp_G import Tryp_T import Tryp_V import argparse #Entry .sort out the arguments pdfExport = False #parser = argparse.ArgumentParser(description='Variant Antigen Profiler - the VAP.') #parser.add_argument('name') #parser.add_argument('-t','-T', action = 'store_true', default = False, help = "Transciptomic Pathway") #parser.add_argument('-p','-P', action = 'store_true', default = False, help = "Export PDFs to HTML directory") #parser.add_argument('strain') #parser.add_argument('Forward_Read_File') #parser.add_argument('Reverse_Read_File') #parser.add_argument('htmlfile') #parser.add_argument('htmlresource') #parser.add_argument('heatmapFile') #parser.add_argument('PCAFile') #parser.add_argument('devheatmapFile') #args = parser.parse_args() #we have numerous parameters.... #hard code it for differnt types? arguments = sys.argv htmldir = arguments[len(arguments)-1] #last argument is always html_resource if not os.path.exists(htmldir): os.mkdir(htmldir) if arguments[1] == 'g_assemble': argdict = {'name':2, 'pdfexport':3, 'kmers':4,'inslen':5, 'covcut':6, 'forward':7, 'reverse':8, 'html_file':9, 'html_resource':10} Tryp_G.assemble(arguments,argdict) if arguments[1] == 'g_contigs': argdict = {'name':2, 'pdfexport':3, 'contigs':4, 'html_file':5, 'html_resource':6} Tryp_G.contigs(arguments,argdict) if arguments[1] == 'transcipt': argdict = {'name':2, 'pdfexport': 3, 'strain': 4, 'forward': 5, 'reverse': 6, 'html_file': 7, 'html_resource': 8} Tryp_T.transcriptomicProcess(arguments,argdict) if arguments[1] == 'v_assemble': argdict = {'name':2, 'pdfexport':3, 'kmers':4,'inslen':5, 'covcut':6, 'forward':7, 'reverse':8, 'html_file':9, 'html_resource':10} Tryp_V.vivax_assemble(arguments,argdict) if arguments[1] == 'v_contigs': argdict = {'name':2, 'pdfexport':3, 'contigs':4, 'html_file':5, 'html_resource':6} Tryp_V.vivax_contigs(arguments,argdict) sys.exit() #if not os.path.exists(args.htmlresource): # os.mkdir(args.htmlresource) #if args.p: # pdfExport = True #if args.t: # print ("Transcriptomic Pathway") # Tryp_T.transcriptomicProcess(args.name, pdfExport, args.strain, args.Forward_Read_File, args.Reverse_Read_File, args.htmlfile, args.htmlresource) #else: #print ('Genomic Pathway') #Tryp_G.test_output(args.name, args.htmlfile, args.htmlresource, args.heatmapFile, args.PCAFile,args.devheatmapFile) # Tryp_G.genomicProcess(args.name, pdfExport, args.Forward_Read_File, args.Reverse_Read_File, args.htmlfile, args.htmlresource)