This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sorting with vcftools for tabix gone wrong

I concatenated two vcf files by scripting, keeping only the header from the first file and adding an extra column in the end to know which file each line is coming from.

I do this:

tabix -p vcf f1PotMut.vcf.gz

and get this:

[ti_index_core] the chromosome blocks not continuous at line 87352, is the file sorted? [pos 131430]

so I do this

vcf-sort f1PotMut.vcf.gz

but the problem remains. I get the same error when trying to tabex my "sorted" file. I somehow feel the vcf-sort command did not work. Any help?

vcftools tabix sort vcf index

1 answer

above you are not showing but for the sake of covering that: are you redirecting the output of the vcf-sort into a different file then are you running tabix on this new file?

Log in to answer this question.