Mercurial > repos > rnateam > graphprot_predict_profile
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)" |