VCF文件中QUAL和GQ的区别

最近开始分析vcf文件, 于是去搜了相关VCF格式解读的博客.

大部分关于这两个指标的解读如下,都是描述质量值的,但也没说具体啥区别.

QUAL：Phred格式(Phred_scaled)的质量值，表示在该位点存在variant的可能性；该值越高，则variant的可能性越大；计算方法：Phred值 = -10 * log (1-p) p为variant存在的概率; 通过计算公式可以看出值为10的表示错误概率为0.1，该位点为variant的概率为90%。

GQ：基因型的质量值(Genotype Quality)。Phred格式(Phred_scaled)的质量值，表示在该位点该基因型存在的可能性；该值越高，则Genotype的可能性越大；计算方法：Phred值 = -10 * log (1-p) p为基因型存在的概率.

于是我专门查了下这两者的区别.下面是相关的解答.

个人理解是: QUAL是指变异位点，而GQ是指特定样品的GT。

QUAL告诉您我们对给定站点中是否存在某种变体充满信心。该变化可以存在于一个或多个样本中。
GQ告诉您我们对分配给特定样品的基因型正确的信心有多大。它仅是第二低的PL，因为它是第二低的PL与最低的PL之间的差（始终为0）。

也就是QUAL是针对该位点的,和样本无关.只是说明这个位点存在变异的可能性.而GQ是和样本相关的,它表示该样本基因型的可能性.

VCF文件中QUAL和GQ的区别相关推荐

2021.06.08|提取、比较各样品vcf文件中snp突变频率
目录摘要环境与方法使用代码分析结果总结摘要接到一个wgs项目,要帮助客户统计vcf文件中snp突变频率,比较两个样品的突变位点.这个工作在上一个项目中是手动处理的,当时参考序列短,突变位 ...
Python 导出手机通讯录文件 VCF 文件中的手机号码
文章目录 Python 导出手机通讯录文件 VCF 文件中的手机号码 1.代码 Python 导出手机通讯录文件 VCF 文件中的手机号码 1.代码 if __name__ == '__main__' ...
python提取文本中的手机号_Python从vcf文件中读取手机号并进行去重操作
文章目录 1. Python代码 file = open('test.vcf', 'r', encoding='utf-8') tels = [] for line in file: line = l ...
C++头文件中cstring和string的区别
1.string与cstring有什么区别是C++标准库头文件,包含了拟容器class std::string的声明(不过class string事实上只是basic_string的typedef) ...
shell高级技巧：提取vcf文件中一个contig
这是一个很小众的需求.大部分变异检测都是基于组装质量比较高的基因组,而不是那种初步拼接的contig. 由于初步拼接的参考序列通常会有成千上万个contig序列,也就导致在VCF的头文件的##cont ...
aspx文件中CodeFile与CodeBehind的区别
在ASP.NET中Web编程时,aspx Web 窗体的@page 指令有三个属性(Inherits.CodeFile.CodeBehind),我们在添加aspx页面时,会默认在代码中使用CodeBe ...
基因数据处理77之从vcf文件中提取某条染色体的数据
1.代码: /*** @author xubo*/ package org.gcdss.cli.vcfimport org.apache.spark.{SparkConf, SparkContext} ...
android 代码布局设置wrap_content,Android布局文件中wrap_content和0dp的区别
Layout下的布局设置Widget宽高的填充形式: (1)match_parent:指占满父容器此时要控件的宽或高等于父容器的宽或高. (2)wrap_content和的用法:指控件的高或宽随内容的 ...
latex bib文件中引用期刊与会议区别（会议不显示booktitle字段的问题）
血的教训- 这个参考文献bibtex有几种域名,有的是会议出版物,有的是article,里面的内容叫法不完全一致. 会议是 @inproceedings,期刊是@article 如果会议用了artic ...

VCF文件中QUAL和GQ的区别

VCF文件中QUAL和GQ的区别相关推荐

最新文章

热门文章