<div dir="ltr">My ultimate scope is to generate a XML file with two sequences: one protein sequence and one DNA sequence, like that:<div><br></div><div>…</div><div>&lt;AAseq&gt;XLLLAINGVTECFTFAAMS&lt;/AAseq&gt;</div><div>&lt;DNAseq&gt;NTTCTCCTGCTTGCCATCAATGGAGTGACAGAGTGTTTCACATTTG&lt;/DNAseq&gt;<br clear="all"><div>…</div><div><br></div><div>But I am afraid that Bio.SeqIO can only keep one Seq object at a time per record.</div><div><br></div><div>Hoping that I am wrong, can some tell me how to do it?</div><div><br></div><div>a0 is my target record:</div><div><br></div><div>SeqRecord(seq=Seq(&#39;XLLLAINGVTECFTFAAMS&#39;, ProteinAlphabet()), id=&#39;sp|O07835|IORA_THEKO&#39;, name=&#39;sp|O07835|IORA_THEKO&#39;, description=&#39;sp|O07835|IORA_THEKO Indolepyruvate oxidoreductase subunit IorA OS=Thermococcus kodakaraensis (strain ATCC BAA-918 / JCM 12380 / KOD1) GN=iorA PE=1 SV=2&#39;, dbxrefs=[])<br></div><div><br></div><div>and b0 is my source record, from which I only want the DNA sequence,  b0.seq = Seq(&#39;NTTCTCCTGCTTGCCATCAATGGAGTGAC...TTG&#39;, DNAAlphabet())</div><div><br></div><div>Any help would be very much appreciated.</div><div><br></div><div>Many thanks in advance,</div><div><br></div><div>Alan</div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Alan Wilter SOUSA da SILVA, DSc<div>Bioinformatician, UniProt</div><div>European Bioinformatics Institute (EMBL-EBI)</div><div>European Molecular Biology Laboratory</div><div>Wellcome Trust Genome Campus</div><div>Hinxton</div><div>Cambridge CB10 1SD</div><div>United Kingdom</div><div>Tel: +44 (0)1223 494588</div></div></div>
</div></div>