# HG changeset patch
# User bornea
# Date 1461598837 14400
# Node ID b451afff1aadd94d688e97658954b8476f10cb5d
# Parent  eb51efc01e3595f647f16f94845ec2ddf57ae7e1
Uploaded

diff -r eb51efc01e35 -r b451afff1aad CRAPomeQuery.py
--- 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"]