Issue with bcftools csq phasing
0
0
Entering edit mode
2.5 years ago

Hi all!

I want to perform phasing on a .vcf file (Human mtDNA sample). For this purpose I run bcftools csq.

My command is as follows:

 bcftools-1.6/bcftools csq -f Human_mito.fasta -g sequence.gff3 -p s sample_fltd.vcf -o trial.out

My gff file is as follows:

##sequence-region NC_012920.1 1 16569
##species https://www.ncbi.nlm.nih.gov/Taxonomy/Browser/wwwtax.cgi?id=9606
NC_012920.1 RefSeq  region  1   16569   .   +   .   ID=NC_012920.1:1..16569;Dbxref=taxon:9606;Is_circular=true;Name=MT;country=United Kingdom: Great Britain;gbkey=Src;genome=mitochondrion;isolation-source=caucasian;mol_type=genomic DNA;note=this is the rCRS;tissue-type=placenta
NC_012920.1 RefSeq  sequence_alteration 10  10  .   +   .   ID=id-NC_012920.1:10..10;Dbxref=dbSNP:879233578;gbkey=variation;replace=T

The error it prints is:

csq.c:715 gff_id_parse] Could not parse the line, "Parent=transcript:" not present: NC_012920.1 RefSeq  exon    577 647 .   +.  ID=exon-TRNF-1;Parent=rna-TRNF;Dbxref=GeneID:4558,HGNC:HGNC:7481,MIM:590070;Note=NAR: 1455;anticodon=(pos:611..613);codons=1;gbkey=tRNA;gene=TRNF;product=tRNA-Phe

Any idea how to solve this error? -Thanks

phasing genotyping • 1.0k views
ADD COMMENT
0
Entering edit mode

How are the contigs in your VCF named? You should also bgzip and tab-index your VCF, as this can resolve some header issues.

ADD REPLY
0
Entering edit mode

No I have not done that. I'll give it a try!

ADD REPLY
0
Entering edit mode

Okay. Please report back when you have time.

ADD REPLY

Login before adding your answer.

Traffic: 2464 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6