Mercurial > repos > bornea > saint_preproc
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() |