skbio - 有没有办法对 skbio 序列进行切片，以便将间隔正确传输到新序列？

Question

我有一个带有注释的 genbank 文件。

我想切片通过读取该文件生成的 seq obj 并将生成的子序列写入新的 genbank 文件。但是要使此功能有用，功能需要随之而来。我想查看该子序列中的引物结合位点、启动子、基因等。

但是切片似乎在设计上放弃了 interval_metadata 属性。为什么是这样？有没有什么办法可以做一个看起来非常标准且非常有用的东西，它有大量的用例？

标签： skbio