bedtools 批量提取snp前后1000bp片段

本文介绍如何利用bedtools从SNP位置提取其前后1000bp的DNA序列。首先需要一个包含SNP染色体、位置信息的.txt文件,然后通过bedtools命令生成sig_sequence文件。
摘要由CSDN通过智能技术生成

准备条件
1 知晓snp位置,位于第 x 条染色体上 12345678 bp
2 生成 .txt 文件,文件格式如下

txt 文件格式

第一列表示 snp 染色体位置
第二列表示 snp 前 500 bp 位置
第三列表示 snp 后 500 bp 位置
要使用 tab 制表符分开
将文件命名(这里命名为sigpointforsequence.txt)
执行下面命令

bedtools getfasta -fi $GENOME(基因组绝对路径) -bed sigpointforsequence.txt -fo sig_sequence

得到 sig_sequence 文件如下

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值