Hello, you can try something like this :

samtools view -F 4 file.bam | awk '{printf $1"\n"}' | sort | uniq -cd

Get the reads mapped, take the read names, sort them, find and count the non unique ones. I let you try to plot your own histogram with these infos.

took it from here