1. 查看gff文件Y染色体基因数目
ensemble下载y然染色体gff注释文件
zcat Homo_sapiens.GRCh38.103.chromosome.Y.gff3.gz | grep -v "#" | cut -f 3 | grep "^gene" | sort |uniq -c
2.第三列类几种
zcat Homo_sapiens.GRCh38.103.chromosome.Y.gff3.gz | grep -v "#" | cut -f 3 | sort | uniq -c | sort -n -k 1
3.匹配exon行再反向输出
zcat Homo_sapiens.GRCh38.103.chromosome.Y.gff3.gz | grep -v "exon" | less -SN
4.匹配UTR或者CDS的行
zcat Homo_sapiens.GRCh38.103.chromosome.Y.gff3.gz | grep -e "UTR" -e "CDS" | head