计算人类基因组外显子区域长度

问题

问题：https://mp.weixin.qq.com/s/xS-h0TQnCwjvgcRA7Ru3NQ

思路：

数据：ftp://ftp.ncbi.nlm.nih.gov/pub/CCDS/current_human/CCDS.current.txt

获取每一个exon的每个碱基坐标。去重重复坐标，然后计数。

## count.py

base_pos = set()   ## 以集合来存储所有非重复位置位置

with open("./CCDS.current.txt", "r") as f:
    fo