[Biopython-dev] Error in SeqFeature.CompoundLocation parsing NCBI efetch format

[Brynjar Smári Bjarnason]

Hello.

I see CompoundLocation is quite new. I am currently using anaconda
(Python 2.7.6 :: Anaconda 1.8.0 (64-bit)) and BioPython 1.62.

I am fetching gi values and using SeqIO to parse them. So far most of
them work but I found one that fail.

Code:

p = Entrez.efetch(db="protein", rettype="gp", retmode="text",id="494379")
seq = SeqIO.read(p,"gb")

Gives error:
ValueError: CompoundLocation should have at least 2 parts

With quite long stack trace and the last one being:

/Bio/SeqFeature.pyc:
       996         if len(self.parts) < 2:
--> 997             raise ValueError("CompoundLocation should have at
least 2 parts")

Any suggestions on how to fix this, and maybe what is different with
this gi from the rest of them (one gi that works: 10342)?

Brynjar