1，从gtf中提取外显子
awk 'BEGIN{OFS="\t";} $3=="exon" {print $1,$4-1,$5}' Iyun.gtf |bedtools sort |bedtools merge -i - > my_exon.bed
- 提取内含子
awk 'BEGIN{OFS="\t";} $3=="transcript" {print $1,$4-1,$5}' Iyun.gtf |bedtools sort |bedtools subtract -a stdin -b my_exon.bed.gz > my_intron.bed
3.网站
[GitHub - davetang/defining_genomic_regions: Define regions in the genome](https://github.com/davetang/defining_genomic_regions)
