Question: [bcf_sync] incorrect number of fields (0 != 5) at 0:0
1
gravatar for david.4.ray
3.1 years ago by
david.4.ray10
United States
david.4.ray10 wrote:

I am attempting to call variants on some genomes and troubleshooting an error I'm getting from bcftools. My commands are as follows:

$SAMTOOLS_HOME/samtools mpileup -Q30 -C50 \
    -uf $REF_HOME/$REFGENOME$INPUT/$SPECIES"_R3_noDup.bam" \
    >$SPECIES"_mPileUp.vcf"

Followed by:

$BCFTOOLS_HOME/bcftools view -bvcg $SPECIES"_mPileUp.vcf" >$SPECIES"_variants_raw.bcf"

I get the following error:

[bcf_sync] incorrect number of fields (0 != 5) at 0:0

I've seen this error mentioned in other posts and most everyone has resolved it by adding -u to the initial mpileup run. However, I've already got that in the initial run.

The mPileup seems to run fine. In fact, I can run the following and generate a consensus.fq with no trouble.

$BCFTOOLS_HOME/bcftools call -c $SPECIES"_mPileUp.vcf" |\
    $VCFTOOLS_HOME/vcfutils.pl vcf2fq >$SPECIES"_consensus.fq"

Any ideas? If you come to my place I will buy you a cookie.

Thanks.

ADD COMMENTlink modified 17 days ago by zx87546.8k • written 3.1 years ago by david.4.ray10

Update:

I ran head and tail on my .vcf file and it doesn't look so good after all.

head yields what you would expect. Standard headers and such.

tail does something very strange. In fact, it appears to be trying to read a binary rather than a text file. No idea where this development is arising.

ADD REPLYlink modified 6 months ago by RamRS20k • written 3.1 years ago by david.4.ray10
0
gravatar for swbarnes2
3.1 years ago by
swbarnes25.0k
United States
swbarnes25.0k wrote:

If you want bcf format, shouldn't you use -g in the mpileup command? Without it, you are making ordinary pileup output, not vcf.

ADD COMMENTlink modified 6 months ago by RamRS20k • written 3.1 years ago by swbarnes25.0k
Please log in to add an answer.

Help
Access

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