grep 練習一

1、需求：根據(jù)染色體位置大小提取indel 信息

我的目的基因在chr4 染色體1349548769-1349551625
先將染色體4號的indel 提取出來
grep 'chr4' indel.vcf > chr4_indel.vcf
然后再提取目的區(qū)域信息
grep '13495[4-5][0-9][0-9][0-9][0-9]' chr4_indel.vcf > ss_indel.txt

取出的信息上下游包含了都多了幾kb 的信息，需要再挑選下，數(shù)據(jù)比較少就直接鼠標選擇一下。

?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
【社區(qū)內容提示】社區(qū)部分內容疑似由AI輔助生成，瀏覽時請結合常識與多方信息審慎甄別。
平臺聲明：文章內容（如有圖片或視頻亦包括在內）由作者上傳并發(fā)布，文章內容僅代表作者本人觀點，簡書系信息發(fā)布平臺，僅提供信息存儲服務。

相關閱讀更多精彩內容

友情鏈接更多精彩內容