> 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/20191230-dui-as-lei-xing-jin-hang-ding-yi.md).

# 20191230对AS类型进行定义

前面通过Blast的方法对可变剪切在不同亚基因组间的保守性进行了分类，然后对每一类基因定义了一个概念，这个概念是紧紧围绕多倍化过程中，可变剪切在二倍中同源基因和四倍体同源基因的变化情况

## 第一类比较保守的AS

***conservation AS pattern***

![conservation AS pattern](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164457.png)

筛选条件

* 匹配长度大于200bp
* 相似度大于90%
* 覆盖度超过85%

代码

```bash
## 在所有基因组中都保守的IR对应的bed
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]>=1&&a[i][2]>=1&&a[i][1]>=1)print i}}' >converse/1_1
## 都不发生IR事件的基因
cat 1_1.bed  2_2.bed  4_2.bed 1_3.bed  2_3.bed  5_1.bed 1_4.bed  3_2.bed  5_2.bed 2_1.bed  4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "Ghir_D" |awk -F "_" '$1~/Ghir/{print $1"_"$2}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_Dt
cat 1_1.bed  2_2.bed  4_2.bed 1_3.bed  2_3.bed  5_1.bed 1_4.bed  3_2.bed  5_2.bed 2_1.bed  4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "Ghir_A" |awk -F "_" '$1~/Ghir/{print $1"_"$2}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_At
cat 1_1.bed  2_2.bed  4_2.bed 1_3.bed  2_3.bed  5_1.bed 1_4.bed  3_2.bed  5_2.bed 2_1.bed  4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "Gor" |awk -F "_" '{print $1}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_D5
cat 1_1.bed  2_2.bed  4_2.bed 1_3.bed  2_3.bed  5_1.bed 1_4.bed  3_2.bed  5_2.bed 2_1.bed  4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "evm" |awk -F "_" '{print $1}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_A2
cat tmp* |sort |uniq |cat - ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold |sort |uniq -u >../converseBed/1_2
## 只在A2与At中发生
awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][0]>=1&&a[i][2]==0)print i}}' >converse/1_3
## 只在D5与Dt中发生
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]==0&&a[i][1]>=1)print i}}' >converse/1_4
```

做成bed文件

```bash
## 在所有基因组中都保守的IR对应的bed
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]>=1&&a[i][2]>=1&&a[i][1]>=1)print i}}' |xargs  -I {} grep {} Dt_query_blast.txt| awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/1_1.bed

## 只在A2与At中发生
awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][0]>=1&&a[i][2]==0)print i}}'|xargs  -I {} grep {} A2_query_blast.txt |awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/1_3.bed

## 只在D5与Dt中发生
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]==0&&a[i][1]>=1)print i}}' |xargs  -I {} grep {} Dt_query_blast.txt| awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/1_4.bed
```

## 第二类多倍化后四倍体中AS发生丢失

***AS loose***

* ***symmetry loose***
* ***At asymmetry loose***
* ***Dt asymmetry loose***

![AS loose](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164455.png)

代码

```bash
## Dt中AS 丢失
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][0]>=1&&a[i][1]>=1)print i}}' >converse/2_1
## At中AS 丢失
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]>=1)print i}}'  >converse/2_2
## At与Dt都丢失
awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][1]>=1&&a[i][2]==0)print i}}' >converse/2_3
```

做成bed文件

```bash
## Dt中AS 丢失
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][0]>=1&&a[i][1]>=1)print i}}' |xargs  -I {} grep {} At_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/2_1.bed

## At中AS 丢失
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]>=1)print i}}'  |xargs  -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/2_2.bed

## At与Dt都丢失
awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][1]>=1&&a[i][2]==0)print i}}'|xargs  -I {} grep {} A2_query_blast.txt | awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq   >converseBed/2_3.bed
```

## 第三类多倍化后四倍体AS重新获得

***AS gain***

* ***symmetry gain***
* ***At asymmetry gain***
* ***Dt asymmetry gain***

![AS gain](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164452.png)

代码

:warning:单个基因组存在的事件得再converseBed目录下运行

```bash
## 只有At 获得新的AS事件
cat 1_1.bed 1_3.bed  2_1.bed 3_2.bed 4_1.bed 4_2.bed 5_2.bed |awk -F "\t" '$5~/Ghir_A/{print ">"$5}'|cat - ../At_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/3_1
## At与Dt 获得同样的AS事件
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][2]==0&&a[i][0]>=1)print i}}'  >converse/3_2
## 只有Dt 获得新的AS事件
cat 1_1.bed 1_4.bed 2_2.bed 3_2.bed 4_1.bed 4_2.bed 5_1.bed|awk -F "\t" '$5~/Ghir_D/{print ">"$5}' |cat - ../Dt_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/3_3
```

