NCBI taxonomy ID (any level) to genomes (assemblies)

¿Necesitas obtener los genomas (ensambles) de todas las especies asociadas a algún taxid (cualquier nivel taxonómico)?, este es tu post. Primero es necesario tener el taxid que nos interesa. En este caso usaré el taxid de una familia de virus llamada Birnaviridae. Si seguimos este link, veremos que a la fecha existen 23 ensambles asociados […]

Transformando formatos de audio

Hace unos días me compartieron unos audios en formato WEBM, los cuales puedo abrir con el reproductor Media Player Classic pero no con mi reproductor favorito Winamp o en otros dispositivos electrónicos. En resumen, necesitaba transformar WEBM a MP3. Revisé en internet y existen distintos programas gratuitos que ofrecen realizar el trabajo pero no me […]

Linux: Ejecutar un comando con el argumento del comando anterior

Has ejecutado un comando en linux (bash) y quieres ejecutar otro comando usando el "argumento" del comando anterior, acá te dejo algunos trucos. La mejor respuesta tomada desde este link. Most likely you'd want to go in the recently created directory. You can do this using next "shortcuts": cd Esc. - type cd and after press Esc followed by . (not in the […]

Tips para procesar archivo fastq

Estos son algunos TIPs para manipular archivos fastq Concatenar por ID El comando "paste" permite "contatenar" el contenido del fastq de cada ID en una sola línea, es decir, las cuatro líneas de cada ID quedarán "unidas" en una sola línea, separadas por tabulador. Obtener lista de IDs Usando el comando "paste", junto con otros […]

Transformar anotación GOs de genes a un archivo GAF

Para realizar un análisis de enriquecimiento funcional necesitamos en muchos casos un archivo GAF (por ejemplo en ontologizer). Acá te explico la forma de obtenerlo si cuentas con la anotación GO de tus genes. 1. ¿Qué es un archivo GAF? El archivo GAF (Gene Association File) es una forma de almacenar la información de anotación funcional de […]

No prodría vivir sin... estos programas

Mi lista de imprescindibles que completaré día a día Cluster: Windows & Tab manager extention for Chrome Grammarly STAMP for microbiome and functional analysis Cytoscape (java problem? Direct the Cytoscape installer to uncompressing folder of jdk-11.0.12_windows-x64_bin.zip java (other link))

Consejos para alineamiento múltiple de secuencias

Estos son mis TIPs para situaciones de alineamiento múltiple de secuencias Orientación de secuencias Algo extremadamente impotante a la hora de realizar un alineamiento múltiple es que todas las secuencias estén en la misma orientación. Cuando son 2 a 3 secuencias, la dirección se puede "ajustar" de forma manual mediante un alineamiento con blast (align2seqs) […]

Setup de datos para STAMP

Quienes conocen STAMP saben lo útil que es para analizar datos de microbioma y similares. Lamentablemente no es muy popular debido a que no es fácil preparar el input de datos. Acá escribiré mis consejos para preparar el input y para usar STAMP. Requisitos Debes tener un Excel con la tabla de abundancia y taxonomía […]

R / RStudio comandos útiles

Lista de comandos salvadores Información de sesiónsessionInfo() Listar paths de librerías.libPaths() Librería del usuarioSys.getenv("R_LIBS_USER") Listado de paquetes instalados y sus ubicacionesp <- available.packages()i <- installed.packages() # goodView(p) Actualización de paquetes BioconductorBiocManager::install()Recomiendo primero instalar todo "a". Algunos paquetes arrojarán error por lo cual en una segunda oportunidad se vuelve a ejeutar el comando pero seleccionas los […]