错误一:
提示
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 \"(.*?)\"/
网友评论