MATLAB Answers

Alignment(3,:) in seqinsertgaps

3 views (last 30 days)
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


Sign in to comment.

Accepted Answer

Tim DeFreitas
Tim DeFreitas on 30 Nov 2018
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
' ||||'
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
madhan ravi on 4 Dec 2018
Bio Ing Sapienza's answer moved here for consistency:

Sign in to comment.

More Answers (0)

Translated by