240 发简信
IP属地:北京
  • 120
    DIY 地热系统

    许多人问,为啥地热系统这么贵,实际上我们不必一定搞贵的,接下来我告诉你如何自制地热系统,而效果如同专业安装的。 你需要的技能:1. 水管安装 2. 挖掘能力(挖坑) 地热系统...

  • P(Genotype|Data)是如何计算的?

    HaplotypeCaller是怎么计算PL的?(VCF解密)

    写在前面 PL是HaplotypeCaller等GATK变异检测软件在sample层面给出的注释,一般记录在VCF文件的FORMAT/sample一栏中。它代表了根据某位点变...

  • htslib 的使用( C++ 处理BAM文件)

    背景 在处理NGS 数据时,尤其是突变的检测。我们需要对bam 文件进行读写。 目前有相对方便的库 pysam 和 htsjdk ,这两个库提供的api 可以方便python...

  • 但是,貌似qualimap 对 比较大的bam, 速度奇慢。

    比对文件BAM/CAM深度检测好用工具:Mosdepth

    做完比对后,最直接的一个问题就是想知道:你的reads比对上了哪些区域没比对上哪些区域,有多少reads比对上特定的区域(该区域比对深度)。传统手艺一般可以使用bedtool...

  • -subject_besthit 此参数确定你能输出 最有的一条吗? 因为之前的blast结果,总是有bug, 总是不能输出最优的一条, 我在github中见到过此issue

    linux BLAST序列比对 (nt/nr库)

    安装本地blast序列比对软件,我们可以搜索一个查询序列定制数据库,例如想研究一个新测序的基因组,或者感兴趣的一组蛋白质序列。有时我们希望把程序插入到一个流程中,例如搜索一个...

  • P(Data | Genotype) 写反了吧, 官网是P(Genotype|Data)

    HaplotypeCaller是怎么计算PL的?(VCF解密)

    写在前面 PL是HaplotypeCaller等GATK变异检测软件在sample层面给出的注释,一般记录在VCF文件的FORMAT/sample一栏中。它代表了根据某位点变...