http://daler.github.io/pybedtools/
用个下面这个
>>> fn = pybedtools.example_filename('test.fa')
>>> BedTool.seq('chr1:2-10', fn)
'GATGAGTCT'
>>> BedTool.seq(('chr1', 1, 10), fn)
'GATGAGTCT'
http://daler.github.io/pybedtools/
用个下面这个
>>> fn = pybedtools.example_filename('test.fa')
>>> BedTool.seq('chr1:2-10', fn)
'GATGAGTCT'
>>> BedTool.seq(('chr1', 1, 10), fn)
'GATGAGTCT'
转载于:https://www.cnblogs.com/nkwy2012/p/6807491.html