[Biopython] 2 seqs in one Bio.SeqIO, is it possible?

[Peter Cock]

On Tue, Mar 17, 2015 at 3:44 PM, Alan <alanwilter at gmail.com> wrote:
> My ultimate scope is to generate a XML file with two sequences: one protein
> sequence and one DNA sequence, like that:
>
> …
> <AAseq>XLLLAINGVTECFTFAAMS</AAseq>
> <DNAseq>NTTCTCCTGCTTGCCATCAATGGAGTGACAGAGTGTTTCACATTTG</DNAseq>
> …
>
> But I am afraid that Bio.SeqIO can only keep one Seq object at a time per
> record.

Correct. You could treat one as the primary sequence, and treat
the second sequence as an annotation of the SeqRecord object?

I would also suggest looking at the Codon-based sequence
objects under Bio.codonalign which might be relevant?

Peter