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

Excluded regions #339

Open
bounlu opened this issue Jul 27, 2023 · 1 comment
Open

Excluded regions #339

bounlu opened this issue Jul 27, 2023 · 1 comment

Comments

@bounlu
Copy link

bounlu commented Jul 27, 2023

I am using delly call with a BAM file with a reference fasta file different than originally used for the alignment. So the chromosomes are not compatible causing the below error:

delly call \
	-g genome.fa \
	-x human.hg38.excl.tsv \
	-o x.vcf \
	x1.deduplicated.sorted.bam \
	x2.deduplicated.sorted.bam

BAM file chromosome chr1_KI270762v1_alt is NOT present in your reference file genome.fa

However, I would expect this shouldn't be a problem since that chromosome is already in the excluded regions file:

$ grep chr1_KI270762v1_alt human.hg38.excl.tsv 
chr1_KI270762v1_alt

How to fix this in the best way?

@tobiasrausch
Copy link
Member

Alignments and delly have to use the same reference genome. This caused many problems in the past and that's why delly is very strict and requires the same reference genome.

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