"split chromosome" after sorting with bcftools
0
0
Entering edit mode
2.4 years ago
dec986 ▴ 380

I am running bcftools in a series of commands

First I run bcftools sort and then

plink --bcf tmp.sort.bcf --allow-extra-chr --make-bed --out merged.eigenstrat.34string.length.vcf.gz

which works, and then

plink --bcf tmp.sort.bcf --allow-extra-chr --recode --out merged.eigenstrat.34string.length.vcf.gz

but then I get an error:

Error: .bim file has a split chromosome.  Use --make-bed by itself to remedy this

this last error should have been remedied by bcftools sort but wasn't.

If I try to run the recode step based on the bed file, I get the error:

Options in effect:
  --allow-extra-chr
  --bed merged.eigenstrat.34string.length.vcf.gz
  --out merged.eigenstrat.34string.length.vcf.gz
  --recode

Error: A full .bed + .bim + .fam fileset is required for this.
For more information, try "plink --help <flag name>" or "plink --help | more".

How can I fix this error? sort didn't fix it, even though it should have.

bcftools plink • 701 views
ADD COMMENT

Login before adding your answer.

Traffic: 1622 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