Mercurial > repos > rnateam > rnacommender

diff rbpfeatures.py @ 4:a609d6dc8047 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/bgruening/galaxytools/tree/rna_commander/tools/rna_tools/rna_commender commit 7ad344d108076116e702e1c1e91cea73d8fcadc4
author rnateam
date Thu, 28 Jul 2016 05:55:25 -0400
parents 8918de535391
children
line wrap: on
line diff
--- a/rbpfeatures.py	Tue Jul 19 10:22:02 2016 -0400
+++ b/rbpfeatures.py	Thu Jul 28 05:55:25 2016 -0400
@@ -1,6 +1,7 @@
 """Compute the RBP features."""
 
 import re
+import sys
 import subprocess as sp
 import uuid
 from os import mkdir
@@ -57,6 +58,9 @@
 
         fasta = fasta_utils.import_fasta(self.fasta)
 
+        if len(fasta) != 1:
+            sys.exit("""Fasta file must contain exactly one sequence.""")
+
         for rbp in sorted(fasta.keys()):
             seq = fasta[rbp]
             text = pfam_utils.sequence_search(rbp, seq)