[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