shell高级技巧：提取vcf文件中一个contig

这是一个很小众的需求。大部分变异检测都是基于组装质量比较高的基因组，而不是那种初步拼接的contig。

由于初步拼接的参考序列通常会有成千上万个contig序列，也就导致在VCF的头文件的##contig=<ID=xxx,length=xxx>部分会有成千上万个contig，将这个文件加载到IGV时, IGV会去解析VCF，这将会是非常缓慢的过程，最好的策略就是只提取其中一个contig查看。

一开始，我想到的是这个方法，直接用bcftools去提取目标contig

bcftools view some_variants.bcf contig_1

但是，你会发现HEADER部分还是会保留原来的信息。怎么办呢？

我们这里会用到一个Shell高级技巧--subshell。举个例子，比如说同时查看一个文件的头10行和后10行.

(head -n 10 ; tail -n10 ) < ~/referece/annotaiton/TAIR10/TAIR10_GFF3_genes.bed

这里的()就是一个subshell, 它的作用就是同时执行()中的所有命令。

用到此处，就是同时完成提取VCF的header部分并只保留目标contig和提取目标区间的所有记录这两个任务，然后保存为vcf.gz格式。

(bcftools view -h some_variants.bcf | grep -v "##contig" | sed "/^##reference/a $(bcftools view -h some_variants.bcf | grep -w "ID=contig_1")" ; bcftools view -H some_variants.bcf contig_1 ) | bcftools view -Oz -o contig_1.vcf.gz

上面是一个非常复杂的shell命令，希望你能看懂。

这里面重复出现了两个变量，一个是some_variants.bcf，一个是contig_1, 我们可以将其用$1和$2代替，然后添加到你的~/.bashrc中，如下所示

# function
subvcf()
{(bcftools view -h $1 | grep -v "##contig" | sed "/^##reference/a $(bcftools view -h $1| grep -w "ID=$2")" ;  bcftools view -H $1 $2 ) | bcftools view -Oz -o $2.vcf.gz && bcftoos index $2.vcf.gz
}

最后用source ~/.bashrc添加刚才的函数到环境变量中，最后就是subvcf some_variants.bcf contig_1 代替之前的长串命令。

shell高级技巧：提取vcf文件中一个contig相关推荐

Tips--利用shell脚本批量提取txt文件中任意字段
利用shell脚本批量提取txt文件中任意字段前言 0. 一个例子 1. cat命令 2. '|'符号与'>'符号 3. grep命令 4. awk命令前言对于测试中出现的log,我们经常 ...
2021.06.08|提取、比较各样品vcf文件中snp突变频率
目录摘要环境与方法使用代码分析结果总结摘要接到一个wgs项目,要帮助客户统计vcf文件中snp突变频率,比较两个样品的突变位点.这个工作在上一个项目中是手动处理的,当时参考序列短,突变位 ...
【Python项目实战】提取.docx文件中的图片并保存到指定的文件夹
文章目录一.需求分析二.系统设计 2.1系统业务流程 2.2系统预览三.系统开发必备 3.1 系统开发环境 3.2文件组织结构四.主函数设计 1.创建窗口 2.创建按钮 3.创建输入框五.函 ...
python读json文件中不同的数据类型_怎么使用python提取json文件中的字段
python中为什么用json有什么作用 python的json模块中如何将变量添加到里面 python的json模块第一个是要打开的文件,第二个是打开的操作,为什么会如果你早认清你在别人心中没那么重 ...
vscode中打开pdf文件_提取pdf文件中的文字
环境说明 windows10系统 python3.6版本安装网上很多说需要安装pdfminer3k和pdfminer3k.six,我尝试了先安装pdfminer3k后安装pdfminer3k.si ...
利用Python提取PDF文件中的文本信息
如何利用Python提取PDF文件中的文本信息日常工作中我们经常会用到pdf格式的文件,大多数情况下是浏览或者编辑pdf信息,但有时候需要提取pdf中的文本,如果是单个文件的话还可以通过复制粘贴来直 ...
如何使用python提取dwg文件中的坐标信息
如果要使用 Python 提取 DWG 文件中的坐标信息,你需要使用专业的 CAD 读取工具,例如 Autodesk AutoCAD.Teigha File Converter 等.这些工具可以将 D ...
提取.c文件中的函数名
shell脚本提取.c生成函数原型写了个脚本主要是为了给自己用的,顺带共享下.作用是提取.c文件中的函数名,生成函数原型. #!/bin/bash sourcesfile=$1 if [[ -f $ ...
提取pdf文件中文字的两种方法
如今,在我们的工作与学习中已经不是单单使用word.Excel等格式文件了,pdf格式的文件已经被广泛地运用到我们的办公室中.大家都知道pdf文件是不可直接编辑与修改的,使用起来有些不便.那么当我们需 ...

shell高级技巧：提取vcf文件中一个contig

shell高级技巧：提取vcf文件中一个contig相关推荐

最新文章

热门文章