How to get a coverage graph in WIG file format directly from an alignment (BAM)

The following command combines the samtools mpileup tool with command line perl to obtain genomic coverage plot from an alignment in BAM format.

samtools mpileup -BQ0 run.sorted.bam | perl -pe '($c, $start, undef, $depth) = split;if ($c ne $lastC || $start != $lastStart+1) {print "fixedStep chrom=$c start=$start step=1 span=1\n";}$_ = $depth."\n";($lastC, $lastStart) = ($c, $start);' | gzip -c > run.wig.gz

Receive updates about NGS articles and trainings

Share this article

Last updated on August 09, 2013

ecSeq is a bioinformatics solution provider with solid expertise in the analysis of high-throughput sequencing data. We organize public workshops and conduct on-site trainings on NGS data analysis.

Would you like to receive updates about our NGS trainings and solutions? Then sign-up for our newsletter