vcf文件的index

GATK在处理vcf文件的时候，需要vcf文件有index文件。要么是*.vcf.idx，要么是.vcf.gz.tbi.

假如一个vcf文件名为 “a.vcf”; .

1. 可以用“a.vcf”生成 “a.vcf.idx”,方法如下：

igvtools index a.vcf #其实很多软件都可以做到

2. 可以用生成“a.vcf.gz.tbi”,方法如下：

bgzip -c a.vcf >a.vcf.gz #这里必须先用bgzip压缩

tabix -p vcf a.vcf.gz

但实际上，现在GATK好像逐渐抛弃了idx格式，更多的建议使用tbi格式。

vcf文件的index相关推荐

python 处理snp的vcf文件，统计snp在基因的intron、exon还是上游、下游还是不在基因及基因附近
1.c处理vcf文件,初步统计snp位置位置信息有5种 down_stream是gene下游2k以内,up_stream 是gene上游2k以内,gene_exon是snp在外显子内,gene_in ...
R语言丨根据VCF文件自动填充对其变异位点并生成序列fa文件
根据VCF文件自动填充对其变异位点并生成序列fa文件首先提出一个问题: 假如有一个重测序结果VCF文件,里面包含了很多个样本在几百个突变位点(snp和iad)的基因型数据,现在想根据这份原始数据,得 ...
实操 | 合并VCF文件的几种方法及注意事项
背景在基因组分析领域的很多不同场景中,需要合并VCF文件. VCF文件.简单来说,就是记录样本基因型的文件.但多数VCF文件不只记录了基因型,也包含有关该基因型的来源的细节. 其它文件.VCF文件 ...
shell高级技巧：提取vcf文件中一个contig
这是一个很小众的需求.大部分变异检测都是基于组装质量比较高的基因组,而不是那种初步拼接的contig. 由于初步拼接的参考序列通常会有成千上万个contig序列,也就导致在VCF的头文件的##cont ...
java利用io流读取文件发生了乱码以及java解析vcf文件
如果文件中含有中文字符,那么这时候直接进行读取,即new FileInputStream(要读取的文件),就会出现了乱码,所以这时候,我们需要设置对应的格式进行读取,所以还需要传递对应的编码,所以将代 ...
Git 索引文件（index file）
这次重点讲解索引文件(index file)的作用. 我们在提交工作时,使用最多的命令就是 git commit -a 了,但是这个将提交你所做的所有工作.其实,如果你了解 commit 的工作机制, ...
pta通讯录排序用python实现,python实现将android手机通讯录vcf文件转化为csv
经常会遇到将手机通讯录导出到电脑并转化为在电脑中可编辑的情况,在网上搜索了很久当前不外乎两种处理方式.1.使用电脑的outlook的通讯簿功能,将手机导出的vcf文件导入到outlook的通讯录中,然 ...
blog.mm index.php,每天一个WordPress文件：index.php
第一个当然是介绍 index.php,WordPress 根目录下面的 index.php 是 WordPress 的最重要的一个文件,根据基本的 htaccess 规则: # BEGIN WordP ...
python提取文本中的手机号_Python从vcf文件中读取手机号并进行去重操作
文章目录 1. Python代码 file = open('test.vcf', 'r', encoding='utf-8') tels = [] for line in file: line = l ...
2021.06.08|提取、比较各样品vcf文件中snp突变频率
目录摘要环境与方法使用代码分析结果总结摘要接到一个wgs项目,要帮助客户统计vcf文件中snp突变频率,比较两个样品的突变位点.这个工作在上一个项目中是手动处理的,当时参考序列短,突变位 ...

vcf文件的index

vcf文件的index相关推荐

最新文章

热门文章