SnpSift 可以分割染色體,截取部分vcf文件。
具體使用方法見:
https://blog.csdn.net/weixin_43569478/article/details/108079138?spm=1001.2101.3001.6650.3&utm_medium=distribute.pc_relevant.none-task-blog-2%7Edefault%7ECTRLIST%7ERate-3-108079138-blog-108079137.pc_relevant_recovery_v2&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2%7Edefault%7ECTRLIST%7ERate-3-108079138-blog-108079137.pc_relevant_recovery_v2&utm_relevant_index=6
但其安裝容易遇到Java版本問題。
可以新建一個conda 環(huán)境,然后安裝SnpSift
conda install -c bioconda snpsift
參考:
https://zhuanlan.zhihu.com/p/476561285
使用方法:
zcat Chr9.SampleRename.maffiltered.vcf.gz | SnpSift filter "( POS > 123456 ) & ( POS < 654321 )" > filtered.vcf
zcat FinalSNP.SampleRename.maffiltered.Chr.vcf.gz | SnpSift intervals my_intervals.bed > LI.vcf
my_intervals.bed 格式
chr start end
最后編輯于 :
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。