Mercurial > repos > bornea > saint_preprocessing

changeset 35:26cc583a4ae4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author bornea
date Thu, 19 May 2016 10:02:48 -0400
parents 05c5844e037b
children bc9c7764cc2f
files SAINT_preprocessing_mq_pep.py
diffstat 1 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/SAINT_preprocessing_mq_pep.py	Wed Apr 27 10:35:22 2016 -0400
+++ b/SAINT_preprocessing_mq_pep.py	Thu May 19 10:02:48 2016 -0400
@@ -128,7 +128,7 @@
         no_error_inter(MaxQuant_input)
         os.rename('prey.txt', sys.argv[5])
     elif prey == 'false':
-        if os.path.isfile('error proteins.txt') == True:
+        if os.path.isfile('./error_proteins.txt') == True:
             no_error_inter(MaxQuant_input)
         pass
     elif prey != 'true' or 'false':
@@ -139,7 +139,7 @@
 
 def get_info(uniprot_accession_in): 
     # Get aa lengths and gene name.
-    error = open('error proteins.txt', 'a+')
+    error = open('./error_proteins.txt', 'a+')
     data = open(fasta_db, 'r')
     data_lines = data.readlines()
     db_len = len(data_lines)
@@ -239,7 +239,7 @@
 
 def no_error_inter(MaxQuant_input):
     # Remake inter file without protein errors from Uniprot.
-    err = readtab("error proteins.txt")
+    err = readtab("./error_proteins.txt")
     bait = readtab(baitfile)
     data = read_MaxQuant(MaxQuant_input).data
     header = read_MaxQuant(MaxQuant_input).header