Foro - Perl en Español

Hola a todos.

Deseo saber si alguien de uds. ha usado algún script en Perl para ejecutar Clustalw y realizar algunas tareas como: cambiar formato de archivos (.pir a msf), realizar alineamientos, etc.

¿Saben dónde puedo encontrar algunos ejemplos?

Les agradeceré mucho cualquier información relacionada.

Gracias.


P.

Si ese programa se puede ejecutar desde la línea de comandos, entonces puedes ejecutarlo desde Perl usando las comillas invertidas (``), las funciones qx(), system() y exec().

Depende un poco de qué resultado quieres obtener. Si solo es ejecutarlo y terminar, pues exec(). Si quieres ejecutarlo y seguir con el programa Perl, pues system(). Y si quieres además obtener la salida devuelta por Clustaw, pues entonces las comillas invertidas o la función qx().

¿Puedes poner un ejemplo?

Me explico:

Tengo instalado clustalw en el sistema, el cual puedo invocar desde la linea de comandos con un solo $ clustalw y que he utilizado para cambiar el tipo formato de un archivo, de fichero.pir crear fichero.msf.

Deseo tener un script en Perl para automatizar este mismo cambio de formato de una carpeta con miles de archivos .pir. El problema es que no sé como abordarlo: dentro de clustalw, para generar un fichero .msf tengo que ir seleccionando de manera interactiva con clustalw las opciones que realizan la tarea de cambio de formato.

Por ejemplo, clustalw me pide:

Código: Seleccionar todo

your choice:

y tecleo el número 1 y un enter, pues en el menú que se despliega también éste corresponde a la opción que permite acceder a cambio de formato. En seguida me aparece

Código: Seleccionar todo

name of file:

y tecleo el nombre del archivo que deseo formatear más un enter. En otro punto me aparece

Código: Seleccionar todo

your choice

y tecleo 4 más un enter, pues éste corresponde dentro del menú desplegado a la opción de formato .msf.

En otro punto

Código: Seleccionar todo

enter name of new file:

y tecleo el nombre que deseo ponerle al archivo que generará en formatdo .msf, que por cierto será el mismo que el original, es decir, fichero.msf.

En otros puntos del proceso tengo que dar unos enter o teclear números que representan, nuevamente, opciones.

No sé si me dí a entender, pero esa es la tarea básicamente.

Gracias.

P.

Leyendo la documentación de clustalw, al final de ella, aparecen los parámetros que puedes indicar en la línea de comandos, por lo que no es necesario que te pregunte.

Intenta ejecutarlo de esa manera, a través de parámetros puestos en la línea de comandos, al lado del ejecutable clustalw, algo así:

Código: Seleccionar todo

clustalw -OUTPUT=GCG -INFILE=fichero.pir (resto de opciones)

Luego, desde Perl, es fácil hacer un bucle recorriendo los ficheros a transformar:

¡Hola!

He averiguado que los siguientes parámetros hacen lo que deseo en clustalw y lo refiero a un archivo en particular:

El comando anterior crea el archivo P0126.msf en el mismo directorio y con el contenido que deseo.

Desearía usar comillas en lugar de system() ya que quiero tener la salida devuelta por clustalw pues a cada uno de los archivos .msf les eliminaré la primera línea.

Los archivos .pir están en una carpeta (llamada PIR_master) y desearía que los que se generen se vayan a la carpeta MSF_master que ya he creado en el mismo directorio.

Gracias

A ver... que no me aclaro... dices que con ese comando te crea un fichero en el mismo directorio. Y al mismo tiempo quieres modificarlo quitándole una línea del principio, que te devuelve la misma ejecución del comando.

Me parece que me he perdido. ¿El comando hace las dos cosas a la vez, crear el fichero y devolverlo por la salida estándar?

Sí, el comando hace esas dos cosas pero ahora a cada archivo .msf que vaya generando a partir de cada uno de los archivos .pir, deseo abrir cada uno para quitarle la primera línea de texto ya que interfiere con lo que deseo hacer con cada uno de ellos.

He estado leyendo y estoy pensando en algo así como

Me he querido referir a un archivo en particular, pero ¿cómo lo hago extensivo a los miles de archivos?

Si quieres editar el fichero que se ha grabado en disco, entonces no necesitas leer nada de la salida estándar devuelta por clustalw.

Una forma de hacerlo sería hacer el bucle, transformar cada uno, y leer el resultado, modificarlo, y escribirlo en el directorio destino.

Algo así (no probado):

Habría un desajuste en el programa si el clustalw crea el fichero en el directorio actual y no en el de trabajo. Yo probaría usando el parámetro OUTFILE de clustalw para indicar dónde tiene que guardar el resultado.