符合过滤条件的snp应该符合以下参数：
QD > 2.0
FS < 60.0
MQ > 40.0
MQRankSum > -12.5
ReadPosRankSum > -8.0
SOR < 3.0
而不是FS > 60.0,SOR > 3.0。对应硬过滤中视为需要过滤的位点的表达式为"QD < 2.0 || MQ < 40.0 || FS > 60.0 || SOR > 3.0 || MQRankSum < -12.5 || ReadPosRankSum < -8.0"。没考虑过为什么你在table的时候结果刚好相反吗？
「简化基因组」如何过滤用GATK分析得到的SNPGATK官方提供了一个SNP过滤的标准，howto-apply-hard-filters-to-a-call-set,如果你按照它的要求来过滤简化基因组中的SNP数据，也就是...
大家好，本周给大家分享的是2021年12月份发表在Current biology上与中国家马体型相关主效基因突变的一篇文章。 文章题目：A single-nucleotide...
pixy用于从vcf变异文件中计算pi（π）,Dxy和Fst的值。 1.使用conda安装pixy github[https://github.com/ksamuk/pixy...
终于超过1k粉丝啦～之前说要小小庆祝一下，害怕直接放文章里会被封文章，所以福利大家可以见评论。前几周忙着答辩毕业事宜没有时间更新，虽然这几天也挺忙哈哈哈哈，不过还是要定期更新...
许多R包都可以画桑基图和轮状图，如 circlize, networkD3, and ggalluvial。这些包都非常好，非常有用。问题来了，吃饱了撑的再弄一个同样的画图包...
aplot(https://github.com/YuLab-SMU/aplot[https://github.com/YuLab-SMU/aplot])是Y叔开发的一个基于...
MOVICS包整合多组学数据进行分析，其安装的过程太让人心酸，大家安装的过程慢慢体会，下面就来学习下这个由“大鱼海棠”老师开发的R包吧。 1.示例数据 brca.tcga.R...
论文是 Global burden of lung cancer attributable to ambient fine particulate matter pollut...
每次有正事的时候就想逃避写个简书，假装自己没有干不正经的事。（明明我还有一些数据需要分析呢，但是毫无分析的动力。）好的，今天我们来介绍一个很酷炫的图：弧形图（Arc Diag...
@斩毛毛 请问如何随机选取20ksnp呢？
admixture 群体结构分析tructure是与PCA、进化树相似的方法，就是利用分子标记的基因型信息对一组样本进行分类，分子标记可以是SNP、indel、SSR。相比于PCA，进化树，群体结构分析可明...
请问小鼠的全基因组重测序需要基因型填充吗？
群体遗传学习笔记-基因型缺失数据的填充Genotype Imputation是在高通量测序中常出现的定义，按照义译就是基因型填充。要真正理解imputation这个概念，我们就需要先理解基因型缺失（genotyp...
请问图一热图的配色是什么呢？感觉很好看
10X单细胞WGCNA分析联合10X空间转录组WGCNA分析各位同学，大家好，今天我们继续WGCNA的分析点，这一次，我们需要联合两种技术，10X单细胞和10X空间转录组，也就是在多组学层面上分析WGCNA得到的模块，关于WGCNA，...
