2019-07-05 macs14 bdg to bw

samtools merge -@ 45 final.bam 1in.bam 2in.final.bam

macs14 -t 2cell-rep1.bam -g hs -n 2cell-rep1 -B --nolambda --nomodel


###bdg to bw

cd /treat/

gunzip *

cat * > all.bdg

bedtools intersect -a all.bdg -b /media/shen/disk1/jjc/reference/hg19/nochr.hg19.chr.limited.bed -wa -f 1.00 > all.bdg.tmp

sort -k1,1 -k2,2n all.bdg.tmp -o sort.all.bdg.tmp

bedGraphToBigWig sort.all.bdg.tmp /media/shen/disk1/jjc/reference/hg19/nochr.hg19.chrom.sizes all.bw

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容