看到一个现金求助的题目:
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%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\n' % ('name', 'TTT', 'TTC', 'TTA', 'TTG', 'CTT', 'CTC', 'CTA', 'CTG', 'ATT', 'ATC', 'ATA', &