I've been using the -s option of samtools "pileup" to output the mapping qualities for each read at each position as a last column in the pileup...