Question: Sort/Index bam prior to peak-calling ChIP-seq & deeptools bigwig visualisation
gravatar for andrew.phipps.00
2.7 years ago by
andrew.phipps.000 wrote:

Hi biostars community,

First off I am quite new to NGS data analysis, and I am trying to find out answers for a few questions.

I have gotten ChIP-seq data from a time-course and disease experiment, that I have trimmed and aligned with bowtie2, giving me .sam files, which I am converting to .bam with SAMtools. I want to check the read coverage for QC (with MultiBamSummary) prior to peak-calling with MACS2 (and would like to visualise with IGV), but I have run into a few options and am not sure which is the right way to go:

Should I use SAMtools to sort and index the bam files prior to peak calling? (and can visualise the sorted.bam & index files in IGV) - is it necessary to sort .bam files prior to peak calling?

Or is it better to use deepTools' Bam coverage to generate bigwig files from the .bam files for visualisation with IGV, and then peak call with the non-sorted/indexed .bam files - saving some processing time?

Thanks for the help and suggestions, Andrew

Edit: I think I answered my own question here - deepTools multibamsummary requires indexed files to will sort and index the files!

chip-seq alignment • 2.2k views
ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by andrew.phipps.000

Coordinate-sorting your alignment is always a good choice as it is necessary for indexing and most downstream tools. You can avoid the sam file by simply piping the alignment directly to SAMtools sort:

bowtie2 (...parameters...) | samtools sort -o sorted.bam -
ADD REPLYlink written 2.7 years ago by ATpoint40k

Thanks for the suggestion! I hadn't thought of that (would have saved me a heap of time...) - I am going to try a few different alignment and peak calling tools, so will keep this in mind for the future! Cheers!

ADD REPLYlink written 2.7 years ago by andrew.phipps.000

You're welcome. There are plenty of shortcuts and little workarounds that you'll learn about once you get more experienced ;-) I strongly encourage ypu to play around with tools and commands if you have a bit of time here and there.

ADD REPLYlink written 2.7 years ago by ATpoint40k

how to index the sort.bam files

ADD REPLYlink written 2.1 years ago by lzhizhao91790

Please use the search function and google for this. The first hit would have redirected you to this post, instructing you to do samtools index in.bam ;-)

ADD REPLYlink written 2.1 years ago by ATpoint40k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1721 users visited in the last hour