Skip to content

segmentation fault when using the -i flag in split-vep #1877

Closed
@sahwa

Description

test.vcf:

##fileformat=VCFv4.2
##FILTER=<ID=PASS,Description="All filters passed">
##contig=<ID=22>
##bcftools_viewVersion=1.17-4-gfe98a6b+htslib-1.17-1-g70df047
##bcftools_viewCommand=view -G -Ob ../data/imputed_v2.0_b38_chr22.bgz; Date=Tue Feb 28 16:28:27 2023
##bcftools_viewCommand=view -Oz /well/ckb/users/aey472/projects/annotate_rsids/data/imputed_v2.0_b38_chr22_nosamples.bgz; Date=Tue Feb 28 16:45:06 2023
##VEP="v107" time="2023-02-28 16:47:53" cache="/well/ckb/users/aey472/data/homo_sapiens/107_GRCh38" ensembl-variation=107.ab78711 ensembl=107.5f39899 ensembl-funcgen=107.0fbd7d5 ensembl-io=107.a473894 1000genomes="phase3" COSMIC="95" ClinVar="202201" HGMD-PUBLIC="20204" assembly="GRCh38.p13" dbSNP="154" gencode="GENCODE 41" genebuild="2014-07" gnomADe="r2.1.1" gnomADg="v3.1.2" polyphen="2.2.2" regbuild="1.0" sift="sift5.2.2"
##INFO=<ID=CSQ,Number=.,Type=String,Description="Consequence annotations from Ensembl VEP. Format: Allele|Consequence|IMPACT|SYMBOL|Gene|Feature_type|Feature|BIOTYPE|EXON|INTRON|HGVSc|HGVSp|cDNA_position|CDS_position|Protein_position|Amino_acids|Codons|Existing_variation|DISTANCE|STRAND|FLAGS|VARIANT_CLASS|SYMBOL_SOURCE|HGNC_ID|CANONICAL|MANE_SELECT|MANE_PLUS_CLINICAL|TSL|APPRIS|CCDS|ENSP|SWISSPROT|TREMBL|UNIPARC|UNIPROT_ISOFORM|GENE_PHENO|SIFT|PolyPhen|DOMAINS|miRNA|AF|AFR_AF|AMR_AF|EAS_AF|EUR_AF|SAS_AF|gnomADe_AF|gnomADe_AFR_AF|gnomADe_AMR_AF|gnomADe_ASJ_AF|gnomADe_EAS_AF|gnomADe_FIN_AF|gnomADe_NFE_AF|gnomADe_OTH_AF|gnomADe_SAS_AF|gnomADg_AF|gnomADg_AFR_AF|gnomADg_AMI_AF|gnomADg_AMR_AF|gnomADg_ASJ_AF|gnomADg_EAS_AF|gnomADg_FIN_AF|gnomADg_MID_AF|gnomADg_NFE_AF|gnomADg_OTH_AF|gnomADg_SAS_AF|MAX_AF|MAX_AF_POPS|CLIN_SIG|SOMATIC|PHENO|PUBMED|MOTIF_NAME|MOTIF_POS|HIGH_INF_POS|MOTIF_SCORE_CHANGE|TRANSCRIPTION_FACTORS">
##VEP-command-line='vep --input_file /well/ckb/users/aey472/projects/annotate_rsids/data/imputed_v2.0_b38_chr22_nosamples.vcf.gz --dir_cache /well/ckb/users/aey472/data --assembly GRCh38 -o ../output/imputed_v2.0_b38_chr22_annot.0Kb.vcf --force_overwrite --vcf --everything --offline --distance 0 --cache'
##bcftools_viewCommand=view /well/ckb/users/aey472/projects/vep_annotation/output/imputed_v2.0_b38_chr22_annot.0Kb.vcf.gz; Date=Thu Mar  2 16:16:48 2023
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO
22      10510497        chr22:10510497:A:G      A       G       .       .       CSQ=G|intergenic_variant|MODIFIER|||||||||||||||rs1457915203||||SNV||||||||||||||||||||||||||||||||||3.319e-05|0|0|0|0|0.0003159|0|0|0|0|0|0.0003159|gnomADg_EAS|||||||||
22      10519134        chr22:10519134:TGATA:T  TGATA   T       .       .       CSQ=-|intergenic_variant|MODIFIER|||||||||||||||rs1188336707||||deletion||||||||||||||||||||||||||||||||||0.0001|0|0|0.0009699|0|0|0|0|4.753e-05|0|0|0.0009699|gnomADg_AMR|||||||||
22      10519402        chr22:10519402:T:G      T       G       .       .       CSQ=G|intergenic_variant|MODIFIER|||||||||||||||rs1193222075||||SNV||||||||||||||||||||||||||||||||||3.128e-05|0|0|0|0|0.000846|0|0|0|0|0|0.000846|gnomADg_EAS|||||||||
22      10519607        chr22:10519607:C:G      C       G       .       .       CSQ=G|intergenic_variant|MODIFIER|||||||||||||||||||SNV|||||||||||||||||||||||||||||||||||||||||||||||||||||||
22      10525190        chr22:10525190:C:G      C       G       .       .       CSQ=G|intergenic_variant|MODIFIER|||||||||||||||rs1267940861||||SNV||||||||||||||||||||||||||||||||||2.493e-05|0|0|0|0|0.0006432|0|0|0|0|0|0.0006432|gnomADg_EAS|||||||||
22      10527415        chr22:10527415:C:T      C       T       .       .       CSQ=T|intergenic_variant|MODIFIER|||||||||||||||rs1321359077||||SNV||||||||||||||||||||||||||||||||||0|0||0|0|0|0|0|0|0|0|0|gnomADg_AFR&gnomADg_AMR&gnomADg_ASJ&gnomADg_EAS&gnomADg_FIN&gnomADg_MID&gnomADg_NFE&gnomADg_OTH&gnomADg_SAS|||||||||
22      10528424        chr22:10528424:G:A      G       A       .       .       CSQ=A|intergenic_variant|MODIFIER|||||||||||||||||||SNV|||||||||||||||||||||||||||||||||||||||||||||||||||||||
22      10531285        chr22:10531285:A:T      A       T       .       .       CSQ=T|intergenic_variant|MODIFIER|||||||||||||||||||SNV|||||||||||||||||||||||||||||||||||||||||||||||||||||||
22      10531362        chr22:10531362:AGATGAAAT:A      AGATGAAAT       A       .       .       CSQ=-|intergenic_variant|MODIFIER|||||||||||||||rs1279372490||||deletion||||||||||||||||||||||||||||||||||2.817e-05|6.702e-05|0|0|0|0|0|0|0|0|0|6.702e-05|gnomADg_AFR|||||||||
bcftools +split-vep -f '%CHROM\t%POS\n' -i 'Gene~"BRCA1"' -d test.vcf

Returns Segmentation fault

Same happens for the -e flag. The segmentation fault doesn't occur when the filtering flags are removed i.e.

bcftools +split-vep -f '%CHROM\t%POS\n'  -d test.vcf
22      10510497
22      10519134
22      10519402
22      10519607
22      10525190
22      10527415
22      10528424
22      10531285
22      10531362
bcftools 1.17-4-gfe98a6b
Using htslib 1.17-1-g70df047

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions