Advertisement

统计CDS中的密码子数目

阅读量:

遇到一道需要现金援助的相关问题:

http://www.timedoo.com/task-id-1194.html

对应的程序代码如下:

复制代码
    #!/usr/bin/env python
    # -*- coding: utf-8 -*-
    __author__ = "Sheng-Wei Ma"
    
    from Bio import SeqIO
    from collections import OrderedDict
    
    records = (r for r in SeqIO.parse('cds.fa', "fasta"))
    count_file = open('cds.count.txt', 'w')
    fre_file = open('cds.frequency.txt', 'w')
    
    count_file.writelines('%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t%s\t

全部评论 (0)

还没有任何评论哟~