Sorting with vcftools for tabix gone wrong
1
1
Entering edit mode
9.6 years ago
axelwilhelm ▴ 120

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 • 7.6k views
ADD COMMENT
2
Entering edit mode
9.6 years ago

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?

ADD COMMENT
0
Entering edit mode

No, should I?

ADD REPLY
1
Entering edit mode

This guy does (bo-peng). Damn tools written by others.

ADD REPLY

Login before adding your answer.

Traffic: 3404 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6