符合過濾條件的snp應(yīng)該符合以下參數(shù)：
QD > 2.0
FS < 60.0
MQ > 40.0
MQRankSum > -12.5
ReadPosRankSum > -8.0
SOR < 3.0
而不是FS > 60.0,SOR > 3.0。對應(yīng)硬過濾中視為需要過濾的位點的表達(dá)式為"QD < 2.0 || MQ < 40.0 || FS > 60.0 || SOR > 3.0 || MQRankSum < -12.5 || ReadPosRankSum < -8.0"。沒考慮過為什么你在table的時候結(jié)果剛好相反嗎？
「簡化基因組」如何過濾用GATK分析得到的SNPGATK官方提供了一個SNP過濾的標(biāo)準(zhǔn)，howto-apply-hard-filters-to-a-call-set,如果你按照它的要求來過濾簡化基因組中的SNP數(shù)據(jù)，也就是...
大家好，本周給大家分享的是2021年12月份發(fā)表在Current biology上與中國家馬體型相關(guān)主效基因突變的一篇文章。 文章題目：A single-nucleotide...
pixy用于從vcf變異文件中計算pi（π）,Dxy和Fst的值。 1.使用conda安裝pixy github[https://github.com/ksamuk/pixy...
終于超過1k粉絲啦～之前說要小小慶祝一下，害怕直接放文章里會被封文章，所以福利大家可以見評論。前幾周忙著答辯畢業(yè)事宜沒有時間更新，雖然這幾天也挺忙哈哈哈哈，不過還是要定期更新...
許多R包都可以畫桑基圖和輪狀圖，如 circlize, networkD3, and ggalluvial。這些包都非常好，非常有用。問題來了，吃飽了撐的再弄一個同樣的畫圖包...
aplot(https://github.com/YuLab-SMU/aplot[https://github.com/YuLab-SMU/aplot])是Y叔開發(fā)的一個基于...
MOVICS包整合多組學(xué)數(shù)據(jù)進(jìn)行分析，其安裝的過程太讓人心酸，大家安裝的過程慢慢體會，下面就來學(xué)習(xí)下這個由“大魚海棠”老師開發(fā)的R包吧。 1.示例數(shù)據(jù) brca.tcga.R...
論文是 Global burden of lung cancer attributable to ambient fine particulate matter pollut...
每次有正事的時候就想逃避寫個簡書，假裝自己沒有干不正經(jīng)的事。（明明我還有一些數(shù)據(jù)需要分析呢，但是毫無分析的動力。）好的，今天我們來介紹一個很酷炫的圖：弧形圖（Arc Diag...
@斬毛毛 請問如何隨機(jī)選取20ksnp呢？
admixture 群體結(jié)構(gòu)分析tructure是與PCA、進(jìn)化樹相似的方法，就是利用分子標(biāo)記的基因型信息對一組樣本進(jìn)行分類，分子標(biāo)記可以是SNP、indel、SSR。相比于PCA，進(jìn)化樹，群體結(jié)構(gòu)分析可明...
請問小鼠的全基因組重測序需要基因型填充嗎？
群體遺傳學(xué)習(xí)筆記-基因型缺失數(shù)據(jù)的填充Genotype Imputation是在高通量測序中常出現(xiàn)的定義，按照義譯就是基因型填充。要真正理解imputation這個概念，我們就需要先理解基因型缺失（genotyp...
請問圖一熱圖的配色是什么呢？感覺很好看
10X單細(xì)胞WGCNA分析聯(lián)合10X空間轉(zhuǎn)錄組WGCNA分析各位同學(xué)，大家好，今天我們繼續(xù)WGCNA的分析點，這一次，我們需要聯(lián)合兩種技術(shù)，10X單細(xì)胞和10X空間轉(zhuǎn)錄組，也就是在多組學(xué)層面上分析WGCNA得到的模塊，關(guān)于WGCNA，...
