Fam_genelist中基因个数超出了所有的鉴定的数据

您好,当我发现我的PAV文件中不同的fam号对应了同一个基因后,我检查我在32个材料中共鉴定了12359个基因,但生成的Fam_gene.list文件中基因总数为12436

我尝试了重新跑perl $scriptsdir/PAV_list.pl -pangenel apple.WRKYgene.pan-genes.list -spid new_id2.txt -unigenel apple.WRKYgene.uniq.list -fam WRKY

awk -F"\t" '{if($0!~/^Fam_N/){for(i=2;i<=27;i++){print $1"\t"$i}}}' PAV_gene.list |sed "s/;$//"|awk -F"\t" '{gsub(";","\n"$1"\t",$0);print}'|awk -F"\t" '{if($2!=""){print}}' > Fam_gene.list

后还是如此。请老师帮忙解答一下,谢谢

attachments-2024-12-HCiQRNlp675a94f86c13d.pngattachments-2024-12-HCiQRNlp675a94f86c13d.png

请先 登录 后评论

1 个回答

每天学习一点点

PAV文件中不同的fam号对应了同一个基因后”,这是什么意思,我没太理解,请提供更详细的说明

请先 登录 后评论
  • 1 关注
  • 0 收藏,133 浏览
  • 树心台 提出于 2024-12-12 15:48

相似问题