分享一个自己写叶绿体基因组注释脚本的方法
叶绿体基因组注释篇-注释脚本研发
B站链接
花费半个多月时间自己写了一个叶绿体基因组注释脚本,流程如上图。一共1141行
写脚本的过程中,又加深了对叶绿体基因组的认识,虽然第一个版本的写好了,使用的过程中还是会报出一些奇奇怪怪的Bug,最近可能会大修一次。
脚本运行大概只需要3S不到的时间
由于自己对叶绿体基因组比较熟悉,做过近200种植物的组装注释,开发这个脚本的时候把自己见过的情况都考虑在内了,所以准确性还是挺高的,有喜欢写脚本的 可以按照上面的思维导图试一试呀。
叶绿体基因组注释需要注意的问题还是挺多的,即使是自己写的脚本后期也需要进行手动较正,注释的方法可以参考上两篇博文
我自己也比较喜欢码代码,如果大家有叶绿体想分析的内容 欢迎评论区留言