Linux练习-gff

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