20191230对AS类型进行定义
Last updated
Last updated
前面通过Blast的方法对可变剪切在不同亚基因组间的保守性进行了分类，然后对每一类基因定义了一个概念，这个概念是紧紧围绕多倍化过程中，可变剪切在二倍中同源基因和四倍体同源基因的变化情况
conservation AS pattern
筛选条件
匹配长度大于200bp
相似度大于90%
覆盖度超过85%
代码
## 在所有基因组中都保守的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文件
## 在所有基因组中都保守的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.bedAS loose
symmetry loose
At asymmetry loose
Dt asymmetry loose
代码
## 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文件
## 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.bedAS gain
symmetry gain
At asymmetry gain
Dt asymmetry gain
代码
:warning:单个基因组存在的事件得再converseBed目录下运行
## 只有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文件
## 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.bedconvergence effect
convergence to gain
convergence to loose
代码
## 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文件
## 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.bedAS事件在不同亚基因组间保守
## 只在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文件
## 只在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文件，然后剔除掉这些就是只在某个基因组中单独存在的剪切事件了。
## 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中都存在的保守事件
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类，对每一类都统计一下对应的事件的数目
##第一大类
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 -lDt中的可变剪切
1_1 1_4 2_2 3_2 3_3 4_1 4_2 5_1At中的可变剪切
1_1 1_3 2_1 3_1 3_2 4_1 4_2 5_2D5中的可变剪切
1_1 1_4 2_1 2_2 2_3 4_2 4_4 5_2A2中的可变剪切
1_1 1_3 2_1 2_2 2_3 4_1 4_3 5_1事件
基因数目
事件数目
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
事件
基因数目
事件数目
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
事件
基因数目
事件数目
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
事件
基因数目
事件数目
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
