Mercurial > repos > abims-sbr > mutcount

annotate scripts/S01b_study_seq_composition_aa.py @ 9:04a9ada73cc4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
author abims-sbr
date Tue, 03 Jul 2018 10:55:46 -0400
parents f1e24200e5ae
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
1 #!/usr/bin/env python
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
2 # -*- coding: ascii -*-
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
3 ## Author: Eric FONTANILLAS
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
4 ## Date: 21.12.10
2
988467f963f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit cf1b9c905931ca2ca25faa4844d45c908756472f
abims-sbr
parents: 1
diff changeset
5 ## Last Version : 12/2017 by Victor Mataigne
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
6 ## Object: Test for compositional bias in genome and proteome as marker of thermal adaptation (comparison between 2 "hot" species: Ap and Ps and two "cold" species: Pg, Pp)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
7
7
f1e24200e5ae planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit b7a3030ea134b5dfad89b1a869db659d72d1145c
abims-sbr
parents: 5
diff changeset
8 import sys,os,shutil,subprocess,string, itertools
2
988467f963f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit cf1b9c905931ca2ca25faa4844d45c908756472f
abims-sbr
parents: 1
diff changeset
9 from functions import simplify_fasta_name, dico
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
10
2
988467f963f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit cf1b9c905931ca2ca25faa4844d45c908756472f
abims-sbr
parents: 1
diff changeset
11 script_path = os.path.dirname(sys.argv[0])
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
12
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
13 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
14 ###### DEF2 ######
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
15 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
16 def base_composition(seq):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
17 count_A=string.count(seq, "A")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
18 count_T=string.count(seq, "T")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
19 count_C=string.count(seq, "C")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
20 count_G=string.count(seq, "G")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
21
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
22
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
23 CG = count_C+count_G
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
24 AT = count_T+count_A
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
25
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
26 AG = count_A+count_G
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
27 TC = count_T+count_C
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
28
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
29 ## 1 ## Search for compositional bias in genome as marker of thermal adaptation: CG vs AT
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
30 ratio_CG_AT=float(CG)/float(AT)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
31
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
32 ## 2 ## Search for compositional bias in genome as marker of thermal adaptation: AG vs TC
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
33 ratio_purine_pyrimidine=float(AG)/float(TC)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
34
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
35 ## 3 ## Nucleotide proportion
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
36 ln = len(seq)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
37 prop_A = float(count_A)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
38 prop_T = float(count_T)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
39 prop_C = float(count_C)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
40 prop_G = float(count_G)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
41
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
42
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
43 return(ratio_CG_AT, ratio_purine_pyrimidine, prop_A, prop_T, prop_C, prop_G)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
44 ##############################################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
45
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
46
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
47 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
48 ###### DEF3 ######
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
49 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
50 def aa_composition1(seq):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
51
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
52 ## 1 ## count occurence of AA
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
53 count_K=string.count(seq,"K")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
54 count_R=string.count(seq,"R")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
55 count_A=string.count(seq,"A")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
56 count_F=string.count(seq,"F")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
57 count_I=string.count(seq,"I")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
58 count_L=string.count(seq,"L")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
59 count_M=string.count(seq,"M")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
60 count_V=string.count(seq,"V")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
61 count_W=string.count(seq,"W")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
62 count_N=string.count(seq,"N")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
63 count_Q=string.count(seq,"Q")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
64 count_S=string.count(seq,"S")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
65 count_T=string.count(seq,"T")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
66 count_H=string.count(seq,"H")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
67 count_Y=string.count(seq,"Y")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
68 count_C=string.count(seq,"C")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
69 count_D=string.count(seq,"D")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
70 count_E=string.count(seq,"E")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
71 count_P=string.count(seq,"P")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
72 count_G=string.count(seq,"G")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
73
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
74
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
75
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
76 ## 2 ## compute relative proportion
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
77 TOTAL=count_K+count_R+count_A+count_F+count_I+count_L+count_M+count_V+count_W+count_N+count_Q+count_S+count_T+count_H+count_Y+count_C+count_D+count_E+count_P+count_G
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
78 if (TOTAL!=0):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
79 ln = TOTAL
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
80
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
81 prop_K=float(count_K)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
82 prop_R=float(count_R)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
83 prop_A=float(count_A)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
84 prop_F=float(count_F)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
85 prop_I=float(count_I)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
86 prop_L=float(count_L)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
87 prop_M=float(count_M)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
88 prop_V=float(count_V)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
89 prop_W=float(count_W)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
90 prop_N=float(count_N)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
91 prop_Q=float(count_Q)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
92 prop_S=float(count_S)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
93 prop_T=float(count_T)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
94 prop_H=float(count_H)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
95 prop_Y=float(count_Y)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
96 prop_C=float(count_C)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
97 prop_D=float(count_D)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
98 prop_E=float(count_E)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
99 prop_P=float(count_P)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
100 prop_G=float(count_G)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
101 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
102 prop_K=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
103 prop_R=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
104 prop_A=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
105 prop_F=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
106 prop_I=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
107 prop_L=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
108 prop_M=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
109 prop_V=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
110 prop_W=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
111 prop_N=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
112 prop_Q=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
113 prop_S=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
114 prop_T=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
115 prop_H=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
116 prop_Y=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
117 prop_C=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
118 prop_D=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
119 prop_E=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
120 prop_P=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
121 prop_G=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
122
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
123
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
124
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
125
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
126 return(prop_K,prop_R,prop_A,prop_F,prop_I,prop_L,prop_M,prop_V,prop_W,prop_N,prop_Q,prop_S,prop_T,prop_H,prop_Y,prop_C,prop_D,prop_E,prop_P,prop_G)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
127
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
128 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
129 ###### DEF4 ######
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
130 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
131 def aa_composition2(seq):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
132
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
133 ## 1 ## count occurence of AA
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
134 count_K=string.count(seq,"K")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
135 count_R=string.count(seq,"R")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
136 count_A=string.count(seq,"A")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
137 count_F=string.count(seq,"F")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
138 count_I=string.count(seq,"I")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
139 count_L=string.count(seq,"L")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
140 count_M=string.count(seq,"M")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
141 count_V=string.count(seq,"V")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
142 count_W=string.count(seq,"W")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
143 count_N=string.count(seq,"N")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
144 count_Q=string.count(seq,"Q")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
145 count_S=string.count(seq,"S")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
146 count_T=string.count(seq,"T")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
147 count_H=string.count(seq,"H")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
148 count_Y=string.count(seq,"Y")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
149 count_C=string.count(seq,"C")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
150 count_D=string.count(seq,"D")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
151 count_E=string.count(seq,"E")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
152 count_P=string.count(seq,"P")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
153 count_G=string.count(seq,"G")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
154
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
155
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
156
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
157 ## 2 ## compute seq length
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
158 TOTAL=count_K+count_R+count_A+count_F+count_I+count_L+count_M+count_V+count_W+count_N+count_Q+count_S+count_T+count_H+count_Y+count_C+count_D+count_E+count_P+count_G
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
159 if (TOTAL!=0):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
160
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
161 ln = TOTAL
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
162 ##3 Famous Hyperthermophile Prokaryotes criterias
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
163
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
164 # 3.1. IVYWREL estimator => positivelly correlated with otpimal growth
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
165 count_IVYWREL = count_I+count_V+count_Y+count_W+count_R+count_E+count_L
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
166 prop_IVYWREL = float(count_IVYWREL)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
167
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
168 # 3.2. ERK estimator (i.e. ERK vs DNQTSHA) => positivelly correlated with optimal growth temperature
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
169 # ERK alone
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
170 count_ERK = count_E + count_R + count_K
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
171 prop_ERK = float(count_ERK)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
172 # DNQTSHA alone
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
173 count_DNQTSH = count_D+count_N+count_Q+count_T+count_S+count_H
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
174 prop_DNQTSH=float(count_DNQTSH)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
175 # ERK vs DNQTSH
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
176 if count_DNQTSH != 0:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
177 ratio_ERK_vs_DNQTSH = float(count_ERK)/float(count_DNQTSH)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
178 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
179 ratio_ERK_vs_DNQTSH=-1
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
180 # EK/QH estimator
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
181 count_EK = count_E+count_K
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
182 count_QH = count_Q+count_H
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
183
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
184 prop_EK = float(count_EK)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
185 prop_QH = float(count_QH)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
186
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
187 if count_QH != 0:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
188 ratio_EK_vs_QH = float(count_EK)/float(count_QH)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
189 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
190 ratio_EK_vs_QH=-1 ## "-1" will indicate the impossibility to compute the ratio (coz the numerator)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
191
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
192 ## 4 ## Mutationnal bias hypothesis => AT rich: favor FYMINK // GC rich: favor GARP
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
193 ## The mutational bias model predict a linear relationship between GARP vs FYMINK ==> so if outliers to that, it means that the excess of GARP or FYMINK are not explained by the mutationnal bias model but by other thing ... selection!!???
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
194 count_FYMINK=count_F+count_Y+count_M+count_I+count_N+count_K
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
195 prop_FYMINK = float(count_FYMINK)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
196
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
197 count_GARP=count_G+count_A+count_R+count_P
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
198 prop_GARP=float(count_GARP)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
199
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
200 ## 5 ## Hydophobicity hypothesis [should INCREASE with thermal adaptation]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
201 ## 5.1. AL
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
202 count_AVLIMFYW = count_A+count_V+count_L+count_I+count_F+count_Y+count_W+count_M
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
203 prop_AVLIMFYW=float(count_AVLIMFYW)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
204 ## 5.2. Only non-aromatic
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
205 count_AVLIM = count_A+count_V+count_L+count_I+count_M
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
206 prop_AVLIM=float(count_AVLIM)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
207 ## 5.3. Only aromatic (have they higher residus volume?? in such case opposite hypothesis based on residu volume, predict DECREASE for these aa in composition)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
208 count_FYW = count_F+count_Y+count_W
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
209 prop_FYW=float(count_FYW)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
210
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
211 ## 6 ## Charged hypothesis => positivelly correlated with optimal growth temperature
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
212 # All charged
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
213 count_RHKDE = count_R + count_H +count_K + count_D + count_E
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
214 prop_RHKDE = float(count_RHKDE)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
215 # Only positive
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
216 count_RHK = count_R + count_H +count_K
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
217 prop_RHK = float(count_RHK)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
218 # Only negative
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
219 count_DE = count_D + count_E
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
220 prop_DE = float(count_DE)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
221
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
222 ## 7 ## Neutral polar hypothesis [should DECREASE with thermal adaptation]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
223 count_STNQ = count_S+count_T+count_N+count_Q
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
224 prop_STNQ=float(count_STNQ)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
225
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
226
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
227 ## 9 ## PAYRE VS MGDS (FONTANILLAS CRITERIA)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
228 ## 9.1 ## Didier's criteria 1 = SMALL / BIG
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
229 count_PAYRE = count_A+count_Y+count_P+count_R+count_E
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
230 prop_PAYRE=float(count_PAYRE)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
231 count_MVGDS = count_V+count_M+count_S+count_G+count_D
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
232 prop_MVGDS=float(count_MVGDS)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
233 if count_MVGDS!= 0:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
234 ratio_PAYRE_vs_MVGDS = float(count_PAYRE)/float(count_MVGDS)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
235 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
236 ratio_PAYRE_vs_MVGDS=-1 ## "-1" will indicate the impossibility to compute the ratio (coz the numerator)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
237
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
238 ## 9.2 ## Didier's criteria 2 = VERY SMALL / BIG
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
239 count_AC = count_A+count_C
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
240 prop_AC=float(count_AC)/float(ln)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
241
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
242 #count_VLIM = count_V+count_L+count_I+count_M
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
243 if count_MVGDS != 0:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
244 ratio_AC_vs_MVGDS = float(count_AC)/float(count_MVGDS)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
245 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
246 ratio_AC_vs_MVGDS=-1 ## "-1" will indicate the impossibility to compute the ratio (coz the numerator)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
247 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
248 count_IVYWREL=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
249 prop_IVYWREL=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
250 count_ERK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
251 prop_ERK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
252 count_DNQTSH=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
253 prop_DNQTSH=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
254 ratio_ERK_vs_DNQTSH=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
255 count_EK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
256 prop_EK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
257 count_QH=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
258 prop_QH=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
259 ratio_EK_vs_QH=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
260 count_FYMINK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
261 prop_FYMINK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
262 count_GARP=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
263 prop_GARP=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
264 count_AVLIMFYW=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
265 prop_AVLIMFYW=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
266 count_AVLIM=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
267 prop_AVLIM=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
268 count_FYW=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
269 prop_FYW=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
270 count_STNQ=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
271 prop_STNQ=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
272 count_MVGDS=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
273 prop_MVGDS=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
274 count_PAYRE=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
275 prop_PAYRE=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
276 count_AC=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
277 prop_AC=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
278 ratio_PAYRE_vs_MVGDS=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
279 ratio_AC_vs_MVGDS=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
280 count_RHKDE=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
281 prop_RHKDE=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
282 count_RHK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
283 prop_RHK=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
284 count_DE=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
285 prop_DE=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
286
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
287 return(count_IVYWREL,prop_IVYWREL,count_ERK,prop_ERK,count_DNQTSH,prop_DNQTSH,ratio_ERK_vs_DNQTSH,count_EK,prop_EK,count_QH,prop_QH,ratio_EK_vs_QH,count_FYMINK,prop_FYMINK,count_GARP,prop_GARP,count_AVLIMFYW, prop_AVLIMFYW,count_AVLIM,prop_AVLIM,count_FYW,prop_FYW,count_STNQ, prop_STNQ, count_MVGDS,prop_MVGDS, count_PAYRE,prop_PAYRE, count_AC,prop_AC, ratio_PAYRE_vs_MVGDS, ratio_AC_vs_MVGDS, count_RHKDE,prop_RHKDE,count_RHK,prop_RHK,count_DE,prop_DE)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
288 #####################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
289
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
290
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
291 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
292 ###### DEF5 ######
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
293 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
294 def aa_properties(fileIN_aaProperties):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
295 next = fileIN_aaProperties.readline() ## JUMP HEADERS
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
296
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
297 bash_aa_properties={}
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
298
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
299 while 1:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
300 next = fileIN_aaProperties.readline()
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
301 if not next:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
302 break
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
303
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
304 S1 = string.split(next, ",")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
305
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
306 aa_name = S1[1]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
307 S2 = string.split(aa_name, "/")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
308 aa_code = S2[1][:-1]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
309
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
310 frequencies = S1[2][:-1]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
311 Residue_Weight = S1[5]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
312 Residue_Volume = S1[6]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
313 Partial_specific_volume = S1[7]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
314 Hydration = S1[8]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
315
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
316 bash_aa_properties[aa_code] = [frequencies,Residue_Weight,Residue_Volume,Partial_specific_volume,Hydration]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
317
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
318 return(bash_aa_properties)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
319
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
320
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
321 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
322 ###### DEF6 ######
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
323 ##################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
324 def sequence_properties_from_aa_properties(seq, bash_properties):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
325
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
326 ## 1 ## count occurence of AA
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
327 count_K=string.count(seq,"K")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
328 count_R=string.count(seq,"R")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
329 count_A=string.count(seq,"A")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
330 count_F=string.count(seq,"F")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
331 count_I=string.count(seq,"I")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
332 count_L=string.count(seq,"L")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
333 count_M=string.count(seq,"M")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
334 count_V=string.count(seq,"V")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
335 count_W=string.count(seq,"W")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
336 count_N=string.count(seq,"N")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
337 count_Q=string.count(seq,"Q")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
338 count_S=string.count(seq,"S")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
339 count_T=string.count(seq,"T")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
340 count_H=string.count(seq,"H")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
341 count_Y=string.count(seq,"Y")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
342 count_C=string.count(seq,"C")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
343 count_D=string.count(seq,"D")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
344 count_E=string.count(seq,"E")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
345 count_P=string.count(seq,"P")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
346 count_G=string.count(seq,"G")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
347
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
348 TOTAL=count_K+count_R+count_A+count_F+count_I+count_L+count_M+count_V+count_W+count_N+count_Q+count_S+count_T+count_H+count_Y+count_C+count_D+count_E+count_P+count_G
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
349
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
350 if (TOTAL!=0):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
351
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
352
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
353 ## 2 ## Compute properties 1: Residue Weight (Mr) (UNIT:Daltons):
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
354
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
355 Total_Residue_Weight = count_K*float(bash_properties["K"][1]) + count_R*float(bash_properties["R"][1]) + count_A*float(bash_properties["A"][1]) + count_F*float(bash_properties["F"][1]) + count_I*float(bash_properties["I"][1]) + count_L*float(bash_properties["L"][1]) + count_M*float(bash_properties["M"][1]) + count_V*float(bash_properties["V"][1]) + count_W*float(bash_properties["W"][1]) + count_N*float(bash_properties["N"][1]) + count_Q*float(bash_properties["Q"][1]) + count_S*float(bash_properties["S"][1]) + count_T*float(bash_properties["T"][1]) + count_H*float(bash_properties["H"][1]) + count_Y*float(bash_properties["Y"][1]) + count_C*float(bash_properties["C"][1]) + count_D*float(bash_properties["D"][1]) + count_E*float(bash_properties["E"][1]) + count_P*float(bash_properties["P"][1]) + count_G*float(bash_properties["G"][1])
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
356 Total_Residue_Volume = count_K*float(bash_properties["K"][2]) + count_R*float(bash_properties["R"][2]) + count_A*float(bash_properties["A"][2]) + count_F*float(bash_properties["F"][2]) + count_I*float(bash_properties["I"][2]) + count_L*float(bash_properties["L"][2]) + count_M*float(bash_properties["M"][2]) + count_V*float(bash_properties["V"][2]) + count_W*float(bash_properties["W"][2]) + count_N*float(bash_properties["N"][2]) + count_Q*float(bash_properties["Q"][2]) + count_S*float(bash_properties["S"][2]) + count_T*float(bash_properties["T"][2]) + count_H*float(bash_properties["H"][2]) + count_Y*float(bash_properties["Y"][2]) + count_C*float(bash_properties["C"][2]) + count_D*float(bash_properties["D"][2]) + count_E*float(bash_properties["E"][2]) + count_P*float(bash_properties["P"][2]) + count_G*float(bash_properties["G"][2])
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
357 Total_Partial_specific_volume = count_K*float(bash_properties["K"][3]) + count_R*float(bash_properties["R"][3]) + count_A*float(bash_properties["A"][3]) + count_F*float(bash_properties["F"][3]) + count_I*float(bash_properties["I"][3]) + count_L*float(bash_properties["L"][3]) + count_M*float(bash_properties["M"][3]) + count_V*float(bash_properties["V"][3]) + count_W*float(bash_properties["W"][3]) + count_N*float(bash_properties["N"][3]) + count_Q*float(bash_properties["Q"][3]) + count_S*float(bash_properties["S"][3]) + count_T*float(bash_properties["T"][3]) + count_H*float(bash_properties["H"][3]) + count_Y*float(bash_properties["Y"][3]) + count_C*float(bash_properties["C"][3]) + count_D*float(bash_properties["D"][3]) + count_E*float(bash_properties["E"][3]) + count_P*float(bash_properties["P"][3]) + count_G*float(bash_properties["G"][3])
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
358 Total_Hydration = count_K*float(bash_properties["K"][4]) + count_R*float(bash_properties["R"][4]) + count_A*float(bash_properties["A"][4]) + count_F*float(bash_properties["F"][4]) + count_I*float(bash_properties["I"][4]) + count_L*float(bash_properties["L"][4]) + count_M*float(bash_properties["M"][4]) + count_V*float(bash_properties["V"][4]) + count_W*float(bash_properties["W"][4]) + count_N*float(bash_properties["N"][4]) + count_Q*float(bash_properties["Q"][4]) + count_S*float(bash_properties["S"][4]) + count_T*float(bash_properties["T"][4]) + count_H*float(bash_properties["H"][4]) + count_Y*float(bash_properties["Y"][4]) + count_C*float(bash_properties["C"][4]) + count_D*float(bash_properties["D"][4]) + count_E*float(bash_properties["E"][4]) + count_P*float(bash_properties["P"][4]) + count_G*float(bash_properties["G"][4])
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
359 else:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
360 Total_Residue_Weight=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
361 Total_Residue_Volume=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
362 Total_Partial_specific_volume=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
363 Total_Hydration=0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
364
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
365 return(Total_Residue_Weight,Total_Residue_Volume,Total_Partial_specific_volume,Total_Hydration)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
366
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
367 ########################################################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
368
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
369
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
370
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
371 ###################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
372 ### RUN RUN RUN ###
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
373 ###################
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
374
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
375 ##Create specific folders
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
376 Path_IN_loci_NUC = "./IN_AA"
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
377 outpath= "./OUT"
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
378 os.makedirs(Path_IN_loci_NUC)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
379 os.makedirs(outpath)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
380
9
04a9ada73cc4 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
abims-sbr
parents: 7
diff changeset
381 infiles = []
04a9ada73cc4 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
abims-sbr
parents: 7
diff changeset
382 with open(sys.argv[2], 'r') as f:
04a9ada73cc4 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
abims-sbr
parents: 7
diff changeset
383 for line in f.readlines():
04a9ada73cc4 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
abims-sbr
parents: 7
diff changeset
384 infiles.append(line.strip('\n'))
04a9ada73cc4 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
abims-sbr
parents: 7
diff changeset
385
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
386 for file in infiles:
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
387 os.system("cp %s %s" %(file, Path_IN_loci_NUC))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
388
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
389 ## 1 ## List taxa
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
390 LT=[]
9
04a9ada73cc4 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
abims-sbr
parents: 7
diff changeset
391 cmd="grep '>' {}".format(sys.argv[1])
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
392 result = subprocess.check_output(cmd, shell=True)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
393 result=result.split('\n')
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
394 for i in result:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
395 sp=i[1:]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
396 if sp !='':
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
397 LT.append(sp)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
398 print LT
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
399
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
400
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
401 ## 2 ## PathIN
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
402 fileIN_properties = open("amino_acid_properties.csv", "r")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
403 Path_IN_loci_AA = "./IN_AA"
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
404 #Path_IN_loci_AA = "02_CDS_No_Missing_Data_aa_CDS_withM"
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
405 Lloci_AA = os.listdir(Path_IN_loci_AA)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
406
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
407 ## 3 ## PathOUT
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
408
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
409 ## 3.1 ## PROT composition
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
410 fileOUT_PROT_ALL=open("./OUT/prot_compositions_All_AA.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
411 fileOUT_PROT_ALL.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
412 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
413 fileOUT_PROT_ALL.write("%s_prop_K,%s_prop_R,%s_prop_A,%s_prop_F,%s_prop_I,%s_prop_L,%s_prop_M,%s_prop_V,%s_prop_W,%s_prop_N,%s_prop_Q,%s_prop_S,%s_prop_T,%s_prop_H,%s_prop_Y,%s_prop_C,%s_prop_D,%s_prop_E,%s_prop_P,%s_prop_G," %(taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
414 fileOUT_PROT_ALL.write("%s_prop_K,%s_prop_R,%s_prop_A,%s_prop_F,%s_prop_I,%s_prop_L,%s_prop_M,%s_prop_V,%s_prop_W,%s_prop_N,%s_prop_Q,%s_prop_S,%s_prop_T,%s_prop_H,%s_prop_Y,%s_prop_C,%s_prop_D,%s_prop_E,%s_prop_P,%s_prop_G" %(LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
415 fileOUT_PROT_ALL.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
416
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
417 ## 3.2 ## PROT IVYWREL
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
418 fileOUT_IVYWREL=open("./OUT/IVYWREL.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
419 fileOUT_IVYWREL.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
420 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
421 fileOUT_IVYWREL.write("%s_count_IVYWREL,%s_prop_IVYWREL," %(taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
422 fileOUT_IVYWREL.write("%s_count_IVYWREL,%s_prop_IVYWREL" %(LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
423 fileOUT_IVYWREL.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
424
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
425 ## 3.3 ## PROT ERK_DNQTSHA
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
426 fileOUT_ERK_DNQTSH=open("./OUT/ERK_DNQTSH.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
427 fileOUT_ERK_DNQTSH.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
428 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
429 fileOUT_ERK_DNQTSH.write("%s_count_ERK,%s_prop_ERK,%s_count_DNQTSH,%s_prop_DNQTSH,%s_ratio_ERK_vs_DNQTSH," %(taxa,taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
430 fileOUT_ERK_DNQTSH.write("%s_count_ERK,%s_prop_ERK,%s_count_DNQTSH,%s_prop_DNQTSH,%s_ratio_ERK_vs_DNQTSH" %(LT[-1],LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
431 fileOUT_ERK_DNQTSH.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
432
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
433 ## 3.4 ## PROT EK_QH
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
434 fileOUT_EK_QH=open("./OUT/EK_QH.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
435 fileOUT_EK_QH.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
436 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
437 fileOUT_EK_QH.write("%s_count_EK,%s_prop_EK,%s_count_QH,%s_prop_QH,%s_ratio_EK_vs_QH," %(taxa,taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
438 fileOUT_EK_QH.write("%s_count_EK,%s_prop_EK,%s_count_QH,%s_prop_QH,%s_ratio_EK_vs_QH" %(LT[-1],LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
439 fileOUT_EK_QH.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
440
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
441 ## 3.5 ## PROT FYMINK_GARP
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
442 fileOUT_FYMINK_GARP=open("./OUT/FYMINK_GARP.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
443 fileOUT_FYMINK_GARP.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
444 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
445 fileOUT_FYMINK_GARP.write("%s_count_FYMINK,%s_prop_FYMINK,%s_count_GARP,%s_prop_GARP," %(taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
446 fileOUT_FYMINK_GARP.write("%s_count_FYMINK,%s_prop_FYMINK,%s_count_GARP,%s_prop_GARP" %(LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
447 fileOUT_FYMINK_GARP.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
448
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
449 ## 3.6 ## PROT AVLIMFYW
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
450 fileOUT_AVLIMFYW=open("./OUT/AVLIMFYW.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
451 fileOUT_AVLIMFYW.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
452 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
453 fileOUT_AVLIMFYW.write("%s_count_AVLIMFYW,%s_prop_AVLIMFYW,%s_count_AVLIM,%s_prop_AVLIM,%s_count_FYW,%s_prop_FYW," %(taxa,taxa,taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
454 fileOUT_AVLIMFYW.write("%s_count_AVLIMFYW,%s_prop_AVLIMFYW,%s_count_AVLIM,%s_prop_AVLIM,%s_count_FYW,%s_prop_FYW" %(LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
455 fileOUT_AVLIMFYW.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
456
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
457 ## 3.7 ## PROT STNQ
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
458 fileOUT_STNQ=open("./OUT/STNQ.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
459 fileOUT_STNQ.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
460 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
461 fileOUT_STNQ.write("%s_count_STNQ,%s_prop_STNQ," %(taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
462 fileOUT_STNQ.write("%s_count_STNQ,%s_prop_STNQ" %(LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
463 fileOUT_STNQ.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
464
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
465 ## 3.8 ## PROT RHKDE
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
466 fileOUT_RHKDE=open("./OUT/RHKDE.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
467 fileOUT_RHKDE.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
468 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
469 fileOUT_RHKDE.write("%s_count_RHKDE,%s_prop_RHKDE,%s_count_RHK,%s_prop_RHK,%s_count_DE,%s_prop_DE," %(taxa,taxa,taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
470 fileOUT_RHKDE.write("%s_count_RHKDE,%s_prop_RHKDE,%s_count_RHK,%s_prop_RHK,%s_count_DE,%s_prop_DE" %(LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
471 fileOUT_RHKDE.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
472
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
473 ## 3.9 ## PROT DIDER CRITERIA
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
474 fileOUT_PAYRE=open("./OUT/PAYRE-MVGDS.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
475 fileOUT_PAYRE.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
476 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
477 fileOUT_PAYRE.write("%s_count_PAYRE,%s_prop_PAYRE,%s_count_AC,%s_prop_AC,%s_count_MVGDS,%s_prop_MVGDS,%s_ratio_PAYRE_vs_MVGDS,%s_ratio_AC_vs_MVGDS," %(taxa,taxa,taxa,taxa,taxa,taxa,taxa,taxa))
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
478 fileOUT_PAYRE.write("%s_count_PAYRE,%s_prop_PAYRE,%s_count_AC,%s_prop_AC,%s_count_MVGDS,%s_prop_MVGDS,%s_ratio_PAYRE_vs_MVGDS,%s_ratio_AC_vs_MVGDS" %(LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1],LT[-1]))
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
479 fileOUT_PAYRE.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
480
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
481 ## 3.10 ## PROT Total residue weight
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
482 fileOUT_TotalResidueWeight=open("./OUT/TotalResidueWeight.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
483 fileOUT_TotalResidueWeight.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
484 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
485 fileOUT_TotalResidueWeight.write("%s_Total_Residue_Weight," %taxa)
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
486 fileOUT_TotalResidueWeight.write("%s_Total_Residue_Weight" %LT[-1])
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
487 fileOUT_TotalResidueWeight.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
488
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
489 ## 3.11 ## PROT Total residue volume
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
490 fileOUT_TotalResidueVolume=open("./OUT/TotalResidueVolume.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
491 fileOUT_TotalResidueVolume.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
492 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
493 fileOUT_TotalResidueVolume.write("%s_Total_Residue_Volume," %taxa)
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
494 fileOUT_TotalResidueVolume.write("%s_Total_Residue_Volume" %LT[-1])
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
495 fileOUT_TotalResidueVolume.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
496
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
497 ## 3.12 ## PROT Total partial specific volume
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
498 fileOUT_TotalPartialSpecificVolume=open("./OUT/TotalPartialSpecificVolume.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
499 fileOUT_TotalPartialSpecificVolume.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
500 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
501 fileOUT_TotalPartialSpecificVolume.write("%s_Total_Partial_Specific_Volume," %taxa)
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
502 fileOUT_TotalPartialSpecificVolume.write("%s_Total_Partial_Specific_Volume" %LT[-1])
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
503 fileOUT_TotalPartialSpecificVolume.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
504
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
505 ## 3.13 ## PROT Total hydratation
1
8de21b6eb110 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 44a89d5eeb82789bfc643b33c11f391281b6374b
abims-sbr
parents: 0
diff changeset
506 fileOUT_TotalHydratation=open("./OUT/TotalHydratation.csv","w")
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
507 fileOUT_TotalHydratation.write("LOCUS,")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
508 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
509 fileOUT_TotalHydratation.write("%s_Total_Hydratation," %taxa)
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
510 fileOUT_TotalHydratation.write("%s_Total_Hydratation" %LT[-1])
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
511 fileOUT_TotalHydratation.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
512
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
513 #####################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
514 ## 4 ## Process Loci
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
515 #####################
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
516 bash_aa_properties = aa_properties(fileIN_properties)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
517
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
518 for locus in Lloci_AA:
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
519 print locus
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
520 path_locus = "%s/%s" %(Path_IN_loci_AA, locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
521 bash = dico(path_locus,LT)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
522
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
523 #print bash
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
524
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
525 fileOUT_PROT_ALL.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
526 fileOUT_IVYWREL.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
527 fileOUT_ERK_DNQTSH.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
528 fileOUT_EK_QH.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
529 fileOUT_FYMINK_GARP.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
530 fileOUT_AVLIMFYW.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
531 fileOUT_STNQ.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
532 fileOUT_RHKDE.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
533 fileOUT_PAYRE.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
534 fileOUT_TotalResidueWeight.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
535 fileOUT_TotalResidueVolume.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
536 fileOUT_TotalPartialSpecificVolume.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
537 fileOUT_TotalHydratation.write("%s," %locus)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
538
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
539 for taxa in LT[0:-1]:
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
540 if taxa in bash.keys():
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
541 seq = bash[taxa]
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
542 prop_K,prop_R,prop_A,prop_F,prop_I,prop_L,prop_M,prop_V,prop_W,prop_N,prop_Q,prop_S,prop_T,prop_H,prop_Y,prop_C,prop_D,prop_E,prop_P,prop_G = aa_composition1(seq) ### DEF3 ###
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
543 count_IVYWREL,prop_IVYWREL,count_ERK,prop_ERK,count_DNQTSH,prop_DNQTSH,ratio_ERK_vs_DNQTSH,count_EK,prop_EK,count_QH,prop_QH,ratio_EK_vs_QH,count_FYMINK,prop_FYMINK,count_GARP,prop_GARP,count_AVLIMFYW,prop_AVLIMFYW,count_AVLIM,prop_AVLIM,count_FYW,prop_FYW,count_STNQ,prop_STNQ, count_MVGDS,prop_MVGDS, count_PAYRE,prop_PAYRE, count_AC,prop_AC, ratio_PAYRE_vs_MVGDS, ratio_AC_vs_MVGDS,count_RHKDE,prop_RHKDE,count_RHK,prop_RHK,count_DE,prop_DE = aa_composition2(seq) ### DEF4 ###
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
544 Total_Residue_Weight,Total_Residue_Volume,Total_Partial_Specific_Volume,Total_Hydration = sequence_properties_from_aa_properties(seq, bash_aa_properties) ### DEF6 ###
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
545
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
546 fileOUT_PROT_ALL.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f," %(prop_K,prop_R,prop_A,prop_F,prop_I,prop_L,prop_M,prop_V,prop_W,prop_N,prop_Q,prop_S,prop_T,prop_H,prop_Y,prop_C,prop_D,prop_E,prop_P,prop_G))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
547 fileOUT_IVYWREL.write("%.5f,%.5f," %(count_IVYWREL, prop_IVYWREL))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
548 fileOUT_ERK_DNQTSH.write("%.5f,%.5f,%.5f,%.5f,%.5f," %(count_ERK,prop_ERK,count_DNQTSH,prop_DNQTSH,ratio_ERK_vs_DNQTSH))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
549 fileOUT_EK_QH.write("%.5f,%.5f,%.5f,%.5f,%.5f," %(count_EK,prop_EK,count_QH,prop_QH,ratio_EK_vs_QH))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
550 fileOUT_FYMINK_GARP.write("%.5f,%.5f,%.5f,%.5f," %(count_FYMINK,prop_FYMINK,count_GARP,prop_GARP))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
551 fileOUT_AVLIMFYW.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f," %(count_AVLIMFYW,prop_AVLIMFYW,count_AVLIM,prop_AVLIM,count_FYW,prop_FYW))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
552 fileOUT_STNQ.write("%.5f,%.5f," %(count_STNQ,prop_STNQ))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
553 fileOUT_RHKDE.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,"%(count_RHKDE,prop_RHKDE,count_RHK,prop_RHK,count_DE,prop_DE))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
554 fileOUT_PAYRE.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f," %(count_PAYRE,prop_PAYRE,count_AC,prop_AC,count_MVGDS,prop_MVGDS,ratio_PAYRE_vs_MVGDS,ratio_AC_vs_MVGDS))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
555 fileOUT_TotalResidueWeight.write("%.5f," %Total_Residue_Weight)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
556 fileOUT_TotalResidueVolume.write("%.5f," %Total_Residue_Volume)
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
557 fileOUT_TotalPartialSpecificVolume.write("%.5f," %(Total_Partial_Specific_Volume))
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
558 fileOUT_TotalHydratation.write("%.5f," % Total_Hydration)
2
988467f963f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit cf1b9c905931ca2ca25faa4844d45c908756472f
abims-sbr
parents: 1
diff changeset
559 else:
4
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
560 fileOUT_PROT_ALL.write("%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s," %("NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
561 fileOUT_IVYWREL.write("%s,%s," %("NA", "NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
562 fileOUT_ERK_DNQTSH.write("%s,%s,%s,%s,%s," %("NA","NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
563 fileOUT_EK_QH.write("%s,%s,%s,%s,%s," %("NA","NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
564 fileOUT_FYMINK_GARP.write("%s,%s,%s,%s," %("NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
565 fileOUT_AVLIMFYW.write("%s,%s,%s,%s,%s,%s," %("NA","NA","NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
566 fileOUT_STNQ.write("%s,%s," %("NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
567 fileOUT_RHKDE.write("%s,%s,%s,%s,%s,%s,"%("NA","NA","NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
568 fileOUT_PAYRE.write("%s,%s,%s,%s,%s,%s,%s,%s," %("NA","NA","NA","NA","NA","NA","NA","NA"))
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
569 fileOUT_TotalResidueWeight.write("%s," %"NA")
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
570 fileOUT_TotalResidueVolume.write("%s," %"NA")
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
571 fileOUT_TotalPartialSpecificVolume.write("%s," %"NA")
5766f80370e7 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 17acd02c547bd1f7661a846661aa99de9743efe9
abims-sbr
parents: 2
diff changeset
572 fileOUT_TotalHydratation.write("%s," %"NA")
5
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
573
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
574 if LT[-1] in bash.keys():
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
575 seq = bash[LT[-1]]
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
576 prop_K,prop_R,prop_A,prop_F,prop_I,prop_L,prop_M,prop_V,prop_W,prop_N,prop_Q,prop_S,prop_T,prop_H,prop_Y,prop_C,prop_D,prop_E,prop_P,prop_G = aa_composition1(seq) ### DEF3 ###
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
577 count_IVYWREL,prop_IVYWREL,count_ERK,prop_ERK,count_DNQTSH,prop_DNQTSH,ratio_ERK_vs_DNQTSH,count_EK,prop_EK,count_QH,prop_QH,ratio_EK_vs_QH,count_FYMINK,prop_FYMINK,count_GARP,prop_GARP,count_AVLIMFYW,prop_AVLIMFYW,count_AVLIM,prop_AVLIM,count_FYW,prop_FYW,count_STNQ,prop_STNQ, count_MVGDS,prop_MVGDS, count_PAYRE,prop_PAYRE, count_AC,prop_AC, ratio_PAYRE_vs_MVGDS, ratio_AC_vs_MVGDS,count_RHKDE,prop_RHKDE,count_RHK,prop_RHK,count_DE,prop_DE = aa_composition2(seq) ### DEF4 ###
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
578 Total_Residue_Weight,Total_Residue_Volume,Total_Partial_Specific_Volume,Total_Hydration = sequence_properties_from_aa_properties(seq, bash_aa_properties) ### DEF6 ###
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
579
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
580 fileOUT_PROT_ALL.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f" %(prop_K,prop_R,prop_A,prop_F,prop_I,prop_L,prop_M,prop_V,prop_W,prop_N,prop_Q,prop_S,prop_T,prop_H,prop_Y,prop_C,prop_D,prop_E,prop_P,prop_G))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
581 fileOUT_IVYWREL.write("%.5f,%.5f" %(count_IVYWREL, prop_IVYWREL))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
582 fileOUT_ERK_DNQTSH.write("%.5f,%.5f,%.5f,%.5f,%.5f" %(count_ERK,prop_ERK,count_DNQTSH,prop_DNQTSH,ratio_ERK_vs_DNQTSH))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
583 fileOUT_EK_QH.write("%.5f,%.5f,%.5f,%.5f,%.5f" %(count_EK,prop_EK,count_QH,prop_QH,ratio_EK_vs_QH))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
584 fileOUT_FYMINK_GARP.write("%.5f,%.5f,%.5f,%.5f" %(count_FYMINK,prop_FYMINK,count_GARP,prop_GARP))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
585 fileOUT_AVLIMFYW.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f" %(count_AVLIMFYW,prop_AVLIMFYW,count_AVLIM,prop_AVLIM,count_FYW,prop_FYW))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
586 fileOUT_STNQ.write("%.5f,%.5f" %(count_STNQ,prop_STNQ))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
587 fileOUT_RHKDE.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f"%(count_RHKDE,prop_RHKDE,count_RHK,prop_RHK,count_DE,prop_DE))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
588 fileOUT_PAYRE.write("%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f,%.5f" %(count_PAYRE,prop_PAYRE,count_AC,prop_AC,count_MVGDS,prop_MVGDS,ratio_PAYRE_vs_MVGDS,ratio_AC_vs_MVGDS))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
589 fileOUT_TotalResidueWeight.write("%.5f" %Total_Residue_Weight)
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
590 fileOUT_TotalResidueVolume.write("%.5f" %Total_Residue_Volume)
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
591 fileOUT_TotalPartialSpecificVolume.write("%.5f" %(Total_Partial_Specific_Volume))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
592 fileOUT_TotalHydratation.write("%.5f" % Total_Hydration)
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
593 else:
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
594 fileOUT_PROT_ALL.write("%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s,%s" %("NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
595 fileOUT_IVYWREL.write("%s,%s" %("NA", "NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
596 fileOUT_ERK_DNQTSH.write("%s,%s,%s,%s,%s" %("NA","NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
597 fileOUT_EK_QH.write("%s,%s,%s,%s,%s" %("NA","NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
598 fileOUT_FYMINK_GARP.write("%s,%s,%s,%s" %("NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
599 fileOUT_AVLIMFYW.write("%s,%s,%s,%s,%s,%s" %("NA","NA","NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
600 fileOUT_STNQ.write("%s,%s" %("NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
601 fileOUT_RHKDE.write("%s,%s,%s,%s,%s,%s"%("NA","NA","NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
602 fileOUT_PAYRE.write("%s,%s,%s,%s,%s,%s,%s,%s" %("NA","NA","NA","NA","NA","NA","NA","NA"))
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
603 fileOUT_TotalResidueWeight.write("%s" %"NA")
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
604 fileOUT_TotalResidueVolume.write("%s" %"NA")
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
605 fileOUT_TotalPartialSpecificVolume.write("%s" %"NA")
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
606 fileOUT_TotalHydratation.write("%s" %"NA")
0ba551449008 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 273a9af69b672b2580cd5dec4c0e67a4a96fb0fe
abims-sbr
parents: 4
diff changeset
607
0
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
608 ## END LINE
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
609 fileOUT_PROT_ALL.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
610 fileOUT_IVYWREL.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
611 fileOUT_ERK_DNQTSH.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
612 fileOUT_EK_QH.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
613 fileOUT_FYMINK_GARP.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
614 fileOUT_AVLIMFYW.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
615 fileOUT_STNQ.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
616 fileOUT_RHKDE.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
617 fileOUT_PAYRE.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
618 fileOUT_TotalResidueWeight.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
619 fileOUT_TotalResidueVolume.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
620 fileOUT_TotalPartialSpecificVolume.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
621 fileOUT_TotalHydratation.write("\n")
78dd6454f6f0 planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit 73670b26c75bb6c1a6332481920f3036314de364
abims-sbr
parents:
diff changeset
622