bio.py 文件源码

def __init__(self, sequence, **kwargs):
        SeqRecord.__init__(self, sequence, **kwargs)

        sequence = upper(sequence)
        g_count = sequence.count('G')
        c_count = sequence.count('C')

        # make GC content as a percent
        content = sum([g_count, c_count]) / float(len(sequence))
        content = content * 100
        self.gc_content = round(content, 0)

        # set the start stop in parent seq
        self.start = int(kwargs.get('id', -1))
        self.stop = self.start + len(sequence) - 1

        # set up score
        self.score = -1
        self.genes = []