1. 序列比对到图基因组 vg map -m short -t 24 -x graph.xg -g graph.gcsa -f Sample40.1.fq.gz -f Sam...
1. 序列比对到图基因组 vg map -m short -t 24 -x graph.xg -g graph.gcsa -f Sample40.1.fq.gz -f Sam...
下载文件夹: /public/home/kongweilong/Biosofts/obsutil/obsutil share-cp https://e-share.obs-w...
https://www.jianshu.com/p/d1d8d59f61fe
确定阈值分为了不同的流派。 第一种方法,Bonferroni矫正 这种方法应用最多,一般是N表示参与分析的SNP的个数,阈值可以是: * 1/N * 0.05/N * 0.0...
1. 代谢组数据示例(test_meta.txt) SampleID mws1401 A13 0.00 A16 0.00 A7 0.00 A17 0.00 A11 0.00 ...
gatk --java-options "-Xmx185g -Djava.io.tmpdir=./tmp" HaplotypeCaller -R /public/home/k...
建立一个.inputrc文件夹 # /etc/inputrc - global inputrc for libreadline # See readline(3readlin...
1. 原因:利用gmaps回帖基因,会损失和错误很多,导致busco结果降低。所以可以直接基于挂载生成的agp文件直接进行转化,不损失任何基因。 #!/usr/bin/env...
Mark一个很经典的三维基因组学学习笔记 参考:https://geekdaxue.co/read/davey-brweo@xfltrm/ncf5i6
1. Hic数据mapping(Align Hi-C data to the assembly, remove PCR duplicates and filter out s...
plink --vcf All.input.vcf --biallelic-only --recode vcf-iid --out input.clean --allow-e...
问题1.因为WGBS测序的原理是亚硫酸盐处理会使未甲基化的胞嘧啶(C)被转化为脱氨基胞嘧啶(T),而甲基化的胞嘧啶(5-methylcytosine,5mC)不受影响。然后通...
1.Bam2fq (by 森哥) 脚本如下: #!/usr/bin/perl =head1 Description ccsbam2fq.pl read BAM file of...
在家目录上新建一个.inputrc文件,内容为: # do not bell on tab-completion #set bell-style none set meta-...
SIFT4G软件默认假设参考基因组上的氨基酸序列是功能正常的,因此它的分析主要关注突变位点上的氨基酸替代对蛋白质功能的可能影响。 SIFT4G会对给定的突变位点进行分析,计算...
利用Canu进行组装 canu -p YueC -d canu genomeSize=3g useGrid=false maxThreads=24 -pacbio-hifi ...
如果两个 VCF 文件中的样本出现的顺序不一样,但具有相同的样本名称,您可以使用 bcftools 的 reheader 命令来调整样本顺序。 bcftools rehead...
1.群体遗传|根据SNP验证自然选择—XPEHH分析 (https://mp.weixin.qq.com/s/1IM6KSCGIq0cUJhKKXVcIw) 2.Human ...
import sys input_file = sys.argv[1] step_len = 100000 result = {} chrom_stats = { "Ch...
Pmap 输入数据文件 col 设置不同染色体中点的颜色 cex 设置点的大小 pch 设置点的形状 band 设置不同染色体之间的间隔 H 设置每个圈的高度...