美文网首页
用samtools+bedtools快速计算不同区域中的变异数量

用samtools+bedtools快速计算不同区域中的变异数量

作者: 大海龟啦啦啦 | 来源:发表于2022-07-08 10:53 被阅读0次

在做生物信息分析时,我们经常会用到滑动窗口统计不同区间内的变异数量,但是对于生物信息初学者(不具备编程能力)存在一定难度,这里有简单便捷实现滑动窗口的方法。

  1. 通过samtools获取基因组中个条染色体的长度
samtools faidx Brapa.fa
  1. 通过bedtools生成区域文件
#生成500kb的区域
bedtools makewindows -g Brapa.fai -w 500000 > region.bed
  1. 输入vcf文件,用bedtools计算每个区域存在的变异个数
bedtools coverage -a region.bed -b aaa.vcf -counts > aaa.vcf_500kb_window.txt

相关文章

  • 用samtools+bedtools快速计算不同区域中的变异数量

    在做生物信息分析时,我们经常会用到滑动窗口统计不同区间内的变异数量,但是对于生物信息初学者(不具备编程能力)存在一...

  • excel中countif函数的使用

    countif函数可以计算某个区域中满足给定条件的单元格数量。 countif(range,criteria) r...

  • 数据透视表高级应用  12月6日

    数据透视表高级应用 12月6日 一、值显示方式 自动计算占比和排名 ①数量占比 用数据透视表快速制作 把数量和金额...

  • 预测变异的效果

    基因组上有很多变异,不同变异的效果不同。除了wet-lab,现在也有很多软件预测。接下来,我打算用SIFT4G、P...

  • 什么叫遗传?疾病是遗传的吗?

    遗传学(Genetics)——自然科学领域中探究生物遗传和变异规律的科学。研究生物的遗传与变异的科学,研究基因的结...

  • 听王老师用7、8的口诀求商的感想

    周二,我们在录播教室听了王老师用7、8的口诀求商。用乘法口诀求商是数学计算中一块重要的基石,它在整个计算领域中起...

  • 重测序分析(8)SV CNV 注释

    结构变异 结构变异(Structural Variation,SV)这种类型比较多,根据结构变异的不同类型,结构变...

  • 类文件结构

    一、概述 在计算机商业领域中,不同的硬件体系结构和不同的操作系统长期并存是必然的结果,所以Java 提出了“Wri...

  • 看板数量的计算

    看板数量的计算 看板是丰田生产方式实现准时化的工具,从手段-目的链的角度来说,准时化是目的,看板是手段。看板用来实...

  • Excel笔记-公式篇countif

    1.计算指定区域中满足指定条件的单元格个数 =COUNTIF(C3:C9,H3) 2.统计某一范围内的单元格数量 ...

网友评论

      本文标题:用samtools+bedtools快速计算不同区域中的变异数量

      本文链接:https://www.haomeiwen.com/subject/ymambrtx.html