Skip to content
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

DragMap fails if --trim_fastq flag is enabled #1445

Open
Ilaria-Pirona opened this issue Mar 20, 2024 · 4 comments
Open

DragMap fails if --trim_fastq flag is enabled #1445

Ilaria-Pirona opened this issue Mar 20, 2024 · 4 comments
Labels
bug Something isn't working

Comments

@Ilaria-Pirona
Copy link

Description of the bug

Good evening,
I am experiencing problems with Sarek's pipeline. When I try to remove the adapters from my WES data, the dragmap tool runs into "Segmentation Fault".
Do you have any idea what the problem might be?

Command line:
nextflow run main.nf
--input sample_sheet_3.csv
--genome GATK.GRCh38
--aligner dragmap
--save_output_as_bam
--step mapping
--outdir './results'
--save_output_as_bam
--wes false
-resume
-profile aws
--trim_fastq

Command used and terminal output

No response

Relevant files

No response

System information

No response

@Ilaria-Pirona Ilaria-Pirona added the bug Something isn't working label Mar 20, 2024
@FriederikeHanssen
Copy link
Contributor

What version are you running? Can you send the error message? there are several issues with dragmap that were never resolved: Illumina/DRAGMAP#47

@Ilaria-Pirona
Copy link
Author

The error is:

line 8: 115 Segmentation fault (core dumped) dragen-os -r dragmap --num-threads 12 -1 0002.11718FA_1.fastp.fastq.gz -2 0002.11718FA_2.fastp.fastq.gz --output-directory output --output-file-prefix 11718FA 2> 11718FA.dragmap.log

and I am using dragmap version:
root@f1bc5cc9b074:/# dragen-os --version
1.2.1

pulled from
quay.io/biocontainers/mulled-v2-580d344d9d4a496cd403932da8765f9e0187774d:5ebebbc128cd624282eaa37d2c7fe01505a91a69-0

@edmundmiller
Copy link
Contributor

Can you confirm it works without the --trim_fastq param?

@Ilaria-Pirona
Copy link
Author

Yes I tried multiple times, and when i added --trim_fastq it failed.

Sarek by default put the arguments
"--disable_adapter_trimming --split_by_lines 200000000 " in the process.

However, "--disable_adapter_trimming" to me is apparently in conflicts with "--detect_adapter_for_pe " which is written inside the module.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

3 participants