# parse uniprot
uniprotToTab /hive/data/outside/uniProtCurrent/  --taxon 128952 eboVirUniprot/
# create a mapping from uniprot prot seq to DNA. The mapping has to be in nucl
# space, so use pslProtCnv
blat -t=dnax -q=prot /hive/data/genomes/eboVir3/eboVir3.2bit eboVirUniprot/uniprot.128952.fa.gz stdout -noHead | sort -k1,1n -r | pslCDnaFilter stdin stdout -minId=0.95 -minCover=0.95 | pslProtCnv > spToeboVir3.psl
# lift the uniprot annotations
uniprotLift eboVirUniprot/128952 eboVir3 spToeboVir3.psl 
# copy over output files
ln -s /cluster/data/eboVir3/bed/uniprot/spMut.bb /gbdb/eboVir3/bbi/spMut.bb
ln -s /cluster/data/eboVir3/bed/uniprot/spAnnot.bb /gbdb/eboVir3/bbi/spAnnot.bb
ln -s /cluster/data/eboVir3/bed/uniprot/spStruct.bb /gbdb/eboVir3/bbi/spStruct.bb
hgBbiDbLink eboVir3 spMut /gbdb/eboVir3/bbi/spMut.bb
hgBbiDbLink eboVir3 spAnnot /gbdb/eboVir3/bbi/spAnnot.bb
hgBbiDbLink eboVir3 spStruct /gbdb/eboVir3/bbi/spStruct.bb
