comparison graphprot_predict_wrapper.py @ 4:4ad83aed5c3c draft

"planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/graphprot commit 15c382f71f1984aecff67ef9d268b992528996b2"
author bgruening
date Thu, 28 Jan 2021 15:06:14 +0000
parents ace92c9a4653
children ddcf35a868b8
comparison
equal deleted inserted replaced
3:ace92c9a4653 4:4ad83aed5c3c
219 c_in_fa = gplib.count_fasta_headers(args.in_fa) 219 c_in_fa = gplib.count_fasta_headers(args.in_fa)
220 assert c_in_fa, "input .fa file \"%s\" no headers found" % (args.in_fa) 220 assert c_in_fa, "input .fa file \"%s\" no headers found" % (args.in_fa)
221 print("# input .fa sequences: %i" % (c_in_fa)) 221 print("# input .fa sequences: %i" % (c_in_fa))
222 # Read in FASTA sequences to check for uppercase sequences. 222 # Read in FASTA sequences to check for uppercase sequences.
223 seqs_dic = gplib.read_fasta_into_dic(args.in_fa) 223 seqs_dic = gplib.read_fasta_into_dic(args.in_fa)
224 # Check for lowercase only sequences, which cause GP to crash.
225 error_mess = "input sequences encountered containing "\
226 "only lowercase characters or lowercase characters in between "\
227 "uppercase characters. Please provide either all uppercase "\
228 "sequences or sequences containing uppercase regions surrounded "\
229 "by lowercase context regions for structure calculation (see "\
230 "viewpoint concept in original GraphProt publication "\
231 "for more details)"
232 if args.ws_pred:
233 bad_ids = gplib.check_seqs_dic_format(seqs_dic)
234 assert not bad_ids, "%s" % (error_mess)
235
224 c_uc_nt = gplib.seqs_dic_count_uc_nts(seqs_dic) 236 c_uc_nt = gplib.seqs_dic_count_uc_nts(seqs_dic)
225 assert c_uc_nt, "no uppercase nucleotides in input .fa sequences. "\ 237 assert c_uc_nt, "no uppercase nucleotides in input .fa sequences. "\
226 "Please change sequences to uppercase (keep in mind "\ 238 "Please change sequences to uppercase (keep in mind "\
227 "GraphProt only scores uppercase regions (according "\ 239 "GraphProt only scores uppercase regions (according "\
228 "to its viewpoint concept)" 240 "to its viewpoint concept)"