提取指定区域基因组序列
- 单个序列,samtools faidx 提?。?br>
samtools faidx reference.fa chr1:1000-1100
- 多个区域序列提取,构建区域bed文件, 利用bedtools提?。?br>
bedtools getfasta -fi reference.fa -bed target.bed -fo output.fa
samtools faidx reference.fa chr1:1000-1100
bedtools getfasta -fi reference.fa -bed target.bed -fo output.fa