command of alignment for general sequences
Ältere Kommentare anzeigen
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 Kommentare
Sven
am 27 Nov. 2011
You'll need to provide an example of what you mean by sequence.
Walter Roberson
am 27 Nov. 2011
Lucio would seem to have given you the answer in http://www.mathworks.com/matlabcentral/answers/21217-fasta-and-blast
huda nawaf
am 28 Nov. 2011
huda nawaf
am 29 Nov. 2011
Antworten (1)
Hin Kwan Wong
am 28 Nov. 2011
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);
Kategorien
Mehr zu Genomics and Next Generation Sequencing finden Sie in Hilfe-Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!