Bismark得到的结果如何导入DSS

# -*- coding:utf-8 -*-

import sys

def run(dirs_path):

    f = open(dirs_path ,'r')

    name = dirs_path.split('/')[-1].split('.')[0]

    w = open('/Users/apple/ERV/'+name+'.txt','w')

    w.write('chr'+'\t'+'pos'+'\t'+'N'+'\t'+'X'+'\n')

    for line in f:

        d = line.strip().split('\t')

        col_3 = int(d[-2])+int(d[-1])

        w.write(d[0]+'\t'+d[1]+'\t'+str(col_3)+'\t'+d[-1]+'\n')

    f.close()

    w.close()

if__name__ =='__main__':

    dirs_path = sys.argv[1]

    run(dirs_path)

这里面唯一需要修改的就是 w = open('/Users/apple/ERV/'+name+'.txt','w')这里，你输入自己的路径就行，其他不动。

用的时候你先把这段代码复制成自己的脚本，后缀.py，比如convert.py

然后用的时候就，python convert.py test.bismark.cov（就是你需要转换的文件）

代码改自：https://zhuanlan.zhihu.com/p/109261289