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