average coverage in target regions

[cmccabe]

I am trying to get the _X coverage in a bam file for target regions in a bed.

find *bam | parallel 'bedtools coverage -hist -abam IonXpress_009_150603.bam -b epilepsy70_medex.bed | grep ^all > output.bam.hist.txt' > all.txt

The files are too large to attach but after the -abam is the bam file to use, after the -b is the terget bed to use with the baits, after the grep^all is the bedTools coverage output and the > all.txt are the results. Is there something missing? Thank you .

EDIT:
Maybe this:

coverageBed -abam IonXpress_009_150603.bam -b epilepsy70_medex.bed -counts > count.txt