符合過濾條件的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月份發(fā)表在Current biology上與中國家馬體型相關主效基因突變的一篇文章。 文章題目：A single-nucleotide...
pixy用于從vcf變異文件中計算pi（π）,Dxy和Fst的值。 1.使用conda安裝pixy github[https://github.com/ksamuk/pixy...
終于超過1k粉絲啦～之前說要小小慶祝一下，害怕直接放文章里會被封文章，所以福利大家可以見評論。前幾周忙著答辯畢業(yè)事宜沒有時間更新，雖然這幾天也挺忙哈哈哈哈，不過還是要定期更新...
許多R包都可以畫?；鶊D和輪狀圖，如 circlize, networkD3, and ggalluvial。這些包都非常好，非常有用。問題來了，吃飽了撐的再弄一個同樣的畫圖包...
aplot(https://github.com/YuLab-SMU/aplot[https://github.com/YuLab-SMU/aplot])是Y叔開發(fā)的一個基于...
MOVICS包整合多組學數據進行分析，其安裝的過程太讓人心酸，大家安裝的過程慢慢體會，下面就來學習下這個由“大魚海棠”老師開發(fā)的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分析各位同學，大家好，今天我們繼續(xù)WGCNA的分析點，這一次，我們需要聯合兩種技術，10X單細胞和10X空間轉錄組，也就是在多組學層面上分析WGCNA得到的模塊，關于WGCNA，...
