grep 练习一

1、需求：根据染色体位置大小提取indel 信息

我的目的基因在chr4 染色体1349548769-1349551625
先将染色体4号的indel 提取出来
grep 'chr4' indel.vcf > chr4_indel.vcf
然后再提取目的区域信息
grep '13495[4-5][0-9][0-9][0-9][0-9]' chr4_indel.vcf > ss_indel.txt

取出的信息上下游包含了都多了几kb 的信息，需要再挑选下，数据比较少就直接鼠标选择一下。

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

推荐阅读更多精彩内容