comparison graphprot_predict_wrapper.py @ 4:2d6de5b769e6 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:05:38 +0000
parents 9a83a84a25a7
children 58ebf089377e
comparison
equal deleted inserted replaced
3:9a83a84a25a7 4:2d6de5b769e6
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)"