[Biopython] Cookbook suggestion
Justin Gibbons
jgibbons1 at mail.usf.edu
Sat Apr 13 16:13:56 EDT 2013
I want to add the following to the cookbook but I am unable to create an
account.
#using SeqIO.write() without holding records in memory.
from Bio import SeqIO
seq_ids=set() #create an empty set to hold the sequence IDs.
indexed_fasta=SeqIO.index(file_path, 'fasta') #Can be searched by sequence
ID but is not held in memory
for seq_record in SeqIO.parse(file_path, 'fasta'):
#Filter according to some critria:
seq_ids.add(seq_record.id)
#write the fasta records to a new file using SeqIO.write()
SeqIO.write([indexed_fasta[seq_id] for seq_id in seq_ids], new_file_path,
'fasta')
So if someone who can edit the cookbook wants to add it feel free to.
Justin Gibbons
More information about the Biopython
mailing list