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
网友评论