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

GenomicsDBImport : Mismatching intervals for input-vcf-files and bed-interval-file #1137

Closed
asp8200 opened this issue Jun 28, 2023 · 0 comments
Labels
bug Something isn't working PR opened

Comments

@asp8200
Copy link
Contributor

asp8200 commented Jun 28, 2023

Description of the bug

There seems to be a mismatch in the GenomicsDBImport process between the intervals for the input-vcf-files and the input-interval-bed-file. Here is an example:

gatk --java-options "-Xmx6553M" GenomicsDBImport \
    --variant HG002.haplotypecaller.chr16_46380683-90228345.g.vcf.gz --variant HG003.haplotypecaller.chr2_89753993-90402511.g.vcf.gz \
    --genomicsdb-workspace-path chr16_46380683-90228345.joint \
    --intervals chr16_46380683-90228345.bed \
    --tmp-dir . \
    --genomicsdb-shared-posixfs-optimizations true --bypass-feature-reader 

More info on Slack:

https://nfcore.slack.com/archives/CGFUX04HZ/p1686902281599399?thread_ts=1686493658.456969&cid=CGFUX04HZ

Command used and terminal output

nextflow run main.nf -profile docker --input /home/ubuntu/test_data/mapped_joint_bam_HG002_HG003_1perc_of_40x.csv --tools haplotypecaller --step variant_calling --joint_germline --outdir results --skip_tools haplotypecaller_filter


### Relevant files

_No response_

### System information

_No response_
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working PR opened
Projects
Archived in project
Development

No branches or pull requests

2 participants