老师,我进行了这个操作,结果如图所示,是不是不对啊,哪里出现问题了呢?
#提取结构变异基因列表
grep "${svgene}" Fam_gene.list|sed "s/${svgene}\t//g"> ${svgene}_gene.list
默认情况下 grep 匹配的是包含子串 "WRKY14" 的内容,所以会同时匹配上WRKY140~WRKY149。可以修改命令grep "\b${svgene}\b" Fam_gene.list|sed "s/${svgene}\t//g"> ${svgene}_gene.list(添加单词边界 \b 限制匹配范围)