是公众号回格式,朋友们,这里回没有用
如果想做植物的转录因子调控网络,我们可能会有以下问题。 1:如何利用有限的叶片原生质体材料获得高质量的TF ChIP-seq数据?2:如何生成可信的无标度转录调控网络图?3:...
ncbi的BLAST比我想象中的要强大,用普通blastn貌似不能将basepire数小于50的序列比对到基因组。像我,想将一些15bp的序列比对到基因组上做统计。其实BLA...
写在前面 两个月前,我推送了一个师妹的投稿,介绍了orthofinder软件的应用场景和具体安装,属于上篇。有上便有下,评论区和后台看到不少朋友在催更下篇。今天即推出,主要内...
大佬您好,方便能看看您生成的raw.matrix的文件格式嘛?我用HiC-pro产出的raw.matrix用您写的脚本计算分辨率:
···
for f in *raw.matrix; do
# echo $f
high_bin=$(cat $f | awk '{x[$1] += $3}END{for (i in x){print i, x[i]}}' | sort -k 1,1n | awk '($2>1000){sum=sum+1}END{print sum}');
total_bin=$(cat ${f/raw.matrix/raw_abs.bed} | wc -l);
awk -v high_bin=$high_bin -v total_bin=$total_bin -v name=$f 'BEGIN{print name,high_bin/total_bin}' >> ./hic_res_suitable.txt ;
done
···
但是得出的结果很怪,越大的分辨率反而,>1000 contact 的 bin比例越少。
···
sample_1000000.matrix 0.819144
sample_100000.matrix 0.958256
sample_10000.matrix 0.905145
sample_15000.matrix 0.93689
sample_20000.matrix 0.951743
sample_250000.matrix 0.915968
sample_25000.matrix 0.960175
sample_30000.matrix 0.965011
sample_500000.matrix 0.869427
sample_50000.matrix 0.970463
···
【HiC】数据如何选择合适分辨率?参考: 比如原位HiC,DLO-HiC 实验流程不同,reads利用率存在差别。但是你的数据量适合多大分辨率,是否满足TAD(40Kb),Loop(10Kb/5Kb)呢? 上...
这个题目是转载别人的,原出处没能找到,原作者看到可以联系删除。 当时在做这个转录组分析之前不太了解这个到底是什么,但是跟着流程完成了,也得到了结果,但是还是有很多问题不太明白...
记录下下载过程,为自己和后人避坑。 1.Conda连接不上镜像源问题 首先是anaconda安装软件或创建环境时遇到的问题。即使换完清华源和其他镜像源以后依旧报错。 尝试了很...
1. 什么是单倍型? 同源染色体:同源染色体,一个来自母本,一个来自于父本。 单倍型:单倍体基因型的简称。遗传学上指在单条染色体上一系列遗传变异位点的组合。 2. 单倍型组装...