Mercurial > repos > abims-sbr > cds_search
comparison scripts/dico.py @ 7:35e39b4128ba draft
planemo upload for repository https://github.com/abims-sbr/adaptsearch commit b7a3030ea134b5dfad89b1a869db659d72d1145c
| author | abims-sbr |
|---|---|
| date | Wed, 28 Feb 2018 10:38:40 -0500 |
| parents | ff98ed7849fa |
| children |
comparison
equal
deleted
inserted
replaced
| 6:c03ba6a7c5f2 | 7:35e39b4128ba |
|---|---|
| 1 import string | 1 import string, itertools |
| 2 | 2 |
| 3 def dico(F1): | 3 def dico(F1): |
| 4 dicoco = {} | 4 dicoco = {} |
| 5 while 1: | 5 with open(F1, "r") as file: |
| 6 next2 = F1.readline() | 6 for name, query in itertools.izip_longest(*[file]*2): |
| 7 if not next2: | 7 if name[0] == ">": |
| 8 break | 8 fasta_name_query = name[:-1] |
| 9 if next2[0] == ">": | 9 Sn = string.split(fasta_name_query, "||") |
| 10 fasta_name_query = next2[:-1] | 10 fasta_name_query = Sn[0] |
| 11 Sn = string.split(fasta_name_query, "||") | 11 fasta_seq_query = query[:-1] |
| 12 fasta_name_query = Sn[0] | 12 dicoco[fasta_name_query] = fasta_seq_query |
| 13 next3 = F1.readline() | |
| 14 fasta_seq_query = next3[:-1] | |
| 15 dicoco[fasta_name_query]=fasta_seq_query | |
| 16 return(dicoco) | 13 return(dicoco) |
