
# manually downloaded Genbank record for KJ660347.2
mkdir -p /hive/data/genomes/eboVir1/bed/genePreds
for i in ../../genbank/*.gb; do f=`basename $i .gb`;  echo s/$f/`basename ../../${f}v*`/g ; done > acc2accVer.sed 

for i in ../../genbank/*.gb
do 
f=`basename $i .gb`
/usr/local/bin/bp_genbank2gff3.pl --split $i
sed 's/iD=/ID=/' $f.gff | gff3ToGenePred stdin  stdout | sed 's/.t01//' | sed -f acc2accVer.sed | sort -nk 4 > $f.gp
done
cat *gp | hgLoadGenePred eboVir2 ncbiGene stdin

# build mafFrames
# genePredToMafFrames eboVir2 /hive/data/genomes/eboVir2/bed/multiz49way/multiz49way.maf frames.tab eboVir2 KJ660347.gp
genePredToMafFrames eboVir2 /hive/data/genomes/eboVir2/bed/multiz49way/multiz49wayStrainNames.maf frames.tab eboVir2 KJ660347.gp
hgLoadMafFrames eboVir2 multiz49wayFrames frames.tab  
