[Biopython] writing cut sites in genbank via biopython

Jocelyne <[email protected]> Wed, 7 Oct 2020 11:48:29 -0700
Newsgroups gmane.comp.python.bio.general
Message-ID <CAKph291PZihfMgGmJ3zb8-39vYqRxrn36D5FY1o1SAwGLqST0w@mail.gmail.com>
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

_______________________________________________
Biopython mailing list  -  [email protected]
https://mailman.open-bio.org/mailman/listinfo/biopython