Sebastian Bassi (aka Bloguero Connor)

Prueba, no molestar!

Probando meter codigo, si funciona pondré algo coherente (o al menos lo intentaré).

from __future__ import with_statement
from Bio import SeqIO, SeqRecord, Seq

#para planillas html de Laura
ain='/mnt/hda2/bio/cln.txt'

cln=set()
with open(ain) as f:
    for line in f:
        cln.add(line[:-1])
        #print line,

#sco esas del fasta y las pongo en otro fasta.
#recorro el fasta
i = 0
fasta69 = []
for record in SeqIO.parse(open('/mnt/hda2/bio/bacs.v247.seq'), "fasta") :
    if record.id in cln:
        fasta69.append(record)

with open('/mnt/hda2/bio/fasta69.txt','w') as outf:
    SeqIO.write(fasta69,outf,"fasta")

    

Etiquetas: python, test

publicadas por Sebastian @ Link to this post8:32 p. m.