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

ERROR in reading reference index. #279

Open
Knnnk opened this issue Apr 11, 2024 · 2 comments
Open

ERROR in reading reference index. #279

Knnnk opened this issue Apr 11, 2024 · 2 comments

Comments

@Knnnk
Copy link

Knnnk commented Apr 11, 2024

Hi there,
When I am using octopus v0.7.4 to call germline SV from 1 bam file without forest model and seq_error_model, it throws an error:

[2024-04-11 17:06:33] <INFO> ------------------------------------------------------------------------
[2024-04-11 17:06:33] <INFO> octopus v0.7.4
[2024-04-11 17:06:33] <INFO> Copyright (c) 2015-2021 University of Oxford
[2024-04-11 17:06:33] <INFO> ------------------------------------------------------------------------
[2024-04-11 17:06:36] <EROR> A user error has occurred:
[2024-04-11 17:06:36] <EROR> 
[2024-04-11 17:06:36] <EROR>     The region Chr1:0-76456565 is after Chr19:0-51544995 in reference
[2024-04-11 17:06:36] <EROR>     index.
[2024-04-11 17:06:36] <EROR> 
[2024-04-11 17:06:36] <EROR> To help resolve this error ensure the region range format is <lhs> to
[2024-04-11 17:06:36] <EROR> <rhs> where lhs occurs before rhs in the reference index.
[2024-04-11 17:06:36] <INFO> ------------------------------------------------------------------------

octopus --version says

octopus version 0.7.4
Target: x86_64 Linux 5.11.0-1028-azure
SIMD extension: AVX512
Compiler: GNU 10.3.0
Boost: 1_74

what I knew is I am using a non-human genome, and the seq names aren't sorted, my command is as follows:

octopus -R /mnt/caller/ref/genome.fasta -I /mnt/caller/rgbam/230.bam -T Chr1 to Chr19 --threads 1 -o 230.vcf.gz
@Knnnk
Copy link
Author

Knnnk commented Apr 11, 2024

my chromosome names are

Chr2    Chr7    Chr13   Chr11   Chr19   Chr15   Chr18   Chr1    Chr4    Chr6    Chr8    Chr9    Chr5    Chr3    Chr10    Chr12   Chr14   Chr16   Chr17

@jelber2
Copy link

jelber2 commented Apr 11, 2024

Seems like you just need to remove -T Chr1 to Chr19 (which, I assume is that you want all 19 chromosomes to be variant-called?)

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