gff文件缺少gene特征

1
awk -F"[\t=]"  'BEGIN{OFS="\t"}{if($3=="mRNA"){print $1,$2,"gene",$4,$5,$6,$7,$8,$9"="$10 ; print $1,$2,$3,$4,$5,$6,$7,$8,$9"="$10"Parent="$10;}else{print $0}}' all.gff3 > all.gff3.new

复制粘贴使用：

1
awk -F"[\t=]"  'BEGIN{OFS="\t"}{if($3=="mRNA"){print $1,$2,"gene",$4,$5,$6,$7,$8,$9"="$10 ; print $1,$2,$3,$4,$5,$6,$7,$8,$9"="$10"Parent="$10;}else{print $0}}'
生物信息学
#生物信息学
gff文件缺少gene特征
https://lixiang117423.github.io/article/gff3gene/
作者
小蓝哥
发布于
2022年10月2日
许可协议