1. 目的:提取十号染色体特定碱基位点的序列
2. 准备文件：基因组文件、特定位点文件

bedtools getfasta -fi Tgra.chr.fa -bed weidian.txt -fo output.fa

3. 结果文件

4. blast检查序列是否提取正确

makeblastdb -in my_sequences.fasta -dbtype nucl -out my_db #建库
blastn -query query.fasta -db nt -out result.txt -evalue 1e-5 -outfmt 6 #核酸比对

5. Bedtools提取序列时，需要让起始位点-1（结论来自睿远论坛）

睿远论坛链接：https://www.kdocs.cn/l/cmP1mhDTHznF

原文链接：Bedtools提取特定片段的序列
https://www.kdocs.cn/l/ca011D4qhpwF