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 comentarios

Sven
Sven el 27 de Nov. de 2011
You'll need to provide an example of what you mean by sequence.
Walter Roberson
Walter Roberson el 27 de Nov. de 2011
Lucio would seem to have given you the answer in http://www.mathworks.com/matlabcentral/answers/21217-fasta-and-blast
huda nawaf
huda nawaf el 28 de Nov. de 2011
thanks walter, i will see the link
sven , I mean any sequence of integers or characters.
There is algorithms of seq. alignment but designed for DNA or protein sequences.
I hope to find for general seq.
thanks
huda nawaf
huda nawaf el 29 de Nov. de 2011
thanks walter,
i did not see the last comment by Lucio.
really thanks

Iniciar sesión para comentar.

Respuestas (1)

Hin Kwan Wong
Hin Kwan Wong el 28 de Nov. de 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);

Categorías

Más información sobre Genomics and Next Generation Sequencing en Centro de ayuda y File Exchange.

Etiquetas

Preguntada:

el 27 de Nov. de 2011

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by