基于liftoff进行基因组注释操作记录

步骤

安装


conda create -n liftoff

conda install liftoff

执行代码


liftoff -g (参考基因组gff3注释文件,如:Gmax_275_Wm82.a2.v1.gene.gff3) -o (输出文件名).gene.gff3 -p 28 -chroms (染色体号对应情况) -u ./unmap.txt(未匹配的日志) ../xxx.fa(需要匹配的基因组) ../Gmax_275_v2.0.fa(参考基因组)

代码参数

-g 参考基因组的注释文件

-o 想要注释的基因组输出的注释文件

-p 线程数

-u 未匹配到的基因日志文件

命令最后两个 target 需要匹配的基因组 reference 参考基因组

-chroms 为了提高基因映射的准确性,这里仅允许相同染色体间的基因进行映射,因此需要指定两个基因组间的染色体对应关系,即-chroms参数,其格式如下:


chrA01,chrA01

chrA02,chrA02

chrA03,chrA03

chrA04,chrA04

chrA05,chrA05

...

Uncluster,Uncluster

Ref

https://github.com/agshumate/Liftoff#readme

https://www.jianshu.com/p/c18f7cac43c6

欢迎大家交流。祝大家代码不报错,科研步步过~
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。