[Biopython] AC lines with multiple accession numbers in Swissprot files
Adam Sjøgren
asjo at koldfront.dk
Wed May 29 09:07:03 UTC 2019
Hi,
When using BioPython to parse a record such as this one:
· https://www.uniprot.org/uniprot/Q12558.txt
where the AC line looks like this:
$ grep ^AC /tmp/Q12558.txt
AC Q12558; Q7LWA9;
the SeqRecord.id attribute contains 'Q12558':
>>> import Bio.SeqIO
>>> [r] = list(Bio.SeqIO.parse("/tmp/Q12558.txt", format="swiss"))
>>> r.id
'Q12558'
How can I access the other value, 'Q7LWA9'?
Best regards,
Adam
--
"Everything needs to change. Adam Sjøgren
And it has to start today." asjo at koldfront.dk
More information about the Biopython
mailing list