Repository 'query_crapome'
hg clone https://toolshed.g2.bx.psu.edu/repos/bornea/query_crapome

Changeset 14:b451afff1aad (2016-04-25)
Previous changeset 13:eb51efc01e35 (2016-04-25) Next changeset 15:0c0abdfad3b7 (2016-04-25)
Commit message:
Uploaded
modified:
CRAPomeQuery.py
b
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"]