[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