第一次写博客，分享一个做的提取基因序列的程序，根据bed文件里的位置信息从基因组里提取序列

源码地址：https://github.com/Liuyuan2018/fastaTools/blob/master/pyGetFasta.py

SRE实战 互联网时代守护先锋，助力企业售后服务体系运筹帷幄！一键直达领取阿里云限量特价优惠。

bed文件通常用来保存注释基因信息，BED文件必须的3列:

  1. chrom - 染色体号
  2. chromStart - feature在染色体上起始位置（其实编号为0）
  3. chromEnd - feature在染色体上末尾位置（不包括此编号）

　　第四列是基因的名称

　　还有些列想了解参考:http://genome.ucsc.edu/FAQ/FAQformat.html#format1

程序依赖 pyfasta模块（https://pypi.org/project/pyfasta/）

安装pyfasta的命令：pip install pyfasta

 

扫码关注我们
微信号：SRE实战
拒绝背锅 运筹帷幄