MATLAB Answers

0

Alignment(3,:) in seqinsertgaps

Asked by Bio_Ing_Sapienza on 30 Nov 2018
Latest activity Commented on by madhan ravi
on 4 Dec 2018
Hi guys. There's someone able to explain me why in seqinsertgaps for the alignment (al) I have to choose the first and the third line?i.e al(1,:)
al(3,:) What should be the general criterion? Thanks a lotmatlabimage.png

  0 Comments

Sign in to comment.

1 Answer

Answer by Tim DeFreitas on 30 Nov 2018
 Accepted Answer

The alignment returned by nwalign is a 3xN character array showing the optimal alignment between the two input protein sequences. The first row is the first input sequence, the middle row shows alignment information, and the last shows the second sequence:
>> al(:,1:10)
ans =
3×10 char array
'QKQEIKMNPN'
' ||||'
'------MNPN'
We then use seqinsertgaps to get the original nucleotide sequence, but with the gap information returned by the alignment, so the inputs are the first and third rows of the alignment character vector.

  1 Comment

madhan ravi
on 4 Dec 2018
Bio Ing Sapienza's answer moved here for consistency:
Great!

Sign in to comment.