braker2的错误

错误一:
提示
ImportError: Bio.Alphabet has been removed from Biopython. In many cases, the alphabet can simply be ignored and removed from scripts. In a few cases, you may need to specify the molecule_type as an annotation on a SeqRecord for your script to work correctly. Please see https://biopython.org/wiki/Alphabet for more information.

经研究,原来是biopython 1.78版本删除了Alphabet 模块

解决方法:
删除1.78,安装1.76

pip3 unintall biopython
pip3 install biopython==1.76

错误二:
提示
Training gene file in genbank format /raid1/wzw/Ensete/braker/braker/train.f.gb does not contain any training genes. Possible known cause: no training genes with sufficient extrinsic evidence support or of sufficient length were produced by GeneMark-ES/ET. If you think this is the cause for your problems, consider running BRAKER with different evidence or without any evidence (--esmode) for training.

genemark-ET的基因预测集是空的,两个可能,一个是没把证书复制,另一个是braker2里面的有一个组件有问题。我是第二个情况。
寻找到filterGenesIn_mRNAname.pl这个文件

$_ =~ m/transcript_id \"(.*)\"/

改为

$_ =~ m/transcript_id \"(.*?)\"/
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

友情链接更多精彩内容