We built the sciClone package for exactly this purpose: https://github.com/genome... It takes inputs of somatic mutations, with readcounts and VAFs, and uses that information to infer subclonal populations in heterogeneous tumors. It also gives you some nice visualization options. - Chris Miller