第七步,deeptools的可视化 具体见:https://mp.weixin.qq.com/s/a4qAcKE1DoukpLVV_ybobA在ChiP-seq 讲解。 首先把bam文件转为bw文件,详情:http://www.bio-info-trainee.com/1815.html cd ~/project/atac/deeptools_result #source activate atac # 由于原本电脑存在deeptools所以就没必要激活了 #ls *.bam ...
cmd=commandArgs(trailingOnly=TRUE);input=cmd[1];output=cmd[2];a=abs(as.numeric(read.table(input)[,1]));png(file=output);hist(a,main="Insertion Size distribution",ylab="Read Count",xlab="Insert Size",xaxt="n",breaks=seq(0,max(a),by=10));axis(side=1,at=seq(0,max(a),by=100...