做成bed文件

```bash
## At与Dt 获得同样的AS事件
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][2]==0&&a[i][0]>=1)print i}}'|xargs  -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq   >converseBed/3_2.bed
```

## 第四类四倍体中两个亚基因组有趋同效应

***convergence effect***

* ***convergence to gain***
* ***convergence to loose***

![convergence effect](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164459.png)

代码

```bash
## D5中不存在
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]>=1&&a[i][0]>=1&&a[i][1]==0)print i}}' >converse/4_1
## A2中不存在
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]>=1&&a[i][2]==0&&a[i][0]>=1)print i}}' >converse/4_2
## 只有A2中有
cat 1_1.bed 1_3.bed 2_1.bed 2_2.bed 2_3.bed 4_1.bed 5_1.bed|awk -F "\t" '$5~/evm/{print ">"$5}'|cat - ../A2_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g'  >../converse/4_3
## 只有D5中有
cat 1_1.bed 1_4.bed 2_1.bed 2_2.bed 2_3.bed 4_2.bed 5_2.bed|awk -F "\t" '$5~/Gor/{print ">"$5}'|cat - ../D5_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/4_4
```

做成bed文件

```bash
## D5中不存在
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]>=1&&a[i][0]>=1&&a[i][1]==0)print i}}' |xargs  -I {} grep {} At_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/4_1.bed

## A2中不存在
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]>=1&&a[i][2]==0&&a[i][0]>=1)print i}}'|xargs  -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq   >converseBed/4_2.bed
```

## 第五类

AS事件在不同亚基因组间保守

```bash
## 只在A2与Dt中保守
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]==0)print i}}' >converse/5_1
## 只在D5与At中保守
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][1]>=1&&a[i][0]==0)print i}}' >converse/5_2
```

做成bed文件

```bash
## 只在A2与Dt中保守
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]==0)print i}}'|xargs  -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/5_1.bed
## 只在D5与At中保守
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][1]>=1&&a[i][0]==0)print i}}'|xargs  -I {} grep {} At_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq  >converseBed/5_2.bed
```

## 求只在单个基因组中的剪切事件

在只某一个基因组中保守的类型，先统计与其他基因组中保守的类型的bed1文件，然后剔除掉这些就是只在某个基因组中单独存在的剪切事件了。

```bash
## At和其他基因组中共有的类型
1_1 1_3 2_1 3_2 4_1 4_2 5_2 
cat 1_1.bed 1_3.bed  2_1.bed 3_2.bed 4_1.bed 4_2.bed 5_2.bed |awk -F "\t" '$5~/Ghir_A/{print ">"$5}'|cat - ../At_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $3"_"$4,$5,$6,$8,$0}' OFS="\t" >At.bed 
## Dt和其他基因组中共有的类型
cat 1_1.bed 1_4.bed 2_2.bed 3_2.bed 4_1.bed 4_2.bed 5_1.bed|awk -F "\t" '$5~/Ghir_D/{print ">"$5}' |cat - ../Dt_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $3"_"$4,$5,$6,$8,$0}' OFS="\t" >Dt.bed
## A2和其他基因组中共有的类型
cat 1_1.bed 1_3.bed 2_1.bed 2_2.bed 2_3.bed 4_1.bed 5_1.bed|awk -F "\t" '$5~/evm/{print ">"$5}'|cat - ../A2_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $2,$3,$4,$6,$0}' OFS="\t" >A2.bed
## D5和其他基因组共有的类型
cat 1_1.bed 1_4.bed 2_1.bed 2_2.bed 2_3.bed 4_2.bed 5_2.bed|awk -F "\t" '$5~/Gor/{print ">"$5}'|cat - ../D5_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $2,$3,$4,$6,$0}' OFS="\t" >D5.bed
```

## A2与At 、D5与Dt间保守的剪切事件

```bash
## A2与At中都存在的保守事件
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]>=1)print i}}' >converse/A2_At_converse.txt
## D5与Dt中都存在的保守事件
awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]>=1)print i}}' >converse/D5_Dt_converse.txt
```

## 统计基因数与事件数

总共有16类，对每一类都统计一下对应的事件的数目

