“PAV文件中不同的fam号对应了同一个基因后”,这是什么意思,我没太理解,请提供更详细的说明
您好,当我发现我的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