Plot coverage per base for BAM file
From your bam-file, you can generate a BED(-graph) of per base coverage with bedtools genomecov (use the -ibam, -d and if you want a BEDgraph also the -bg flags).

There are multiple R packages for bioconductor which read and visualise these, including sushi.