[python] sequence logo 그리는 법

import logomaker
# 예제용 무작위 서열 정보 생성
import random
def generate_random_9mer_sequence(n,length=9):
    '''
    this function generate random 9mer peptide sequence 
    '''
    aas=list('ACDEFGHIKLMNPQRSTVWY')
    # Generate n random 9-mer peptide sequences
    peptides = set()  # Use a set to store unique peptides
    while len(peptides) < n:
        peptide = ''.join(random.choice(aas) for _ in range(length))
        peptides.add(peptide)
    return list(peptides)
    
random_seq=generate_random_9mer_sequence(n=1000,length=9)

# PSSM 생성
pssm=logomaker.alignment_to_matrix(sequences=random_seq,to_type='information',pseudocount=0.1)
logomaker.Logo(pssm)