Mercurial > repos > bornea > query_crapome

changeset 14:b451afff1aad draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bornea
date Mon, 25 Apr 2016 11:40:37 -0400
parents eb51efc01e35
children 0c0abdfad3b7
files CRAPomeQuery.py
diffstat 1 files changed, 4 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/CRAPomeQuery.py	Mon Apr 25 11:30:53 2016 -0400
+++ b/CRAPomeQuery.py	Mon Apr 25 11:40:37 2016 -0400
@@ -123,7 +123,7 @@
                         filt.append(j)
                         flag +=1
         if flag == 0: # if protein is not present in CRAPome database then add it
-            filt.append(["\t", "\t", i[0], "Invalid identifier / gene not available","\t","\t","\t"])
+            filt.append(["\t", "\t", i[0], "Invalid identifier / gene not available"])
     total = 0 # Experiment counter
     query = []
     for i in filt: # Create CRAPome file as list
@@ -135,7 +135,6 @@
             ave = []
             total = len(i[3:]) # calculate total experiments
             for j in i[3:]:
-                print j
                 if j != '0':
                     ave.append(int(j)) # calculate Ave.SC on only experiments with ID
                     cnt+=1
@@ -149,6 +148,9 @@
         else:
             temp.append(i[2]) # append accession
             temp.append(i[3])
+            temp.append("\t")
+            temp.append("\t")
+            temp.append("\t")
         query.append(temp) # final query results
 
     header = ["User Input","Mapped Gene Symbol","Num of Expt. (found/total)","Ave SC","Max SC"]