This is a test version of Biostars. For the public version, visit https://www.biostars.org.
two way get analysis data to get the obvious range difference in IGV

Hi all, I do ChIP-seq for TF. And here I get tdf file through two different way. The up one i use bowtie2 and macs2 to analysis fastq file and get bedgraph file then use igvtool transfer to tdf file(the range from 0 -200). the bottom one I use whole homer analysis to get. The peak location look like same, but the peak range has big difference(the range form 0 -100). Did some know which one should be the right one ? And why this range different can happen? Thanks in adavance.

chip-seq igv peak range

Thanks! Can you see the picture now? I can not see the link by myself.

do not know how to upload the figure in my computer

Right! Thank you for your help!

1 answer

Different tools, different outputs, mainly because (probably) some default parameters cause different behaviours. I recommend the latest version of deeptools bamCoverage to have things nice and reproducible:

bamCoverage -bs 1 --normalizeUsing CPM --bam in.bam -o file_CPM.bigwig -e <int>

This will give you a high-resolution browser track at base-pair resolution (-bs 1, the lower the bin size the smoother the peaks will look, but the program runs longer and files are bigger), normalized to the count per million (CPM). You can optionally and an integer value to -e to elongate the short reads to the average fragment size (if you have that information, based on gel or Bioanalyzer).

Log in to answer this question.