2022-06-22VCF中提取指定区域的SNP

awk ' $2>42746030 {print}' R04_R4.vcf | awk '$ 2<42767716 {print $2"\t"$ 4"\t"$5}' >R04_R4_1.vcf

©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

提取目标区域的vcf文件
话不多说，直接上解决办法 vimhouxuan.txt 1110000 #用TAB分割 1100019000...
宗肃書阅读 866评论 0赞 0
2021-01-04根据SNP验证自然选择—XPEHH分析
转载自https://mp.weixin.qq.com/s/1IM6KSCGIq0cUJhKKXVcIw 背景选...
麦冬花儿阅读 14,178评论 3赞 26
2022-03-03 snp gff提取注释信息
查看第9列有哪些注释信息： awk 'BEGIN{FS=OFS="\t"} {split($9, a, ";");...
生信圈阅读 523评论 0赞 0
2022-06-08修改vcf文件样本名称，提取指定vcf样本
1.修改vcf样本名使用bcftools的 rehead指令 bcftools reheader -s samp...
汪大山阅读 8,787评论 6赞 6
实践：awk只提取正则匹配的指定字符串内容
1. 先过滤原始文件中不需要的内容——含有"@_@rm"的行，使用ag 的-v选项过滤掉： ag -v "@_@r...
五大RobertWu伍洋阅读 7,987评论 0赞 0

赞1赞

赞赏

手机看全文