Tag: treat_pileup

ChIPseq gapped and braod peaks file are empty

ChIPseq gapped and braod peaks file are empty 0 Hello, I am using macs2 for peak calling , on ChIPseq data using the following command. I have done the normalization using the drosophila genome on my data. macs2 callpeak -t /Normalized_samples_ip/_IP_unique_rmdup_sampled.bam \ -c _001_unique_rmdup_sampled.bam -f BAM –gsize 3.0e9 –broad –broad-cutoff…

Continue Reading ChIPseq gapped and braod peaks file are empty

1 is not found in chromosome sizes file

1 is not found in chromosome sizes file 0 I am trying to convert bedGraph to bigWig file with bedGraphToBigWig. but I always get “1 is not found in chromosome sizes file”. I used the following script macs2 callpeak -c SRRc_sorted.bam -t SRRe_sorted.bam -B –nomodel –extsize 200 –SPMR -q 0.01…

Continue Reading 1 is not found in chromosome sizes file

Calculate scores using certain method by comparing a bedGraph file from treatment and a file from control representing local bias

bdgcmp(1) Calculate scores using certain method by comparing a bedGraph file from treatment and a file from control representing local bias SYNOPSIS bdgcmp <-t TREATMENT.BEDGRAPH> <-c CONTROL.BEDGRAPH> <-o OUTPUT.BEDGRAPH> [-m METHOD] DESCRIPTION Calculate scores using certain method by comparing a bedGraph file from treatment and a file from control representing…

Continue Reading Calculate scores using certain method by comparing a bedGraph file from treatment and a file from control representing local bias

difference between treat_pileup and bdgcmp fold enrichment tracks macs2

difference between treat_pileup and bdgcmp fold enrichment tracks macs2 0 Hello, I created bigwig file from a treat_pileup.bdg file generated by macs2 and also used treat_pileup.bdg and control_lambda.bdg with macs2 bdgcmp. Here is my codes; macs2 callpeak -t sample.bam -c sample_input.bam -g hs -f BAM -q 0.001 –bdg –outdir /folder…

Continue Reading difference between treat_pileup and bdgcmp fold enrichment tracks macs2