parse_primers causes some initialization problem

[manulera]

Below code to reproduce

from pydna.primer import Primer
from pydna.parsers import parse_primers
from Bio.SeqFeature import SimpleLocation
from pydna.dseqrecord import Dseqrecord

primer = Primer('AAAAA')
dseqrecord_primer = Primer(Dseqrecord('AAAAA', id='primer'))
parsed_primer, *_ = parse_primers('>primer\nAAAAA')

loc = SimpleLocation(1, 3)

loc.extract(primer).seq.rc()  # Does NOT give error
loc.extract(dseqrecord_primer).seq.rc()  # Does NOT give error
loc.extract(parsed_primer).seq.rc()  # Gives error

Error:

    loc.extract(parsed_primer).seq.rc()  # Gives error
    ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
AttributeError: 'Seq' object has no attribute 'rc'

[BjornFJohansson]

Probably because a biopython Seq is returned. They have only a ._reverse_complement() method, no .rc().

type(loc.extract(parsed_primer).seq)
Out[24]: Bio.Seq.Seq 

[manulera]

Yes, but I was not able to pinpoint where this happens. I did not look too hard.