python学习——把计算GC含量的代码封装成函数

把代码封装成函数的好处是可以重复使用该段代码，并且会使代码结构清晰

例如要计算chr1以及chr2染色体的GC含量，代码如下：

# 将代码封装为函数并重复使用，例如计算染色体的GC含量
chr1 = 'ATCGCGATCCGAATGCAGGAatcagagataccgaNNNNNNttAATGAC'
chr2 = 'agataCGATGNNGATcCCGGGTTATAGCTAGCNNNGGG'
def get_gc_content(chr):
    chr = chr.upper() #也可是使用.lower()把大写转换成小写计算
    count_a = chr.count('A')
    count_t = chr.count('T')
    count_c = chr.count('C')
    count_g = chr.count('G')
    gc_content = (count_g + count_c) / (count_a + count_c + count_g + count_t)
    return gc_content

chr1_gc_content = get_gc_content(chr1) # 计算chr1的GC含量
chr2_gc_content = get_gc_content(chr2) # 计算chr2的GC含量

print('chr1 gc content : {:.2%}'.format(chr1_gc_content))
print('chr2 gc content : {:.2%}'.format(chr2_gc_content))