annotate Vap.py @ 30:6669fd407dc9 draft default tip

Uploaded
author johnheap
date Fri, 07 Jun 2019 11:07:05 -0400
parents 26ec953069b3
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
20
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
1 """
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
2 * Copyright 2018 University of Liverpool
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
3 * Author: John Heap, Computational Biology Facility, UoL
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
4 * Based on original scripts of Sara Silva Pereira, Institute of Infection and Global Health, UoL
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
5 *
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
6 * Licensed under the Apache License, Version 2.0 (the "License");
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
7 * you may not use this file except in compliance with the License.
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
8 * You may obtain a copy of the License at
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
9 *
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
10 * http://www.apache.org/licenses/LICENSE-2.0
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
11 *
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
12 * Unless required by applicable law or agreed to in writing, software
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
13 * distributed under the License is distributed on an "AS IS" BASIS,
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
14 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
15 * See the License for the specific language governing permissions and
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
16 * limitations under the License.
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
17 *
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
18 """
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
19 #import subprocess
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
20 #import re
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
21 import os
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
22 import sys
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
23 #import pandas as pd
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
24 #import numpy as np
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
25 #import seaborn as sns
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
26 #import matplotlib.pyplot as plt
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
27 #from matplotlib.mlab import PCA
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
28 import Tryp_G
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
29 import Tryp_T
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
30 import Tryp_V
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
31 import Tryp_V_T
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
32 import argparse
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
33 #Entry .sort out the arguments
11
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
34
20
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
35 pdfExport = False
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
36 #parser = argparse.ArgumentParser(description='Variant Antigen Profiler - the VAP.')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
37 #parser.add_argument('name')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
38 #parser.add_argument('-t','-T', action = 'store_true', default = False, help = "Transciptomic Pathway")
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
39 #parser.add_argument('-p','-P', action = 'store_true', default = False, help = "Export PDFs to HTML directory")
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
40 #parser.add_argument('strain')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
41 #parser.add_argument('Forward_Read_File')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
42 #parser.add_argument('Reverse_Read_File')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
43 #parser.add_argument('htmlfile')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
44 #parser.add_argument('htmlresource')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
45 #parser.add_argument('heatmapFile')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
46 #parser.add_argument('PCAFile')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
47 #parser.add_argument('devheatmapFile')
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
48 #args = parser.parse_args()
11
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
49
20
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
50 #we have numerous parameters....
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
51 #hard code it for differnt types?
11
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
52
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
53
20
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
54 arguments = sys.argv
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
55 htmldir = arguments[len(arguments)-1] #last argument is always html_resource
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
56 if not os.path.exists(htmldir):
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
57 os.mkdir(htmldir)
11
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
58
20
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
59 if arguments[1] == 'g_assemble':
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
60 argdict = {'name':2, 'pdfexport':3, 'kmers':4,'inslen':5, 'covcut':6, 'forward':7, 'reverse':8, 'html_file':9, 'html_resource':10}
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
61 Tryp_G.assemble(arguments,argdict)
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
62 if arguments[1] == 'g_contigs':
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
63 argdict = {'name':2, 'pdfexport':3, 'contigs':4, 'html_file':5, 'html_resource':6}
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
64 Tryp_G.contigs(arguments,argdict)
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
65 if arguments[1] == 'transcipt':
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
66 argdict = {'name':2, 'pdfexport': 3, 'strain': 4, 'forward': 5, 'reverse': 6, 'html_file': 7, 'html_resource': 8}
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
67 Tryp_T.transcriptomicProcess(arguments,argdict)
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
68 if arguments[1] == 'v_assemble':
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
69 argdict = {'name':2, 'pdfexport':3, 'kmers':4,'inslen':5, 'covcut':6, 'forward':7, 'reverse':8, 'html_file':9, 'html_resource':10}
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
70 Tryp_V.vivax_assemble(arguments,argdict)
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
71 if arguments[1] == 'v_contigs':
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
72 argdict = {'name':2, 'pdfexport':3, 'contigs':4, 'html_file':5, 'html_resource':6}
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
73 Tryp_V.vivax_contigs(arguments,argdict)
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
74 if arguments[1] == 'v_transcript':
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
75 argdict = {'name':2, 'pdfexport': 3, 'refFastq': 4, 'forward': 5, 'reverse': 6, 'html_file': 7, 'html_resource': 8}
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
76 Tryp_V_T.transcriptomicProcess(arguments,argdict)
11
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
77
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
78
20
26ec953069b3 Uploaded
johnheap
parents: 11
diff changeset
79 sys.exit()
11
c4463e84a7ea Uploaded
johnheap
parents: 0
diff changeset
80