Saturday, November 3, 2012

Frequency of Bases

#output frequency of bases for a given sequence
def freq_print(sequence):
    sequence = sequence.upper()
    A = float(sequence.count('A'))
    T = float(sequence.count('T'))
    G = float(sequence.count('G'))
    C = float(sequence.count('C'))
    dict_freq={"A":100.0*(A/len(sequence)), "T":100.0*(T/len(sequence)), "G":100.0*(G/len(sequence)), "C":100.0*(C/len(sequence))}
    return dict_freq

No comments:

Post a Comment