Bcftools unsorted positions?
I am trying to sort a particular vcf file using bcftools for use for another software tool (MuSE). I tried running bcftools sort file.vcf > newfile.vcf, which ran without errors, but I get this message from the software tool: [E::hts_idx_push] Unsorted positions on sequence #24: 13474350 followed by 9950583. I'm not sure where these numbers are coming from, since I tried running grep 13474350 file.vcf and grep 9950583 file.vcf but don't see anything there. Am I missing something?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
samtools index error - "Unsorted positions on sequence" even though file is sorted
written by Geoffrey •
-
MethylDackel Error running on HPC server
written by magicridge •
-
samtools index error
written by pt.taklifi
-
Calling variants on reads with MAPQ=0 on HaplotypeCaller or bcftools mpileup
written by Anupriya Dalmia •
-
Merging multiple VCF files into one with bcftools
written by nadya •
-
Error while attempting to create a VennDiagram with bcftools stats
written by caro-ca •
-
How to fetch sequences from Proteinortho5 output containing all test species and no duplication in…
written by saadleeshehreen
-
How to specify the sort based on name in samtools sort?
written by anc.informatics •
-
[E::hts_idx_push] chromosome blocks not continuous :Indexing of gff3 files using tabix error
written by chevivien
-
What are the best resources or tools to sort an unsorted human gVCF
written by Irene@Sequencing.com
version of bcftools ? is there a dictionary ('##contig' lines in the header) ,
Try grepping for 13474351 and 9950584: VCF coordinates are 1-based while bcftools uses 0-based coordinates internally.
@Pierre Yes there is a dictionary, and it's version 1.9.
@chrchang Tried it but didn't output anything either...
Ok, just took a quick look at the relevant htslib source code (hts.c line 1851 in the current develop branch) and it does convert back to 1-based coordinates when printing an error message.
You may need to post an example file and command which can be used to reproduce the error to get useful help at this point.
This is really odd. Why would
bcftools sortthrow an error that says that the input file is unsorted? Is it not the job ofbcftools sortto do the sorting? I have a feeling that the file you're trying to sort and the file you'regrep-ing are not the same file.Are you piping stuff to
bcftools sort, perhaps?Sorry, I can see how the post is confusing. Edited for clarification.
Do you see any results when grep-ing for those positions in
newfile.vcf? Does MuSE say anything else?Did you get the
file.vcffile by subsetting a different VCF file using a bed file perhaps? If so, overlapping regions in the BED file could cause duplicate entries occurring in different positions in the resultant VCF file.