热点新闻
基于liftoff进行基因组注释操作记录
2023-07-16 12:15  浏览:8341  搜索引擎搜索“手机低淘网”
温馨提示:信息一旦丢失不一定找得到,请务必收藏信息以备急用!本站所有信息均是注册会员发布如遇到侵权请联系文章中的联系方式或客服删除!
联系我时,请说明是在手机低淘网看到的信息,谢谢。
展会发布 展会网站大全 报名观展合作 软文发布

步骤

安装

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

欢迎大家交流。祝大家代码不报错,科研步步过~
发布人:c274****    IP:61.145.45.***     举报/删稿
展会推荐
让朕来说2句
评论
收藏
点赞
转发