vg version: v1.41.0 "Salmour",
command line:
./vg construct -r genomic.fna -v bufo.12samples.vcf.gz -t 24 -m 32 > bufo.vg
error:
[tabix++] the index file is older than the vcf file. Please use '-f' to overwrite or reindex.
Describe:
I should clarify that the index for vcf is .csi insdead of .tbi because the bufo chromosomes are big so that can not be stored in tbi index. And I'm pretty sure that bufo.12samples.vcf.gz.csi is younger than bufo.12samples.vcf.gz.
Question:
How can I get rid of this error?
Thanks,
Maxine
1 answer
re-index
bcftools index --tbi --force bufo.12samples.vcf.gz
or
tabix -p vcf -f bufo.12samples.vcf.gz
tbi index just cannot build even force created.
bcftools index --tbi --force bufo.12samples.vcf.gz
[E::hts_idx_check_range] Region 961848..611782282 cannot be stored in a tbi index. Try using a csi index
index: failed to create index for "bufo.12samples.vcf.gz"
tabix -p vcf -f bufo.12samples.vcf.gz
[E::hts_idx_check_range] Region 961848..611782282 cannot be stored in a tbi index. Try using a csi index
tbx_index_build failed: bufo.12samples.vcf.gz
ah sorry, I was too fast ! is vg able to handle large chromosomes ?
I don't think so - they seem to check for a TBI or CSI file (https://github.com/vgteam/vg/blob/2c36d2cd08bef668eca2a8a64ba51abc037c905b/src/index_registry.cpp#L561) before the error message though, so I don't know what's happening.
OP is using a version released on 6-Jun-2022, and this line of code existed before then (whether or not it was executed is not something I can check). Only the devs can answer the question now, and they ask for posts to be opened here on biostars.
OP created a new issue on the GitHub repo. It seems like vg does not (and does not plan to in the immediate future) support CSI indexes.
Log in to answer this question.
-
sniffles process cannot end with warnings
written by Maxine •
-
VCF-merge fails due to tabix not producing .tbi files
written by brallen •
-
passing vcf.gz and vcg.gz.tbi files in tabix comand while doing range search
written by mandecent.gupta •
-
[Error] getting fasta file from 1000genomes
written by Mari •
-
Reindex A Vcf File / How To Handle Cosmic Vcf Files
written by secretjess
Please use concise titles. Add detailed error messages in the body of the post. I've changed it for you this time.
Just so you're 100% sure, can you run:
Yes
This might be a bug in their software - they may have not yet accounted for CSI indices. I think you should open a GitHub issue on their repository and point to this post.
A previous issue indicates that the software thinks a non-existent tbi file is somehow older than the VCF file - this is probably what's happening in your case, given that .tbi does not exist. vg needs to check for a .tbi or a .csi that's created after the VCF.
Sure! Thank you.