command of alignment for general sequences
이전 댓글 표시
hi, is there a command of sequence alignment for any sequence generally ? I did not mean swalign or nwalign, these just for biological sequences. if not are there anyone wrote a code of sequence alignment for non biological sequences in matlab?
thanks
댓글 수: 4
Sven
2011년 11월 27일
You'll need to provide an example of what you mean by sequence.
Walter Roberson
2011년 11월 27일
Lucio would seem to have given you the answer in http://www.mathworks.com/matlabcentral/answers/21217-fasta-and-blast
huda nawaf
2011년 11월 28일
huda nawaf
2011년 11월 29일
답변 (1개)
Hin Kwan Wong
2011년 11월 28일
What sequence? If you mean any two numerical sequences which are correlated. Their alignment can be done by using cross correlation and finding where the peak is in terms of the delay, then shifting one of the sequence by the delay.
xcorr = ifft(fft(signalA).*conj(fft(signalB));
delay = find((xcorr==max(xcorr))==1);
AlignedA = circshift(signalA,1-delay);
카테고리
도움말 센터 및 File Exchange에서 Genomics and Next Generation Sequencing에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!