2024-03-02从gff文件中获取多区域基因ID

awk -F "\t" 'BEGIN {OFS="\t"}$3=="gene" {print $1,$4,$5,$9}' final_sorted_with_promoter.change.gff3  >final_sorted_with_promoter.change.gene
bedtools intersect -a region.xls -b final_sorted_with_promoter.change.gene -wa -wb | awk -F ";" '{print $1}'> EAR_SR3.gene.name

注意,region.xls和 final_sorted_with_promoter.change.gene文件均为bed格式。
bed格式前3列为必须有的,分别为染色体,起始位置,终止位置。


image.png

由于 region.xls 我在excel表格整理出来的,格式是dos系统的,会报错。因此需要在linux下打开后用set ff=unix改成linux系统的。

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容