基于前面已经分好的类进行统计

由于A基因组与D基因组间的差异

统计一下A基因组中内含子保留的长度与D基因组中内含子保留的长度的比较

## Dt与D5间保守
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}}' |awk -F "_" '{print "Dt_D5",$(NF-1)}' OFS="\t" >intronR_length/genome_length.txt
## At与A2间保守
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}}'|awk -F "_" '{print "At_A2",$(NF-1)}' OFS="\t" >>intronR_length/genome_length.txt