bcftools
bcftools copied to clipboard
Symbolic alleles not supported when using consensus _2
Hello!
Im using bcftools 1.15.1 :
bcftool consensus -H 1 -f ref.fa hifi_on_primary_phased.vcf.gz > hap1.fasta
error:
Symbolic alleles other than <DEL>, <*> or <NON_REF> are currently not supported, e.g. "<DUP>" at ptg000001l:19654749.
Please use filtering expressions to exclude such sites, for example by running with: -e 'ALT~"<.*>"'
grep: zgrep 19654749 hifi_on_primary_phased.vcf.gz
ptg000001l 19654749 pbsv.INS.DUP.4265 T <DUP> . PASS SVTYPE=DUP;END=19654820;SVLEN=71 GT:AD:DP:SAC 1/1:2,19:21:2,0,11,8
I read a similar error message in previous versions program, but the error was fixed ...
However, in the late version the same error...
The error message is correct, <DUP>
is not supported and never has been. I cannot remember what the behavior used to be, whether the program silently ignored them or crashed, but they never been used.
I will mark this as an enhancement.
I was confused by this message https://github.com/samtools/bcftools/issues/1176 ...