把read比对到基因组之后,需要提取唯一比对来进行下一步分析。 bowtie2和HISAT2 都没有只输出唯一比对的命令,所以需要对比对结果sam文件进行提取,可以通过sam...
把read比对到基因组之后,需要提取唯一比对来进行下一步分析。 bowtie2和HISAT2 都没有只输出唯一比对的命令,所以需要对比对结果sam文件进行提取,可以通过sam...
以下只是我自己学习LMM所做的笔记,并不保证完全正确。 基于R的线性混合效应模型(linear mixed effects models)分析 why 为什么要采用线性混合效...
什么是Ka/Ks The ratio of the number of nonsynonymous substitutions per nonsynonymous site ...
R包vegan的Mantel tests Mantel tests是确定两组距离测度矩阵(而非两组变量矩阵)之间相关性的相关性测试方法,用于判断一个矩阵中的样本距离与另一矩阵...
欢迎关注"生信修炼手册"! GATK4 对于体细胞突变和生殖细胞突变的检测分别给出了对应的pipeline: Germline SNPs+Indels Somatic SNV...
最近工作中遇到需要比对两个不同的VCF文件,看看两个VCF文件中哪些变异位点是相同的,哪些是不同的。为了寻找解决办法,于是求救Google,发现了一篇非常不错的教程。于是今天...
1.箱线图 Box plots 1.1 语法 geom_boxplot( mapping = NULL, data = NULL, stat = "boxplot", pos...
昨天看了gatk的官网,从2018年发布正式版的4.0.0开始,到现在已经更新到4.1.8,在速度和准确度上都有了大幅的提升。gatk4除了整合picard软件之外,在使用上...
Cold signaling in plants: Insights into mechanisms and regulation ##Xiao yu Guo ; Dongf...