> 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/20191231-tong-ji-fen-xi.md).

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

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

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

```bash
## 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
```

[![intronR长度比较](https://s2.ax1x.com/2019/12/31/llHK1A.md.png)](https://imgchr.com/i/llHK1A)