bcftools icon indicating copy to clipboard operation
bcftools copied to clipboard

Symbolic alleles not supported when using consensus _2

Open sashulkaSh opened this issue 2 years ago • 2 comments

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...

sashulkaSh avatar Aug 17 '22 12:08 sashulkaSh

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.

pd3 avatar Aug 17 '22 15:08 pd3

I was confused by this message https://github.com/samtools/bcftools/issues/1176 ...

sashulkaSh avatar Aug 18 '22 08:08 sashulkaSh