因为我想根据基因在染色体上位置获得基因序列，但是发现有一个基因搜索不到序列。其它都可以搜到。
但是这个基因的CDS序列又能通过命令获得。
根据输入的reg.bed文件内容提取对应固定区域（region）序列:seqtk subseq in.fa reg.bed > out.fa
准备bed文件：[基因所在染色体 基因起始位置 基因终止位置]
参考笔记：https://www.omicsclass.com/article/2190
