> For the complete documentation index, see [llms.txt](https://zpliu.gitbook.io/booknote/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://zpliu.gitbook.io/booknote/ke-bian-jian-qie/shu-ju-chu-li-liu-cheng/20200112-gen-ju-ji-yin-ir-bao-shou-cheng-du-dui-ji-yin-jin-hang-go-fu-ji.md).

# 根据IR的保守程度对基因进行GO富集分析

| 类型         | 事件数目 | 基因数  |
| ---------- | ---- | ---- |
| At 与 A2保守  | 6474 | 2387 |
| Dt 与 D5 保守 | 7733 | 2814 |
| At与Dt保守    | 6418 | 2494 |
| A2与D5保守    | 6229 | 2418 |

## 多倍化后只在At上保守

在选择压的作用下At基因保持一定的保守性，而Dt基因组则可以随意的剪切

```bash
awk '$5==0&&$7>0{print $0}' 1
awk '$5==0&&$7>0{print $1"\
n"$3}' 1|xargs  -I {} grep {} ~/genome_data/Ghirsutum_genome_HAU_v1.1/Gh_Noscagenes_GO_V3.annot  >GO/At_conserve.go
```

## 只在Dt上保守的

```bash
awk '$7==0&&$5>0{print $0}' 1
awk '$7==0&&$5>0{print $0}' 1|xargs  -I {} grep {} ~/genome_data/Ghirsutum_genome_HAU_v1.1/Gh_Noscagenes_GO_V3.annot  >GO/Dt_conserve.go
```

## At与Dt中都保守

```bash
awk '$7==1&&$5==1||($7=="NA"&&$5=="NA"){print $0}' 1
awk '$7==1&&$5==1||($7=="NA"&&$5=="NA"){print $1"\n"$3}' 1|xargs  -I {} grep {} ~/genome_data/Ghirsutum_genome_HAU_v1.1/Gh_Noscagenes_GO_V3.annot >GO/At_Dt_conserve.go
```

## At与Dt中都不保守

```bash
awk '$7==0&&$5==0{print $1"\n"$3}' 1
```