Mercurial > repos > charles_s_test > seqsero2
annotate libs/Initial_functions.py @ 16:ed79f56d2184 draft default tip
planemo upload commit a4fb57231f274270afbfebd47f67df05babffa4a-dirty
| author | charles_s_test |
|---|---|
| date | Thu, 30 Nov 2017 16:58:25 -0500 |
| parents | 6895de35a263 |
| children |
| rev | line source |
|---|---|
|
0
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
1 #!/usr/bin/env python |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
2 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
3 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
4 def To_list(L): |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
5 import string |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
6 New_list=[] |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
7 for x in L: |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
8 x1=x[:-1] |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
9 x1=string.atoi(x1) |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
10 New_list.append(x1) |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
11 return New_list |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
12 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
13 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
14 def Uniq(L): #return the uniq list and the count number |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
15 Old=L |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
16 L.sort() |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
17 L = [L[i] for i in range(len(L)) if L[i] not in L[:i]] |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
18 count=[] |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
19 for j in range(len(L)): |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
20 y=0 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
21 for x in Old: |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
22 if L[j]==x: |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
23 y+=1 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
24 count.append(y) |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
25 return (L,count) |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
26 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
27 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
28 def Parse_seros_in_genome_trakr(L): #return the sero names in the sra_result.xlsx, the next step is usually "Uniq" in above |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
29 names2=[] |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
30 for x in L: |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
31 if "serovar" in x: |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
32 key_word=x.split("serovar")[1].split("_")[1] #the seronames |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
33 if key_word!="str." and key_word!="group": #to eliminate some "serovar_str." and "serovar_group" |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
34 if key_word in ["I","II","III","IIIa","IIIb","IV","VI","B"]: #the serovar is behind those letters |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
35 if x.split("serovar")[1].split("_")[2]!="str.": |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
36 names2.append(x.split("serovar")[1].split("_")[2]) |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
37 else: |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
38 names2.append(x.split("serovar")[1].split("_")[1]) |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
39 return names2 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
40 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
41 |
|
6895de35a263
planemo upload commit 844a891e4eaf732830043204ac636907eefb011d-dirty
charles_s_test
parents:
diff
changeset
|
42 |
