def get_seq_GC(seq):
    """
    From a nucleotide sequence return the GC content
    """
    seq = seq.lower()
    n_C = seq.count('c')
    n_G = seq.count('g')
    GC_content = (n_C + n_G) / float(len(seq))
    return GC_content
