Question: Issue with bcftools csq phasing
gravatar for shubhra.bhattacharya
20 months ago by
shubhra.bhattacharya120 wrote:

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

genotyping phasing • 722 views
ADD COMMENTlink modified 20 months ago • written 20 months ago by shubhra.bhattacharya120

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 REPLYlink written 20 months ago by Kevin Blighe61k

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

ADD REPLYlink written 20 months ago by shubhra.bhattacharya120

Okay. Please report back when you have time.

ADD REPLYlink written 20 months ago by Kevin Blighe61k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1980 users visited in the last hour