把read比對到基因組之后,需要提取唯一比對來進行下一步分析逻炊。 bowtie2和HISAT2 都沒有只輸出唯一比對的命令互亮,所以需要對比對結(jié)果sam文件進行提取,可以通過sam...
把read比對到基因組之后,需要提取唯一比對來進行下一步分析逻炊。 bowtie2和HISAT2 都沒有只輸出唯一比對的命令互亮,所以需要對比對結(jié)果sam文件進行提取,可以通過sam...
以下只是我自己學習LMM所做的筆記余素,并不保證完全正確豹休。 基于R的線性混合效應模型(linear mixed effects models)分析 why 為什么要采用線性混合效...
R包vegan的Mantel tests Mantel tests是確定兩組距離測度矩陣(而非兩組變量矩陣)之間相關(guān)性的相關(guān)性測試方法,用于判斷一個矩陣中的樣本距離與另一矩陣...
歡迎關(guān)注"生信修煉手冊"桨吊! GATK4 對于體細胞突變和生殖細胞突變的檢測分別給出了對應的pipeline: Germline SNPs+Indels Somatic SNV...
最近工作中遇到需要比對兩個不同的VCF文件,看看兩個VCF文件中哪些變異位點是相同的屏积,哪些是不同的医窿。為了尋找解決辦法,于是求救Google炊林,發(fā)現(xiàn)了一篇非常不錯的教程姥卢。于是今天...
1.箱線圖 Box plots 1.1 語法 geom_boxplot( mapping = NULL, data = NULL, stat = "boxplot", pos...
昨天看了gatk的官網(wǎng),從2018年發(fā)布正式版的4.0.0開始渣聚,到現(xiàn)在已經(jīng)更新到4.1.8独榴,在速度和準確度上都有了大幅的提升。gatk4除了整合picard軟件之外奕枝,在使用上...
Cold signaling in plants: Insights into mechanisms and regulation ##Xiao yu Guo ; Dongf...
背景篇 MITE屬于II類非自主轉(zhuǎn)座因子棺榔,并且在真核生物中存在大量的拷貝。 MITE長度大約是500bp隘道,并且擴增速度非持⑿快。2013年一篇發(fā)表在NAR的文章"P-MITE:...
BSA雖然不是我最早接觸的高通量數(shù)據(jù)分析項目(最早的是RNA-seq)谭梗,但是卻是我最早獨立開展的數(shù)據(jù)分析項目忘晤, 我曾經(jīng)專門寫過一篇文章介紹如何使用SHOREMap做擬南芥的E...
deepTools 是一套基于python開發(fā)的工具闰蛔,適用于有效處理分析高通量測序數(shù)據(jù),可用于ChIP-seq, RNA-seq 或 MNase-seq图柏。 #1. deepT...
根據(jù)相似性原理,序列相似爆办,功能相似难咕,所有功能注釋無非是用比對工具將輸入序列比對到數(shù)據(jù)庫序列课梳,再將輸入ID對應數(shù)據(jù)庫ID距辆,進一步對應到功能條目的關(guān)系余佃。 數(shù)據(jù)庫要么建到本地,要么...
寫在前面 能干什么背犯? 學習任何一個工具或是文章中的一類炫酷圖表時坏瘩,我覺得首先需要能明白:這個工具/這類圖表能干什么?這個工具/圖的展示方式是最好的么漠魏,還有哪些不足之處倔矾。然后再...
BLAST 是由美國國立生物技術(shù)信息中心(NCBI)開發(fā)的一個基于序列相似性的數(shù)據(jù)庫搜索程序。BLAST是“局部相似性基本查詢工具”(Basic Local Alignmen...
GOplot 包通過封裝好的函數(shù)可視化基因功能分析結(jié)果柱锹。 #1. 安裝 #2. GOplot 內(nèi)置數(shù)據(jù) ##2.1 腦和心臟內(nèi)皮細胞的轉(zhuǎn)錄組數(shù)據(jù) 來自文章Nolan et a...