```bash
##第一大类
wc -l 1_1 
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 1_1 |sort |uniq |wc -l
wc -l 1_2
wc -l 1_3 
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 1_3 |sort |uniq |wc -l
wc -l 1_4 
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 1_4 |sort |uniq |wc -l
##第二大类
wc -l 2_1
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 2_1 |sort |uniq |wc -l
wc -l 2_2
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 2_2 |sort |uniq |wc -l
wc -l 2_3
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 2_3|sort |uniq |wc -l
##第三大类
wc -l ../converseBed/At.bed
cut -f5 ../converseBed/At.bed|awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}'  |sort |uniq |wc -l
wc -l 3_2
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 3_2 |sort |uniq |wc -l
wc -l ../converseBed/Dt.bed
cut -f5  ../converseBed/Dt.bed |awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}'|sort |uniq |wc -l
## 第四大类
wc -l 4_1
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 4_1 |sort |uniq |wc -l
wc -l 4_2
awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 4_2 |sort |uniq |wc -l
wc -l ../converseBed/A2.bed
cut -f5 ../converseBed/A2.bed |awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' |sort |uniq |wc -l
wc -l ../converseBed/D5.bed
cut -f5 ../converseBed/D5.bed |awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' |sort |uniq |wc -l
```

## intronR

| 事件   | 基因数目 | 事件数目  |
| ---- | ---- | ----- |
| 1\_1 | 840  | 1489  |
| 1\_2 | 8071 |       |
| 1\_3 | 1122 | 1869  |
| 1\_4 | 1438 | 2247  |
| 2\_1 | 651  | 977   |
| 2\_2 | 600  | 861   |
| 2\_3 | 1043 | 1541  |
| At   | 3858 | 7683  |
| 3\_2 | 1038 | 1571  |
| Dt   | 3924 | 8053  |
| 4\_1 | 579  | 876   |
| 4\_2 | 825  | 1245  |
| A2   | 4087 | 8451  |
| D5   | 5066 | 10217 |

## ExonS

| 事件   | 基因数目  | 事件数目 |
| ---- | ----- | ---- |
| 1\_1 | 38    | 42   |
| 1\_2 | 15390 |      |
| 1\_3 | 103   | 108  |
| 1\_4 | 235   | 254  |
| 2\_1 | 34    | 35   |
| 2\_2 | 38    | 39   |
| 2\_3 | 83    | 87   |
| At   | 712   | 800  |
| 3\_2 | 66    | 68   |
| Dt   | 661   | 740  |
| 4\_1 | 22    | 23   |
| 4\_2 | 86    | 88   |
| A2   | 656   | 759  |
| D5   | 1406  | 1648 |

## AltA

| 事件   | 基因数目  | 事件数目 |
| ---- | ----- | ---- |
| 1\_1 | 137   | 167  |
| 1\_2 | 12517 |      |
| 1\_3 | 212   | 262  |
| 1\_4 | 423   | 477  |
| 2\_1 | 90    | 110  |
| 2\_2 | 62    | 69   |
| 2\_3 | 172   | 198  |
| At   | 1805  | 2312 |
| 3\_2 | 367   | 421  |
| Dt   | 1777  | 2359 |
| 4\_1 | 81    | 91   |
| 4\_2 | 218   | 261  |
| A2   | 1555  | 2126 |
| D5   | 2127  | 2809 |

## AltD

| 事件   | 基因数目  | 事件数目 |
| ---- | ----- | ---- |
| 1\_1 | 82    | 95   |
| 1\_2 | 13275 |      |
| 1\_3 | 169   | 216  |
| 1\_4 | 284   | 331  |
| 2\_1 | 66    | 85   |
| 2\_2 | 44    | 50   |
| 2\_3 | 125   | 146  |
| At   | 1671  | 2122 |
| 3\_2 | 316   | 364  |
| Dt   | 1714  | 2243 |
| 4\_1 | 56    | 62   |
| 4\_2 | 148   | 181  |
| A2   | 1291  | 1688 |
| D5   | 1550  | 1934 |

## 对同源基因间的剪切取交集

Dt中的可变剪切

```bash
1_1 1_4 2_2 3_2 3_3 4_1 4_2 5_1
```

At中的可变剪切

```bash
1_1 1_3 2_1 3_1 3_2 4_1 4_2 5_2
```

D5中的可变剪切

```bash
1_1 1_4 2_1 2_2 2_3 4_2 4_4 5_2
```

A2中的可变剪切

```bash
1_1 1_3 2_1 2_2 2_3 4_1 4_3 5_1
```