cd /hive/data/genomes/wuhCor1/bed/ncbiGene
 
# get Genbank entry NC_045512.2.gbk
curl -s  "https://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nucleotide&id=NC_045512.2&rettype=gb&retmode=txt"; > NC_045512.2.gbk
genbankToBigGenePred NC_045512.2.gbk NC_045512v2 ../../chrom.sizes tab bb/
 
# these are not necessary anymore, genbankToBigGenePred now includes them
#bbedToBigBed -allow1bpOverlap -tab -type=bed12+ -as=as prev.bed ../../chrom.sizes bb/genes.bb 
#bedToBigBed -allow1bpOverlap -tab -type=bed12+ -as=as prev.bed ../../chrom.sizes bb/genes.bb   -extraIndex=name

tawk '{print $4,   $21, $23, $22, $22}' prev.bed | sed 's/\.[0-9]*$//' > ncbiSearch.txt
ixIxx ncbiSearch.txt ncbiSearch.ix ncbiSearch.ixx
ln -s `pwd`/ncbiSearch.ix /gbdb/wuhCor1
ln -s `pwd`/ncbiSearch.ixx /gbdb/wuhCor1
