步骤
安装
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