JoseBlanca / vcf2fasta

Create sample sequences from a reference genome and a VCF file

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

I get a new problem

zhouyanli11 opened this issue · comments

Reading VCF file [Osjap10g14600.1.gff] ... done
Reading VCF file [japonica.combine.vcf.gz] ... done
Reading VCF file [all.chrs.con.fa] ... done
Ploidy is: 2
No phased genotypes found on first variant. Treating as "unphased"
Concatenating all [ CDS ]
Setting output directory to: vcf2fasta_CDS
vcf2fasta_CDS exists. Do you want to proceed? [y|n]: y
Total number of genes found: 1
Traceback (most recent call last):
File "/home/zhouyanli/vcf2fasta/vcf2fasta.py", line 481, in
main()
File "/home/zhouyanli/vcf2fasta/vcf2fasta.py", line 119, in main
sequences = getSequences(gff, gene, args.feat, args.blend, ref, vcf, ploidy, phased, samples)
File "/home/zhouyanli/vcf2fasta/vcf2fasta.py", line 216, in getSequences
alleles,max_len = getAlleles(rec, ploidy)
File "/home/zhouyanli/vcf2fasta/vcf2fasta.py", line 264, in getAlleles
max_len = max([ len(i) for i in segregating if i is not None ])
ValueError: max() arg is an empty sequence

The vcf file is too big,