This is a test version of Biostars. For the public version, visit https://www.biostars.org.
From VCF to BEAGLE format with multiple chromosomes

I wish to convert a VCF file in the BEAGLE format with all chromosomes contained in the VCF file with vcftools. Unfortunately, it seems there is no option to use vcftools with all the chromosome, at least one chromosome has to be specified in the command.

vcftools --gzvcf GVCF_SNPs.vcf.gz \
        --out beagle_file.PL \
        --BEAGLE-PL \
        --chr chromosome_ID

If the chromosome is omitted, an error message shows up:

Error: Require a chromosome (--chr) when outputting Beagle likelihoods.

Would anyone know why it is the case? Or if there is any way around to process all the chromosomes?

vcf vcftools beagle gvcf
• 2,966 views • link
written by Begonia_pavonina • Begonia_pavonina

Hi there,

I'm trying to do exactly the same thing but not working. Did you figure it out?

Thanks!

• 0 views • link
written by marcela.herrerasarrias • marcela.herrerasarrias

Hello,

Yes, but it's not the same thing. I want to use my snp matrix from GATK and not genotype likelihoods from angsd. Did you manage to convert your vcf file to beagle binary using all chromosomes at the same time?

Thanks!

• 0 views • link
written by marcela.herrerasarrias • marcela.herrerasarrias

Unfortunately not, I think I have made some trials with PLINK, but didn't follow up on it. Good luck!

• 0 views • link
written by Begonia_pavonina • Begonia_pavonina

0 answers

No answers yet.

Log in to answer this question.

More posts like this