首页 > 其他分享 >变异等位基因频率VAF计算方法

变异等位基因频率VAF计算方法

时间:2024-11-12 15:44:55浏览次数:3  
标签:变异 VAF 等位基因 reads 测序 位点 计算方法

在生物信息学中,变异等位基因频率(Variant 等位基因频率,VAF)是指在一个种群中,某个特定变异位点上变异等位基因出现的频率,常见的计算 VAF 的方法如下:

1. 基于测序深度的计算方法

  • 原理:通过统计变异位点的测序深度以及变异等位基因的支持 reads 数来计算 VAF。测序深度是指在该位点上所有覆盖到的 reads 数量,而变异等位基因的支持 reads 数是指与参考基因组不同的等位基因的 reads 数量
  • 计算公式:VAF = 变异等位基因的支持reads数量/ 测序深度
  • 例子:假设在某个变异位点上,测序深度为 100X,其中有 20 条 reads 支持变异等位基因,则 VAF = 20 / 100 = 0.2

2. 基于基因型数据的计算方法

  • 原理:如果已经有了个体的基因型数据,可以直接根据基因型来计算 VAF。对于二倍体生物,每个个体在每个位点上有两个等位基因。
  • 计算公式
    纯合变异:VAF=1.0
    杂合变异:VAF=0.5
  • 例子:若一个个体在某位点的基因型为 AA,则该位点无变异,VAF = 0;若基因型为 Aa,则为杂合变异,VAF = 0.5;若基因型为 aa,则为纯合变异,VAF = 1。

3. 基于群体数据的计算方法

  • 原理:在群体遗传学研究中,通常会对多个个体进行测序或基因分型,以了解某个变异在群体中的分布情况。此时,可以通过统计群体中所有个体在该变异位点上的等位基因情况来计算群体水平的 VAF。
  • 计算公式:VAF = 群体中所有个体的变异等位基因数 / 群体中所有个体的等位基因总数
  • 例子:假设在一个包含 100 个个体的群体中,某个变异位点上,所有个体的等位基因总数为 200 个,其中变异等位基因数为 40 个,则VAF = 40 / 200 = 0.2。

标签:变异,VAF,等位基因,reads,测序,位点,计算方法
From: https://blog.csdn.net/LittleComputerRobot/article/details/143672426

相关文章

  • 用Javafx开发定时器
    选中小时分钟秒代码附上:packagecom.example.javafx03;importjavafx.application.Application;importjavafx.fxml.FXMLLoader;importjavafx.scene.Parent;importjavafx.scene.Scene;importjavafx.scene.image.Image;importjavafx.stage.Stage;importjav......
  • JavaFx项目打包
    JavaFx项目打包成exe可执行程序使用Java8以上的版本编写JavaFx项目的时候,需要额外下载JavaFx并且配置环境,那么打包的时候也需要进行对应的配置。网上很多方法都打包不成功,所以这里记录一下我成功打包的过程。步骤一:在IDEA上打包成jar包打包jar包打开File>ProjectStru......
  • WIFI理论吞吐量计算方法和WIFI网卡选购
    WIFI理论吞吐量计算方法和WIFI网卡选购1计算公式 Wi-Fi理论吞吐量=(符号位长×码率×数据子载波数量×空间流数量)÷子载波传输时间2 各计算参数说明2.1符号位长一个Symbol能承载的bit数量​11a/g11N11AC11AX最大调制方式64QAM64QAM256QAM1024QA......
  • JavaFX+JavaCV实现批量视频处理及批量生成视频开发笔记--003,批量视频混剪功能设计与代
    我要使用JavaFX+JavaCV实现一个桌面应用,可以打包成Windows和Mac的桌面应用。实现的功能是:批量视频混剪。具体操作是:在界面上选择一个文件夹或多个视频文件,对文件夹中的所有视频文件(仅.mp4格式)或者选中的文件进行处理,随机截取原视频中指定长度的视频片段(如5秒),拼接成多个新的......
  • 容斥原理计算方法
    一些条件,都要满足为什么容斥问题会有一套专门的计算方法?其实容斥问题是一种常见子集答案总和的信息,常见的求解方法为DP。在求解过程中往往需要利用之前已经有了的信息,尝试整体转移,以优化时间复杂度。定义根据定义式:\(\sum_{S\subsetT}(-1)^{|S|}f_S\)进行计算,复杂度\(\ma......
  • C#联合Visionpro编程学习记录,视觉中需要考虑旋转中心工况的计算方法探讨
    一、考虑旋转中心的工况解法,1,视觉中引导定位或者对位贴合时,机械手或者xyzr轴上手爪中心和末端轴中心不同轴时,就要考虑旋转中心问题;2,如果设备的CT要求没有很苛刻,可以采用2次拍照的方案解决,1次拍照后纠偏角度,然后在纠正角度后的位置2次拍照纠正x、y偏差; 看下图:第一次拍照得到红......
  • java-快速将普通main类变为javafx类,并加载自定义fxml
    java-快速将普通main类变为javafx类,并加载自定义fxml前提步骤1.普通类继承Application2.实现main方法3.写一个controller4.写一个fxml文件5.写start方法加载fxml6.具体代码7.运行即可前提使用自带javafx的jdk,这里使用的是jdk1.834,当然你可以使用其他的可行......
  • 开源大模型占GPU显存计算方法
    运行大模型GPU占用计算公式:\(M=\frac{(P*4B)}{32/Q}*1/2\)M:以GB标识的GPU内存P:模型中的参数数量,例如一个7B模型有70亿参数4B:4个字节,表示用于每个参数的字节32:4个字节中有32位Q:应该用于加载模型的位数,例如16位、8位、4位1.2:表示在GPU内存中加载其......
  • 课前准备---单细胞数据检测SNV(变异、插入、缺失、等位基因连锁)
    作者,EvilGenius单细胞检测变异的分析已经分享了很多,全部发的高分文章。日前外显子的课程已经安排上了,但是不一定能上,可能外显子分析对大家来讲不太重要吧,但更可能是大家都会分析,不管上不上,先好好备课吧(正好也偷个懒)。不过提醒大家一句,正因为做的人少,才是机会,表达信息的分析......
  • ADC采样时间、采样周期、采样频率计算方法
    一、概述ADC转换就是输入模拟的信号量转换成数字量。读取数字量必须等转换完成后,完成一个通道的读取叫做采样周期。采样周期一般来说=转换时间+读取时间,而转换时间=采样时间+12.5个时钟周期。采样时间是你通过寄存器告诉STM32采样模拟量的时间,设置越长越精确二、STM32ADC......