001、生成基因名称的列表
awk -F "\t" '$3 == "gene" && $NF ~ /gene=/ {print $NF}' chr1.gff | sed 's/\(.*\)\(gene=[^;]\+\)\(.*\)/\2/' | sort | uniq > gene.list
002、
Linux 中 shell 脚本实现根据gff统计每一个基因的转录本数目
发布时间 2023-05-16 22:37:33作者: 小鲨鱼2018
001、生成基因名称的列表
awk -F "\t" '$3 == "gene" && $NF ~ /gene=/ {print $NF}' chr1.gff | sed 's/\(.*\)\(gene=[^;]\+\)\(.*\)/\2/' | sort | uniq > gene.list
002、