根据nanopolish的结果统计不同概率下甲基化的数目

1.nanopolish能call出5mc的结果文件。如下图所示：

2.现在统计methylated_frequency在0到1之间每隔0.1分一个阶段 统计个数。脚本如下：

#!/usr/bin/env python 
# -*- coding: UTF-8 -*-

import re,sys,os

num1,num2,num3,num4,num5,num6,num7,num8,num9,num10=0,0,0,0,0,0,0,0,0,0
fa=open("nanopolish结果文件路径"+sys.argv[1]+"/merge.tsv", 'r')    ##传入的参数是样本编号
fb=open("输出文件路径"+sys.argv[1]+"/jieguo.txt", 'a')
for line in fa.readlines():
	str1=line.split("\t")[-2].strip()
	#print(str1)
	if 0<float(str1)<=0.1:
		num1+=1
	elif 0.1<float(str1)<=0.2:
		num2+=1
	elif 0.2<float(str1)<=0.3:
		num3+=1
	elif 0.3<float(str1)<=0.4:
		num4+=1
	elif 0.4<float(str1)<=0.5:
		num5+=1
	elif 0.5<float(str1)<=0.6:
		num6+=1
	elif 0.6<float(str1)<=0.7:
		num7+=1
	elif 0.7<float(str1)<=0.8:
		num8+=1
	elif 0.8<float(str1)<=0.9:
		num9+=1
	elif 0.9<float(str1)<=1:
		num10+=1
result=sys.argv[1]+'\t'+str(num1)+'\t'+str(num2)+'\t'+str(num3)+'\t'+str(num4)+'\t'+str(num5)+'\t'+str(num6)+'\t'+str(num7)+'\t'+str(num8)+'\t'+str(num9)+'\t'+str(num10)+'\n'
fb.write(result)
#print(num1)
#print(int(num2))
fa.close()
fb.close()