你为什么要替换你的库文件的内容呢,你直接用库文件做注释就行了,脚本里这里也没有让你修改ID:
老师,我正在做泛基因家族分析,结构变异与基因注释部分,#转换gff3到genePred 格式并转换转录本格式至基因ID,尝试了不同的指令,均未能成功修改 sed -i "s/\\.[0-9]*$//" ${ref}_refGene.txt ,之后在excel中完成修改,结果如下:之后,#建立注释库
retrieve_seq_from_fasta.pl --format refGene --seqfile ${ref}.fa ${ref}_refGene.txt --out ${ref}_refGene.fa ,在此过程中出现许多错误并得到最终的${ref}_refGene.fa,请老师帮忙看看是否正常