gff

001、数据和脚本 [root@PC1 test2]# ls GCF_001704415.1_ARS1_genomic.gff record.sh 002、脚本 [root@PC1 test2]# cat record.sh ## 脚本内容 #!/bin/bas ## step1: eliminat ......

001、 [b20223040323@admin2 test]$ ls a.gff [b20223040323@admin2 test]$ cat a.gff region 1 pseudogene 2 transcript 3 exon 4 pseudogene 5 transcript 6 ex ......

001、生成基因名称的列表 awk -F "\t" '$3 == "gene" && $NF ~ /gene=/ {print $NF}' chr1.gff | sed 's/\(.*\)\(gene=[^;]\+\)\(.*\)/\2/' | sort | uniq > gene.list 002 ......