首先要下载并且得到人类基因组的外显子坐标记录文件
这里我用的参考基因组版本仍然是hg19，所以去CCDS数据库里面下载对应版本，并且格式化成BED文件。
wget ftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/archive/Hs37.3/CCDS.20110907.txt cat CCDS.20110907.txt |perl -alne '{/[(.*?)]/;next unless $1;$exons=$1;$exons=~s/\s//g;$exons=~s/-/\t/g;print "$F[0]\t$" foreach split/,/,$exons;}' >hg19exon.bed