泛基因家族-结构分析

老师,我进行了这个操作,结果如图所示,是不是不对啊,哪里出现问题了呢?

#提取结构变异基因列表

grep "${svgene}" Fam_gene.list|sed "s/${svgene}\t//g"> ${svgene}_gene.list

attachments-2025-10-QBbJro9S68fe0180e0dd3.png

请先 登录 后评论

1 个回答

rzx

默认情况下 grep 匹配的是包含子串 "WRKY14" 的内容,所以会同时匹配上WRKY140~WRKY149。
可以修改命令
grep "\b${svgene}\b" Fam_gene.list|sed "s/${svgene}\t//g"> ${svgene}_gene.list
(添加单词边界 \b 限制匹配范围)

请先 登录 后评论
  • 1 关注
  • 0 收藏,26 浏览
  • 生信小白 提出于 1天前

相似问题