本篇內容主要講解“如何使用Beagle進行基因型填充”,感興趣的朋友不妨來看看。本文介紹的方法操作簡單快捷,實用性強。下面就讓小編來帶大家學習“如何使用Beagle進行基因型填充”吧!
創(chuàng)新互聯(lián)建站專業(yè)為企業(yè)提供賓川網站建設、賓川做網站、賓川網站設計、賓川網站制作等企業(yè)網站建設、網頁設計與制作、賓川企業(yè)網站模板建站服務,十載賓川做網站經驗,不只是建網站,更提供有價值的思路和整體網絡服務。
軟件的官網如下
https://faculty.washington.edu/browning/beagle/beagle.html
相比同類軟件,該軟件有兩個顯著的優(yōu)勢,第一個就是內存消耗小,第二個運行速度快。對于數(shù)據(jù)量如此大的reference panel數(shù)據(jù),提出了一種新的數(shù)據(jù)存儲格式bref3
。文章中比較了不同格式消耗的內存大小,結果如下
可以看到,bref3格式的內存消耗最小。關于運行時間的比較結果如下
可以看到Beagle 5.0版本的運行時間最短。該軟件采用java語言進行開發(fā),安裝簡單,直接下載jar文件即可?;居梅ㄈ缦?/p>
java \
-jar beagle.24Aug19.3e8.jar \
ref=chr1.1kg.phase3.v5a.b37.bref3 \
map=plink.GRCh48.map \
gt=input.chr1.vcf.gz \
out=out.chr1.gt
ref
參數(shù)指定reference panel的數(shù)據(jù),支持vcf和bref3兩種格式,官網提供了1000G的數(shù)據(jù)供下載,鏈接如下
http://bochet.gcc.biostat.washington.edu/beagle/1000_Genomes_phase3_v5a/
map
參數(shù)指定參考基因組的連鎖圖譜,是一個可選參數(shù),官網提供了human對應的連鎖圖譜,鏈接如下
http://bochet.gcc.biostat.washington.edu/beagle/genetic_maps/
到此,相信大家對“如何使用Beagle進行基因型填充”有了更深的了解,不妨來實際操作一番吧!這里是創(chuàng)新互聯(lián)網站,更多相關內容可以進入相關頻道進行查詢,關注我們,繼續(xù)學習!