Acesso
remoto a servidores e introdução ao linux
Configurar
servidor (TAREFA)
Sistema
operacional:
- CentOS (distribuição gratuita do RedHat Enterprize)
- Fedora
- Ubuntu
Como se
conectar a partir de Windows:
- Troca de arquivos: Secure
Shell
- Conexão direta SSH: Putty
Comandos
Linux:
- Veja aqui
Entrando
numa conta no servidor:
- Abra o programa Putty
- Complete o campo Host name com maracatu.icb.ufmg.br e mude Protocol para SSH
- Complete o campo Saved sessions com maracatu
- Clique o botão Save
- De agora em diante clicando
maracatu o terminal abre
- Logue como curso, pegue a senha na aula
- Vamos exercitar alguns comandos
- $mkdir
seunome (não use acentos, maiúsculas, espaços)
- $ls
– lista o conteúdo
- $cd seunome
- $mkdir
teste
- $ls
- $cd .. (olha os dois pontos
ai gente!!!)
Editor
de texto online (vi):
- $ls
/home/treinamento
- $more /home/treinamento/aldolase.cds
- $cd seunome
- $cp
/home/treinamento/aldolase.cds . [opag] pergunte na aula sobre o “dedinho do bioinformata”
- $ls
- Nesse momento vc já ta se virando bem, falta só um poder: editar o aquivo lá no servidor
- $vi aldolase.cds
- Para entrar no INSERT MODE
digite “i”
- Troque o nome da seqüência para
>teste
- Pra salvar tecle ESC + dois
pontos (:) e x!
- $more aldolase.cds
- $mv aldolase.cds aldolase.hsa
- Belo treino!
- Agora pegue uma seqüência de mioglobina no NCBI (veja como na aula)
- Send to File
- Abra com bloco de notas
- Digite no terminal: $vi mioglobina
- Copie do bloco de notas
- Botão da direita no terminal
cola!!!!!!!!
- Sai como antes EXC + : + x!
- Agora vc
sabe tudo!
O melhor
de todos os comandos: grep
- $ls
/home/treinamento/db
- $more /home/treinamento/db/hsa.ORESTES (aperte
“q” para sair)
- $cat /home/treinamento/db/hsa.ORESTES | grep “>” |
more
- $cat /home/treinamento/db/hsa.ORESTES
| grep “>” –c
- Pronto agora vc sabe contar seqüências num arquivo
BLAST no NCBI: http://www.ncbi.nlm.nih.gov/
- Clique em BLAST e em nucleotide BLAST
- Mude para Others (nr etc) – acham que todo mundo só
analisa humanos...
- Mude para Somewhat similar sequences (blastn)
- $more aldolase.hsa selecione com o mouse e cole na janela do BLAST
- Rode o programa, explore a
coluna score, clique nela, olhe o E-value
- E-value é o número de alinhamentos
melhores que apareceriam sem homologia, para sentir isso, submeta cinco
linhas de bases aleatórias no formato FASTA:
>seunome
acgatcgatcgatcgatcgatcgtagctacgtacg....
Analisando
ORESTES:
- BLAST: $blastall
–p programa –i sequencias –d base –e 1e10 –F F
–b 1 –a 4 –m 8 –o saida_m8 &
- p=programa: tblastn,
etc.
- i=entrada: query
- d=database: tem que estar formatada
- e=E-value: o número de alinhamentos iguais ou melhores que pode ser obtido
sem homologia
- F: low
complexity filter,
- b: o número de alinhamentos
reportados para cada query
- a: número de processadores,
use –a 4
- m: é o tipo de saída, teste –m 8
e –m 9. –m 0 é a saída web
- o: é o nome que vc quer para a saída
- &: manda o job
para o background
- Quando acabar de
colocar pra rodar de um “top” para acompanhar.
- $blastall
–p blastn –i aldolase.hsa –d /home/treinamento/db/hsa.ORESTES –e 1e-10 –F F –b 1000000000 –a 4 –o saida &
- $blastall
–p blastn –i aldolase.hsa –d /home/treinamento/db/hsa.ORESTES –e 1e-10 –F F –b 1000000000 –a 4 –m 8 –o
saida_m8 &
- Compare as duas saídas.
- $cat
saída_m8 | awk ‘{print
$2}’ | sort | uniq | wc
Analisando
a glicólise (idem acima com –m 8)
- $blastall
–p tblastn –i /home/treinamento/db/hsa.glicolise –d /home/treinamento/db/hsa.ORESTES –e 1e-10 –F F –b 1000000000 –a 4 –m 8 –o
saida_m8_glicolise &
- $cat
saída_m8_glicolise | awk
‘{print $1}’ | sort | uniq –c | more
E-mail Miguel