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

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容