对使用misa得到叶绿体基因组SSR的结果进行注释
关于叶绿体基因组的文章中,SSR的鉴定是一件常做的分析,文章中往往会列出每个ssr所在的基因组的位置信息,如:在某基因处,在内含子区域,或者在基因间区,那么这些信息是怎么得到的呢?
诸如此类的信息:
SS1 | p2 | (TA)6 | 12 | 9662 | 9673 | IGS(trnR-UCU,atpA) |
SS1 | p1 | (A)11 | 11 | 12336 | 12346 | atpF(intron) |
SS1 | p2 | (TA)5 | 10 | 19614 | 19623 | rpoC2 |
手填的话是可以实现的,但是比较麻烦,如果数量多的话是个很大的工作量。所以可以使用脚本代替我们做这件事情。
实现原理很简单,只需要记录genbank中基因的位置,然后把该信息映射到每个ssr所在的位置就可以了。
需要脚本的小伙伴可加群获取哦~ 群号:636121790