use Nucleotide mismatch frequency from bam files as a score to visualize

[al2na]

PAR-CLIP or BS-seq relies on nucleotide mismatch frequencies as genome-wide scores. We can use pileup function from samtools to extract those frequencies to a scoreMatrix object???