I'm seeking recommendations for analyses & dealing with homopolymers? Has anyone used a software script for efficient analyses and generation of "methylation plots" directly from IonTorrent reads?
I was initially writing it in R and did the simulation work in R to demonstrate the principal, but then I swapped to Python. I'm using the Python pysam package which is a wrapper for samtools. In this instance it's much more efficient and flexible than R.