建议用课程里面的,ATAG软件包处理一下,再用这个脚本统计, 不然可能由于GFF文件没有排序,格式不标准导致出差错; #保留蛋白编码基因agat_sp_filter_feature_by_attribute_value.pl --gff $database/homo_protein/osa.gff.gz --attribute biotype --value protein_coding -t '!' -o Oryza_sativa.protein_coding.gff3#...
回答于 2023-08-14 09:25
输入文件都截图看看的,文件之间染色体编号不一致吧,你检查一下;
回答于 2023-08-14 09:22