
import os, string 

printsdb = open('/users/spyros/folder1/python/PRINTSmotifs/prints41_1.kdat', 'r')
lines = printsdb.readlines()

# find PRINTS name entries
for line in lines:
    while line.startswith('gc;'):
        protname = line.lstrip('gc;')
        protnames = (protname)
        if not protnames:
            print('error in creating tuple') # check if tuple is true or false 
        #print(protnames)
        break
    
query = input("search a protein: ")
query = query.upper() 
if query in protnames: 
    print("\nDisplaying Motifs")
else:
    print("\nentry not in database")
    
# Parse motifs 
def extract_motifs(query):
    motif_id = ()
    motif = ()
    while query in lines:  ####for query, get motif_ids and motifs 
        while line.startswith('ft;'):
            motif_id = line.lstrip('ft;')
            motif_ids = (motif_id)
            #print(motif_id) 
            while line.startswith('fd;'):
                motif = line.lstrip('fd;')
                motifs = (motif)
            #print(motif)
            return motif_id, motif 
        
if __name__ == '__main__':
    final_motifs = extract_motifs('query') 

    

    



    



















    



    







    

    

    






    
