본문 바로가기
2023 상반기/바이오 컴퓨팅

multiple sequence alignment & 병합

by concho 2023. 5. 2.

1. Merge the pairwise alignment into the multiple sequence alignment.

  PSA의 x MSA의 x  
1.  not '-' not '-' Align this column
2. not '-'  '-' all PSA gap
3.  '-' not '-' all MSA gap
4.  '-'  '-' Align this column

 

PSA

PSA0
x YWQLIK
y D-EHQH
PSA1
x YWQLIK
z H--LAK
PSA2
x -Y-WQLIK
w AAWHL--

MSA

x YWQLIK
y D-EHQH
z H--LAK
   

PSA1

x YWQLIK
z H--LAK

MSA2

x YWQLIK
y D-EHQH
z H--LAK

PSA2

x -YWQLIK
w AAWHL--

MSA3

x -YWQLIK
y -D-EHQH
z -H--LAK
w  

 

 


x (center) ATCGATCG 
y ATTCGCCG
z ATCGAATT
w AACGTCGG

PSA

x AT-CGATC-G
y ATTCG--CCG
   
x ATCGA-TCG
z ATCGAATT-
   
x ATCGATC-G
w AACG-TCGG

new PSA

MSA
x1: AT-CGATC-G
y: ATTCG--CCG 

 

x1(M)과 x2(P)를 비교해서 z or P에 테이블대로 gap 삽입 혹은 그대로 추가 or MSA에 gap추가

 

MSA
x1: AT-CGA-TC-G gap 추가
y: ATTCG---CCG gap 추가
z: AT-CGAATT-- 추가

x1(M)와 x3(P)를 비교해서 w or P에 테이블대로 gap 삽입 혹은 그대로 추가 or MSA에 gap추가

 

MSA
x: AT-CGA-TC-G

y: ATTCG---CCG

z: AT-CGAATT--

w: AA-CG--TCGG

 

MSA

x1 AT-CGA-TC-G
             3
y ATTCG---CCG
             3
z AT-CGAATT--
1121 113 11
w AA-CG--TCGG

 

댓글