2021年3月13日
如何从vcf文件中批量提取一系列基因的SNP位点？
摘要： 需求 客户的一个简单需求： 我有一批功能基因位点，想从重测序的群体材料中找到这些位点，如何批量快速获得？ 示例文件 gene.txt test.vcf 代码实现 run.sh cat $1 |while read gene chr from to do #echo $chr $from $to if 阅读全文
posted @ 2021-03-13 23:04 生物信息与育种 阅读(5833) 评论(0) 推荐(1) 编辑