[Biopython] Writting fasta file

xyz mitlox at op.pl
Wed Aug 25 07:22:24 EDT 2010


Hello,
the following code saves a multiple fasta file in a list and after 
reverse complement it writes the content back in a fasta file.

   in_handle = open(opts.inputFasta, "r")
   records = list(SeqIO.parse(in_handle, "fasta"))

   for record in records:
       record.seq = record.seq.reverse_complement()

   out_handle = open(opts.outputFasta, "w")
   SeqIO.write(records, out_handle, "fasta")

How is it possible to do it without a list?

Thank you in advance.


More information about the Biopython mailing list