ClustalW y Omega de EBI y Blast de NCBI, ambos alinean secuencias globalmente. Smith Waterman de EBI alinea secuencias localmente, pero trabaja con solo dos secuencias. ¿Cómo puedo alinear más de 2 secuencias localmente?
Tienes que hacerlo por parejas (no hay otra opción).
Puedes hacer eso con BLAST
o incluso SSEARCH
(Smith-Watermann). Si tiene un fasta
archivo de sus secuencias, haga esto:
makeblastdb -dbtype nucl -in <yourfastafile> -out <yourdbindex>
BLAST
con su <yourfastafile>
entrada y <yourdbindex>
como base de datos con los parámetros deseados. Esto hará una alineación por pares entre todas las secuencias.Se puede utilizar una estrategia similar para SSEARCH
.
Incluso el agua (en EBI en sí) es una herramienta de alineación local para la alineación de secuencias múltiples. También en BLAST obtienes bla2seq, que también se usa para la alineación por pares. En cualquiera de estas herramientas, solo tiene que ingresar secuencias FASTA (cualquier número que desee) o cualquier secuencia con el signo> y ejecutarla contra el tema (base de datos o su secuencia requerida). Si desea ver los resultados de toda la alineación en una pantalla, una herramienta como BioEdit sería ideal. Espero que esto sea útil, no soy un experto en bioinformática, así que siéntete libre de corregirme :)
Yo uso Jalview, que es un pequeño programa escrito en Java. Es fácil de usar y puede hacer muchas cosas. Está basado en servicios web y tiene 7 servicios diferentes integrados solo para la alineación de secuencias. Puede cargar sus secuencias desde muchos formatos de archivo u obtenerlas por ID de referencia. Es realmente una pequeña navaja suiza. Realmente recomiendo echarle un vistazo. Está disponible en http://www.jalview.org/ .
EDITAR:
Es posible que me haya perdido la parte de alineación local antes, pero investigué un poco y descubrí que MULAN (herramienta de visualización de conservación y alineación local de secuencia múltiple) puede ser lo que está buscando.
jalview
no es para alineaciones locales.
usuario11667
WYSIWYG
blast
en la línea de comandos (terminal UNIX o símbolo del sistema de Windows).blast
no tiene una GUI o incluso si la hay, nunca la he usado. Es mucho mejor escribir comandos (ahorra RAM y también te permite algo de automatización)terdón
WYSIWYG
terdón