cytoBand.txt 人类染色体区带文件参考本人文章:
生信技能30 - 获取CNV开始位置和结束位置所在的染色体区带
获取参考基因组各染色体长度参考本人文章:
生信软件8 - bedtools进行窗口划分、窗口GC含量、窗口测序深度和窗口SNP统计
1. 获取染色体开始至结束区带长度
# 存储区带信息文件
cytoband_file = '/path/cytoBand.txt'
# 存储染色体长度文件
length_file = '/path/genome_length.txt'
# 读取基因组染色体长度文件&