[Biopython] writing cut sites in genbank via biopython

Jocelyne jocelyne at gmail.com
Wed Oct 7 18:48:29 UTC 2020


Hello,
I'm trying to write a feature with location "1^2" (a cut site between the
first and second base).

If I parse a file that contains a feature with location "1^2", I get back a
FeatureLocation(ExactPosition(1), ExactPosition(1))

However, if I write a file with `FeatureLocation(ExactPosition(1),
ExactPosition(1))` , I get "1^1"

I have tried:
* FeatureLocation(1,1)   -> 1^1
* FeatureLocation(1,2)   -> 2
* FeatureLocation(ExactPosition(1), ExactPosition(1))   -> 1^1
* FeatureLocation(ExactPosition(1), ExactPosition(2))   -> 2
* FeatureLocation(FeatureLocation(BetweenPosition(0,0,1),
BetweenPosition(1,1,2)))  -> (1^2)..(1^2)

I'm a bit at a loss at this point. Does anyone know how to "1^2" to a
genbank file?

Any help is appreciated.

Thank you!

Jocelyne
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.open-bio.org/pipermail/biopython/attachments/20201007/a0daa532/attachment.htm>


More information about the Biopython mailing list