Bedtools coverage histogram



I am following this tutorial in order to plot che coverage depth for my target exomes: Visualize coverage for targeted NGS (exome) experiments. I also want to produce a plot similar to the post.

So I am using the following command line for all my bam files:

bedtools coverage -hist -abam {} -b input.bed | grep ^all > {}.hist.all.txt

I get an output where depth of coverage is either 0 or 1. But in the plot from the post the values are from 0 to 400. If i change bam files as -b and bed file as -a the output is empty. Why is that and what am I doing wrong?

all     0       697333613       802606707       0.5315176
all     1       105273094       802606707       0.4684824



Source link