8000 Incomplete pre-processing when running from "--step prepare_recalibration"? · Issue #1875 · nf-core/sarek · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

Incomplete pre-processing when running from "--step prepare_recalibration"? #1875

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
DHelix opened this issue Apr 28, 2025 · 0 comments
Open
Assignees
Milestone

Comments

@DHelix
Copy link
DHelix commented Apr 28, 2025

Hi,

First of all, thanks so much for developing this pipeline!
I'm very new to variant calling and nf-core/sarek, so apologies in advance for any naive questions.

I'm working on normal-tumor pairs and have aligned BAM files with duplicates already marked.

When I started the pipeline from --step markduplicates, it finished without any errors and seemed to complete all pre-processing steps from marking duplicates to recalibration.
However, when I started the pipeline from --step prepare_recalibration, it also finished without any errors but seemed to stop before recalibration.

Here are the details from when I ran with --step markduplicates

nextflow run nf-core/sarek \
-r 3.5.1 \
-profile singularity \
-c sys.config \
-params-file params_ct01chr7p_mkDup.yaml \
--step markduplicates

params_ct01chr7p_mkDup.yaml

input: '/pipeline_test/nfcore_sarek/sample_sheet/ss_ct01chr7p_mkDup.csv'
outdir: '/pipeline_test/nfcore_sarek/output/ct01chr7p_mkDup'
genome: 'GATK.GRCh38'

ss_ct01chr7p_mkDup.csv

patient,status,sample,lane,bam,bai
test,0,normal,L1,/pipeline_test/data/ct01_chr7p/normal_chr7p.bam,/pipeline_test/data/ct01_chr7p/normal_chr7p.bam.bai
test,1,tumor,L1,/pipeline_test/data/ct01_chr7p/tumor_chr7p.bam,/pipeline_test/data/ct01_chr7p/tumor_chr7p.bam.bai

Here are the details from when I ran with --step prepare_recalibration

nextflow run nf-core/sarek \
-r 3.5.1 \
-profile singularity \
-c sys.config \
-params-file params_ct01chr7p_prepReCal.yaml \
--step prepare_recalibration

params_ct01chr7p_prepReCal.yaml

input: '/pipeline_test/nfcore_sarek/sample_sheet/ss_ct01chr7p_prepReCal.csv'
outdir: '/pipeline_test/nfcore_sarek/output/ct01chr7p_prepReCal'
genome: 'GATK.GRCh38'

ss_ct01chr7p_prepReCal.csv

patient,status,sample,lane,bam,bai
test,0,normal,L1,/pipeline_test/data/ct01_chr7p/normal_chr7p.bam,/pipeline_test/data/ct01_chr7p/normal_chr7p.bam.bai
test,1,tumor,L1,/pipeline_test/data/ct01_chr7p/tumor_chr7p.bam,/pipeline_test/data/ct01_chr7p/tumor_chr7p.bam.bai

Thanks so much for your help in advance!
Best,
Dan

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants
0