728x90
반응형
def blosum62_score():
'''
이 기능은 blosum62 score matrix 행렬과 20 x 20]
이를 dictionary 형태로 바꾼 것을 주는 기능이다.
#return
return[0] = 24 x 24 blosum62 matrix
return[1] = 24 x 24 blosum62 dictionary {('amino','acid'):score,......}
'''
import Bio
from Bio.Align import substitution_matrices
mat=substitution_matrices.load("BLOSUM62")
alphabet=list(mat.alphabet)
# DataFrame
df=pd.DataFrame(mat,index=alphabet,columns=alphabet)
# dictionary
dic={}
for aa in alphabet:
for bb in alphabet:
dic[(aa,bb)]=df.loc[aa,bb]
return df,dic
728x90
반응형
'Bioinformatics(생정보학)' 카테고리의 다른 글
Atchley amino acid factor 테이블 (0) | 2024.05.24 |
---|---|
fasta 파일의 특정 영역을 N으로 마스킹 (masking) 하기 (0) | 2024.05.14 |
[de novo assembly] cap3 설치법 (0) | 2024.04.29 |
Samtools 설치법 (0) | 2024.04.23 |
pMHC1 binding benchmark running script (0) | 2024.03.18 |