I am currently trying to analyze whole genome sequencing data in lung cancer and have split the samples (80 in total) into 2 groups - those with doubled genomes and those without. I would like to see what genes are disproportionately mutated in the doubled genome samples compared to the non-doubled genome samples (and vice versa). I have the mutations for each sample in MAF format. Is there a package that can easily perform this type of analysis? If not, how should I approach this?