seqkit安装+genomeTools安装+awk

# 下载并解压(以 Linux 为例)
wget https://github.com/shenwei356/seqkit/releases/download/v2.6.0/seqkit_lin...
tar -zxvf seqkit_linux_amd64.tar.gz

# 将 seqkit 移动到可执行路径(如 /usr/local/bin)
sudo mv seqkit /usr/local/bin/

# 验证安装
seqkit version
 

 

###awk命令替换pep序列中间的*并记录位置到replaced_positions.log

awk '/^>/ {print; next} {gsub(/\*$/, ""); orig=$0; gsub(/\*/, "X"); print; if (orig != $0) print ">" FILENAME ": replaced * at positions: " > "replaced_positions.log"}' input.fasta > processed.fasta