Usando o RGS

-
Antes de rodar o RGS, carreguei um arquivo de contigs: 
contigs=fastaread('contigs_frag_k19_opt.fa');
Carregando e rodando o RGS: 
load dbstruct
anotado = sila(dbstruct,contigs,'HGF',true);
an2 = annotmat2fascov(anotado);

Comentários

Postar um comentário

Postagens mais visitadas deste blog

Jellyfish script

-
Eu precisei criar um histograma dos k-mers para diferentes tamanhos. O jellyfish surpreende em performance e facilidade de uso. Mas mesmo assim, precisamos de vários passos. Segue um script para facilitar o uso dele. O script precisa de dois parâmetros: o nome do arquivo fasta/fastq e o prefixo do arquivo de resultados. Ex: run.sh reads.fasta jf_kmer_ #!/bin/bash for kmerSize in {13..29..2} do    echo "kmer size $kmerSize - creating kmer-map..."    jellyfish count -m $kmerSize -o output -c 16 -s 500000000 -t 60 $1    arq_out=`ls -1 output* | wc -l`    banco=out_$kmerSize    if [ $arq_out -gt 1 ]; then       echo "merging"       jellyfish merge -o $banco output\*       rm -rf output*       sleep 5    else       echo "renaming..."       auxFileName=`ls -1 output*`       mv $auxFileName $banco    fi   ...
Leia mais

Conversão do encode do MariaDB para atender o moodle 3.8

-
A nova versão do moodle precisa do encode utf8mb4 para codificação os emoticons, após a atualização da versão 3.2 para a 3.8, tive que rever as tabelas no mariadb. Uma consulta simples para verificar a codificação das tabelas: select table_schema,table_name,table_collation from information_schema.TABLES where TABLE_SCHEMA='NOME_DO_BANCO'; E a rotina para converter as mais de 360 tabelas do moodle... DELIMITER // CREATE PROCEDURE convert_encode() BEGIN    DECLARE _stmt CHAR(255);    DECLARE cmd CHAR(255);    DECLARE done BOOL DEFAULT 0;    DECLARE c1 CURSOR FOR         SELECT CONCAT("ALTER TABLE ", TABLE_SCHEMA, '.', TABLE_NAME,         ' CONVERT TO CHARACTER SET "utf8mb4" COLLATE "utf8mb4_bin";')         FROM INFORMATION_SCHEMA.TABLES WHERE TABLE_SCHEMA='NOME_DO_BANCO' ;    DECLARE CONTINUE HANDLER FOR NOT FOUND SET done = 1;    OPEN c1; ...
Leia mais

O GBParsy é uma biblioteca para realizar o parser de arquivos GenBank para o Python

-
O GBParsy é uma biblioteca para realizar o parser de arquivos GenBank para o Python. No Ubuntu para realizar a instalação do módulo eu segui as instruções http://code.google.com/p/gbfp/wiki/HowToInstall até o ponto em que a instrução "$ python setup.py build " produziu o erro de ausência do arquivo "python.h". Para resolvê-la, sudo apt-get install python-dev site do projeto: http://code.google.com/p/gbfp/ Artigo: http://www.biomedcentral.com/1471-2105/9/321
Leia mais