Issue with getPloidy

[rancilhac]

Hi,

I'm trying to run vcf2fasta.py to extract CDS from a large vcf. The GFF, fasta reference and vcf load fine, but then I get this error:

(vcf2fasta) [los@rackham3 fasta_CDS]$ more fasta_CDS.err
Traceback (most recent call last):
  File "/home/los/software/vcf2fasta/vcf2fasta.py", line 515, in <module>
    main()
  File "/home/los/software/vcf2fasta/vcf2fasta.py", line 94, in main
    ploidy = getPloidy(vcf)
  File "/home/los/software/vcf2fasta/vcf2fasta.py", line 490, in getPloidy
    return p[0]

My data is diploid and phased. Any idea why this happens, and how to fix it?

Thanks a lot,
Loïs