By Gengxin, 31 March, 2026
Forums

存在问题:计算红豆杉中基因对之间的Ka.Ks的值全都是NaN

原因:每条蛋白质序列末尾都有X

解决:去掉蛋白序列末尾的 X

错误:...PILLX  正确:...PILL

sed -i '/^>/! s/X$//' Twallichiana.final.pep
# /^>/!:只对非标题行(序列行)操作
# s/X$//:替换行尾 X 为空