#!/bin/bash
sudo apt-get install exonerate
if [ ! -e Homo_sapiens.GRCh37.58.dna.chromosome.22.fa ]; then
  wget ftp://ftp.ensembl.org/pub/current_fasta/homo_sapiens/dna/Homo_sapiens.GRCh37.58.dna.chromosome.22.fa.gz
fi
if [ ! -e uniprot_sprot.0.fasta ]; then
  wget http://trac.nchc.org.tw/grid/raw-attachment/wiki/jazz/10-08-03/uniprot_sprot.0.fasta
fi
gunzip Homo_sapiens.GRCh37.58.dna.chromosome.22.fa.gz
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'