工具鏈安裝 LACHESIS依賴boost庫(kù),且經(jīng)過測(cè)試1.52~1.67版本的boost庫(kù)可以順利編譯LACHESIS,版本過高會(huì)導(dǎo)致報(bào)錯(cuò),而低版本的boost庫(kù)在高版本的...
創(chuàng)建環(huán)境 編輯tremolo.yml文件,內(nèi)容如下: 創(chuàng)建conda環(huán)境 啟動(dòng)環(huán)境 安裝RaGOO 安裝TrEMOLO 修改環(huán)境 安裝samtools 修改TrEMOLO 測(cè)...
在如何用WGDI進(jìn)行共線性分析(上)[http://www.itdecent.cn/p/a50548e81ac0], 我們基于blastp的結(jié)果繪制了點(diǎn)圖,之后用 -icl...
0.部分參考https://doi.org/10.1093/molbev/msae239[https://doi.org/10.1093/molbev/msae239] 用了...
1. 群體演化歷史分析的案例: 群體演化包括有效群體大小(Ne)變化、基因流,遷徙,分化等,會(huì)對(duì)等位基因頻率產(chǎn)生顯著影響,塑造了現(xiàn)有遺傳多樣性的模式和水平。群體演化、遺傳漂變...
偶然間發(fā)現(xiàn)了下面這個(gè)圖的代碼和數(shù)據(jù) 這個(gè)圖的代碼和數(shù)據(jù)的鏈接 https://github.com/nrennie/tidytuesday/blob/main/2021/07...
1. Summary 基因家族擴(kuò)張收縮分析可細(xì)分為六部分: 獲取每個(gè)基因?qū)?yīng)的最長(zhǎng)編碼區(qū)轉(zhuǎn)錄本 OrthoFinder聚類基因家族 系統(tǒng)發(fā)育樹推斷 物種分歧時(shí)間推斷 基因家族...
前言:在對(duì)miRNA進(jìn)行靶標(biāo)預(yù)測(cè)分析時(shí),需要特定物種的轉(zhuǎn)錄本作為靶向的數(shù)據(jù)庫(kù)。通常對(duì)miRNA預(yù)測(cè)以3'UTR區(qū)域?yàn)橹鳎@就需要對(duì)轉(zhuǎn)錄本數(shù)據(jù)庫(kù)進(jìn)行UTR區(qū)域的提取。我在之前寫...
目的:從測(cè)序數(shù)據(jù)(clean reads) 到樣本的GVCF文件。利用bwa,samtools,和GATK 獲得每個(gè)樣本的GVCF文件。
mcmctree需要3個(gè)輸入文件 多序列比對(duì)的phy格式文件(注意:這個(gè)phy和一般軟件的phy不太一樣)使用腳本把多序列比對(duì)的fa格式轉(zhuǎn)為phy格式 準(zhǔn)備包含2個(gè)物種之間的...
本次推送是文獻(xiàn)分享22的對(duì)應(yīng)內(nèi)容。我與生信,公眾號(hào):我與生信文獻(xiàn)分享22:泛基因組解析柑橘亞科進(jìn)化以及柑橘果實(shí)中檸檬酸積累的關(guān)鍵基因[http://mp.weixin.qq....
當(dāng)我們進(jìn)行群體分析時(shí),獲得vcf文件后,可以根據(jù)變異位點(diǎn)對(duì)這些樣本進(jìn)行PCA分析,現(xiàn)簡(jiǎn)單介紹 1、軟件安裝 2、簡(jiǎn)單操作 本次使用train.vcf.gz 作為例子 第一步:...