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

Hi,

I need to call somatic variants from a BAM file of cancer panel.

Can anyone please suggest any suitable tool for calling the variants and generate a VCF file.

Thank You

bam ngs variants cancer

2 answers

"Suitable" is very context-dependent, are you working with WGS or WES data? short o long reads? Check the Mutect2 pipelines from GATK and you must look for other variant callers in the literature.

It is an Ampliseq sequencing with custom bed file

https://github.com/nf-core/sarek

Sarek is a workflow designed to detect variants on whole genome or targeted sequencing data. Initially designed for Human, and Mouse, it can work on any species with a reference genome. Sarek can also handle tumour / normal pairs and could include additional relapses.

The pipeline is built using Nextflow, a workflow tool to run tasks across multiple compute infrastructures in a very portable manner. It comes with Docker containers making installation trivial and results highly reproducible.

Log in to answer this question.