#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
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