提取区间内的基因

1、创建区间文件   regions.bed

2、bedtools 提取区间内信息

bedtools intersect -a Triticum_aestivum.IWGSC.59.gtf -b regions.bed > genes_in_regions.gff

3、提取基因ID

awk '$3 == "gene" { print $0 }' genes_in_regions.gff | awk -F '\t' '{print $1,$4,$5,$9}' | awk -F ';' '{print $1}' > select_id

id 带有gene_id和引号，需要替换