This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Variant calling results in polyploidy data

Hi,

I ran CRISP to do the variant calling in some Drosophila pooled samples and when I try to filter for minimum allele frequency with vcftools I get the following error:

VCFtools - 0.1.14
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
    --vcf new_variantcalls.vcf
    --maf 0.15
    --out berg

After filtering, kept 14 out of 14 Individuals
Error: Polyploidy found, and not supported by vcftools: YHet:49

Does anybody know how to fix this?

I run the grep command Pierre suggested and this is the output is here: https://jpst.it/MUgG

next-gen pool-seq crisp vcftools
• 4,797 views • link
written by mufernando • mufernando

Thanks for the heads up. I'll do that.

• 1 views • link
written by mufernando • mufernando

That is not possible as no one answered my questions.

• 1 views • link
written by mufernando • mufernando

I've moved the 'comments' to 'answers' for the two first items.

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

1 answer

looking at the vcftools source code: https://github.com/kiwiroy/vcftools/blob/master/cpp/vcf_entry_setters.cpp#L97

        if (pos != pos2)
LOG.error("Polyploidy found, and not supported by vcftools: " + CHROM + ":" + header::int2str(POS));

you have a strange genotype in CHROM= YHet at POS = 49. (something that doesn't look like '([0-9]+/[0-9]+)'. I suppose you have more than two alleles for a Genotype

please, show us the output of 

 grep -v "^#"  new_variantcalls.vcf | cut -f 10- | cut -d ':' -f 1 | sort | uniq

Does anybody know how to fix this?

change the way your VCF is called. ...

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

and as you can see, this is not a diploid genotype. You have

0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0

ploidy=48

while vcftools expects something like a diploid genotype ( e.g: 0/0 ) to calculate the MAF

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Would you say this is a problem with the variant caller or with sth that I performed before (e.g., alignment, local realignment)? If something could have happened before how can I check?

• 1 views • link
written by mufernando • mufernando

Would you say this is a problem with the variant caller

sorry, I don't know "CRISP".

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Hi I am getting the same error of "Error: Polyploidy found, and not supported by vcftools "6:6292579"

The output of

grep -v "^#" chr6.vcf | cut -f 10- | cut -d ':' -f 1 |sort | uniq

is giving : " 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 0|0 "

which means it is a diploid genotype only. What could be the reason for the error then?

• 0 views • link
written by shrutishreyajha • shrutishreyajha

Log in to answer this question.

More posts like this