计算基因组外显子长度

 

下载基因组外显子信心

网站 ftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/current_human/

wget ftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/current_human/CCDS.current.txt

运行下列代码 得到外显子大约36M

import re

import os

from collections import OrderedDict

from operator import itemgetter

 

os.chdir( '/Users/yangqin/Desktop/biotree/python/' )

 

exonLength = 0 # 外显子长度尚未计算

overlapExons = OrderedDict()

with open