序列GC含量计算

blackjimmy

               
               
                import string
def Compute_GC (seq):
    GC = 0
    total =0
    for i in range(0, len(seq)):
        if string.upper(seq) in ['A','T','G','C']:
            total = total +1
        if string.upper(seq) in ['G', 'C']:
            GC = GC+1
    return float(GC)/total*100
seq = 'atcgAATTTCCGGGGTTAATCAGNNNtgGT'
print 'GC percent of %s %2.2f%%' % (seq,Compute_GC(seq))


本文来自ChinaUnix博客，如果查看原文请点：http://blog.chinaunix.net/u1/33851/showart_1890704.html