This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Somatic mutations per gene

Hi,

I have a vcf file (from Mutect2) - list of somatic mutations and their genomic locations & other info.

Does anyone know that there is a R package that can calculate the number of somatic mutations per gene?

For example, map the somatic mutations in the vcf file to any one gene based on their location and count how many mutations in each gene?

I can do it myself manually, but if there's a package that can do that, I can save time.

Thanks..

r gene mutect2 vcf somatic mutation
• 2,086 views • link
updated by Alewa Alewa • written by zephyr_falcon zephyr_falcon

As far as counting, I don't know of any, but snpEff can annotate SNPs, mapping them to genes and qualifying them as "HIGH", "MEDIUM", or "LOW" severity. You can then filter for SNPs that are meaningful to your project.

From there you could load the data into R, convert the gene data into a factor, then run the summary command to count the instances of each gene.

• 0 views • link
written by shawn.w.foley shawn.w.foley

1 answer

@zephyr_falcon maybe Maftools the R package could be useful (i could be wrong though); https://bioconductor.org/packages/release/bioc/html/maftools.html

Maftools accepts a standard .maf file as input but there is a way to convert vcf2maf; https://bioconductor.org/packages/release/bioc/vignettes/maftools/inst/doc/maftools.html#2_generating_maf_files

• 0 views • link
written by Alewa Alewa

Log in to answer this question.

More posts like this