comparison SAINT_preprocessing_v6_mq_pep.py @ 22:729aac64ce43 draft

Uploaded
author bornea
date Thu, 19 Nov 2015 13:38:31 -0500
parents 9c3d2b1d0c1c
children
comparison
equal deleted inserted replaced
21:9e0a894d2676 22:729aac64ce43
15 # 1) infile: maxquant "Samples Report" output 15 # 1) infile: maxquant "Samples Report" output
16 # 2) baitfile: SAINT formatted bait file generated in Galaxy 16 # 2) baitfile: SAINT formatted bait file generated in Galaxy
17 # 3) prey: Y or N for generating a prey file (requires internet connection) 17 # 3) prey: Y or N for generating a prey file (requires internet connection)
18 ####################################################################################### 18 #######################################################################################
19 mq_file = sys.argv[1] 19 mq_file = sys.argv[1]
20 cmd = r"Rscript /home/bornea/galaxy_moffitt_dev/tools/Moffitt_Tools/bubblebeam/pre_process_protein_name_set.R " + str(mq_file) 20 ins_path = sys.argv[8]
21 names_path = str(ins_path) + r"uniprot_names.txt"
22 cmd = r"Rscript "+ str(ins_path) +"pre_process_protein_name_set.R " + str(mq_file) + " " + str(names_path)
21 os.system(cmd) 23 os.system(cmd)
22 24
23 infile = "./tukeys_output.txt" #maxquant "Samples Report" output 25 infile = "./tukeys_output.txt" #maxquant "Samples Report" output
24 prey = sys.argv[2] # Y or N 26 prey = sys.argv[2] # Y or N
25 fasta_db = sys.argv[3] 27 fasta_db = sys.argv[3]
26 if fasta_db == "None": 28 if fasta_db == "None":
27 fasta_db = "/home/bornea/galaxy_moffitt_dev/tools/Moffitt_Tools/bubblebeam/SwissProt_HUMAN_2014_08.fasta" 29 fasta_db = str(ins_path) + "SwissProt_HUMAN_2014_08.fasta"
28 make_bait= sys.argv[6] 30 make_bait= sys.argv[6]
29 31
30 def bait_create(baits, infile): 32 def bait_create(baits, infile):
31 #Takes the Bait specified by the user and makes them into a Bait file and includes a check to make sure they are using valid baits. 33 #Takes the Bait specified by the user and makes them into a Bait file and includes a check to make sure they are using valid baits.
32 baits = make_bait.split() 34 baits = make_bait.split()