| 1 | #!/bin/bash |
|---|
| 2 | sudo apt-get install exonerate |
|---|
| 3 | if [ ! -e Homo_sapiens.GRCh37.58.dna.chromosome.22.fa ]; then |
|---|
| 4 | wget ftp://ftp.ensembl.org/pub/current_fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.58.dna.chromosome.22.fa.gz |
|---|
| 5 | fi |
|---|
| 6 | if [ ! -e uniprot_sprot.0.fasta ]; then |
|---|
| 7 | wget http://trac.nchc.org.tw/grid/raw-attachment/wiki/jazz/10-08-03/uniprot_sprot.0.fasta |
|---|
| 8 | fi |
|---|
| 9 | gunzip Homo_sapiens.GRCh37.58.dna.chromosome.22.fa.gz |
|---|
| 10 | exonerate uniprot_sprot.0.fasta Homo_sapiens.GRCh37.58.dna.chromosome.22.fa --model protein2genome --showsugar false --showvulgar false --showalignment false --ryo 'RESULT %S %pi %ql %tl %g %V\n' |
|---|
