鉴定DRMs区域

计算某个区域平均甲基化水平

• 每个已经甲基化的C的甲基化水平： 支持的reads/(支持read+不支持read)

• 所有甲基化的C/比上该区域C的数目

## 计算每个区域中已经甲基化的C的甲基化程度
~/software/bedtools2-2.29.0/bin/intersectBed  -loj -a /data/cotton/zhenpingliu/DNA_methlation_rawData/TM-1/CpG_context_TM1_qs_rep1_count.txt_binom.txt_fdr.bed_end_sorted -b ./IR.txt  >test.bed
## 计算该区域C的个数。求平均值