Hi @ariadnaaterrades ! Can you please post some more details: command run, samplesheet, any other configs/param files you are using, and the nextflow log?

Hi @FriederikeHanssen , thanks for your quick reply.

I've been using mainly the default parameters with GRCh38 genome version , here is the command run:

cd /DNA/Batch_01
srun --mpi=none --mem=40G --nodes=4 --ntasks-per-node=2 --partition=long --pty bash -i

module load Java/13.0.2
module load singularity/3.7.3

./nextflow run nf-core/sarek --input ./samplesheet_10102022.csv --outdir /DNA/Batch_01/ --genome GATK.GRCh38 -profile singularity --wes --intervals /DNA/xgen-exome-hyb-panel-v2-targets-hg38.bed --tools mutect2,freebayes,cnvkit,strelka

The sample sheet looks like this:

<title></title>
<meta name="generator" content="LibreOffice 6.4.7.2 (Linux)"/>
<style type="text/css">
	body,div,table,thead,tbody,tfoot,tr,th,td,p { font-family:"Liberation Sans"; font-size:x-small }
	a.comment-indicator:hover + comment { background:#ffd; position:absolute; display:block; border:1px solid black; padding:0.5em;  } 
	a.comment-indicator { background:red; display:inline-block; border:1px solid black; width:0.5em; height:0.5em;  } 
	comment { display:none;  } 
</style>

And the thing is that there is no error in the execution_report_2022-10-11_13-00-30.html file...all the processes seem to be complete correctly. I did end up having output files, but not the somatic ones. For that reason I decided to look at the strelka tool information and I realized that it might not be possible to run tumor-only somatic variant calling.

Hey! Apologies for the late reply, I was OOO, but back now. So strelka is run with the following command for tumor-only samples:

https://github.com/nf-core/sarek/blob/master/modules/nf-core/modules/strelka/germline/main.nf

So with the "germline" module which is the only one we have available now as you correctly noted. However, there should still be the respective output present for you sample names. If that is the case, then the output is as expected and we can close the issue. The somatic ones can only be called with paired samples to my knowledge.

Hi!

As you mention, I ended up having germline results instead of somatic ones. However, I decided to run strelka among other methods to perform tumor-only somatic variant calling because in your website indicates that it would be possible to do so. https://nf-co.re/sarek/3.0.2/parameters

As you said, according to strelka's paper: https://academic.oup.com/bioinformatics/article/28/14/1811/218573 it is not possible to run it but when I saw your website information I thought that you may have implemented some changes in order to be capable of running it.

I think that changing the website information could prevent other people from misunderstanding.

Thank you so much!