GATK的HaplotypeCaller 應(yīng)該是目前最常用的變異檢測(cè)軟件邻眷，尤其是在人類(lèi)基因組上。不過(guò)HaplotypeCaller的速度相對(duì)于其他軟件唠雕，例如bcftools,...
發(fā)簡(jiǎn)信
IP屬地：陜西
-
「GATK 4」如何提高HaplotyperCaller的效率 -
差異剪接分析軟件rMATS結(jié)果文件解讀(v4.1.2)第一次做差異剪接分析港谊，rMATS做完后一頭霧水骇吭，查閱資料整理了一下結(jié)果文件。 我使用的版本是v4.1.2歧寺，用conda安裝的燥狰，軟件的安裝和使用就不細(xì)說(shuō)了棘脐，網(wǎng)上已經(jīng)有很多帖子，...
-
感謝作者的分享
生成geno2go那一步有點(diǎn)慢龙致，用這個(gè)來(lái)生成geno2go的表會(huì)快很多
gos_list <- function(x){
the_gos <- str_split(x[2], ",", simplify = FALSE)[[1]]
df_temp <- data.frame(GID = rep(x[1], length(the_gos)),
GO = the_gos,
EVIDENCE = rep("IEA", length(the_gos)))
return(df_temp)
}
gene2gol <- apply(as.matrix(gos),1,gos_list)
gene2gol_df <- do.call(rbind.data.frame, gene2gol)熱愛(ài)大自然的小和尚 評(píng)論自生信 | 構(gòu)建物種的OrgDb -
生信基因功能分析工具：Orthofinder使用教程
最近Orthofinder2開(kāi)始陸續(xù)更新荆残，文章已經(jīng)投發(fā)到biorkix上，在網(wǎng)上搜了一圈净当，關(guān)于Orthofinder的使用的中文教程幾乎是空白的内斯，這周就借此機(jī)會(huì)和大家一起來(lái)學(xué)...
-
你好，想問(wèn)下你的問(wèn)題解決了嗎像啼，我做vcftools和plink都出現(xiàn)了跟你差不多的問(wèn)題俘闯，但兩個(gè)出現(xiàn)的問(wèn)題又不完全一樣，首先生成的ped文件應(yīng)該都沒(méi)問(wèn)題忽冻，但vcftools生成的map文件第一列為0真朗，也是跟你一樣的錯(cuò)誤沒(méi)法識(shí)別染色體，但第二列顯示正常的snp位置僧诚，而plink生成的map文件則是第一列正常遮婶，第2列為0。這樣生成的文件沒(méi)法用eigensoft進(jìn)行pca分析湖笨，而用plink卻可以獨(dú)立完成pca分析旗扑，不知這個(gè)問(wèn)題是否跟你一樣？
2021-03-17 在linux上將vcf文件轉(zhuǎn)plink的格式bed,bim,fam
我現(xiàn)在的問(wèn)題是這樣的vcf文件轉(zhuǎn)plink的格式方法一vcftools 出錯(cuò)這樣 方法二plink 都是一樣的結(jié)果 就無(wú)語(yǔ)之前我的文件有這些 現(xiàn)在是這些 然后看到了這些 然后...
-
