使用python进行序列的GC含量的统计

在python的IDLE中新建一个文件用于存储：命名为gcgc.py

import sys
with open('E:/cds1.fa','r') as f:#输入文件
    seq=''#用于存放序列
    for line in f:
        line = line.strip('\n')#去除\n符号
        if line.startswith('>'):#如果以>开头，去除>
            name = line[1:]
        else:#序列大写，放在一起
            line = line.upper()
            seq = seq+line
#统计GC含量
g = seq.count('G')
c = seq.count('C')
seq_len = len(seq)#计算序列全长
gc_content = (g+c)/seq_len
print('%s total sequence length :%d' %(name,seq_len))
print('%s GC content: %f' %(name,gc_content))