split-clusters.py 文件源码

python
阅读 33 收藏 0 点赞 0 评论 0

项目:zika-pipeline 作者: zibraproject 项目源码 文件源码
def main():
        records = SeqIO.to_dict(SeqIO.parse(open(sys.argv[1]), 'fasta'))

    reader = csv.DictReader(sys.stdin, dialect="excel-tab")
    clusters = list(reader)

    groups = set([c['group'] for c in clusters])

    for group in groups:
        print "cluster%s\t%s-cluster%s" % (group, sys.argv[1], group)
        with open('%s-cluster%s' %(sys.argv[1], group), 'w') as fout:
            SeqIO.write([records[i['node']] for i in clusters if i['group'] == group], fout, 'fasta')
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号