On Wed, 2010-07-21 at 10:47 +0200, Pierre-Yves wrote: > s = Seq(seq.seq[start:stop], generic_dna) > seq_out = SeqRecord(s, id = row[col_name]) I found the solution: s = Seq(str(seq.seq[start:stop]), generic_dna) Sorry for the noise, Pierre