fastaから配列を取り出す

$ samtools faidx genome.fasta chr1:15,000-50,000
こんな感じでfastaファイルから取り出したい配列を指定して取り出すことが可能。
ある遺伝子を特定したとして、その上流プロモーターの配列をfastaから取り出したいなどの用途で使える。

また
$ samtools faidx genome.fasta genecode1 genecode2 genecode3....
とすると該当するgenecodeのfastaを集めてくれる。