Error with bcftools

I tried to run the following command and got an error stating

bcftools: ccall.c:237: update_bcf1: Assertion `call->unseen==nals-1' failed. Aborted (core dumped)

This is my command:

samtools mpileup -C 0 -A -B -d 10000 -v -u -f hg19/hg19lite_ucsc_chrM.fa MB_1294_N_P0061.R1_sort.bam | /usr/bin/bcftools call -O v -v -c -n 0.05 -p 1 -A -o MB_1294_N_P0061.R1.vcf

How can I solve this problem?

Thanks!

bcftools samtools

• 1,943 views • link

updated 3.3 years ago by Ram

• written 10.3 years ago by cassandratzh •

Hello cassandratzh!

It appears that your post has been cross-posted to another site: http://stackoverflow.com/questions/36830782

This is typically not recommended as it runs the risk of annoying people in both communities.

• 1 views • link

written 10.3 years ago by Pierre Lindenbaum

0 answers

No answers yet.

More posts like this

bcftools get allele abundance
written 3.2 years ago by Fernando •

I'm using bcftools to extract variants from a bam file, but I have reference data that tells me whether the patient is homozygous or heterozygous. …
Error in running GenomicsDBConfigException
written 4.6 years ago by shivangi.agarwal800

Hi Guys, I am running GATK for variant calling using bam files for two samples (example). I first created gvcf using following commands: > java …
VCF samtools
written 4.7 years ago by Lucía •

Hello, I am having trouble when doing variant calling with samtools. I am getting only the header an no variants. If I would instead use …
change missing genotype in bcf file from 0/0 to ./. to run plink
written 5.4 years ago by rah •

I have a lot of bam files for which i need to call the genotype for, to create bcf files for and then run Plink. …
view variants in IGV between multiple samples.
written 6.2 years ago by akshitha4 •

Hello, I am trying to use 12 fastq file of tomato chloroplast against a reference genome. I am trying to use the samtools mpileup command …
bcftools call always assigns a genotype even if there are no reads?
written 6.5 years ago by Candah

Hi, I am variantcalling my 96 multi sampled bam file with samtools/bcftools however I obtain some strange results. When I use the following command: bcftools …
Genotyping with samtools/bcftools
written 8.1 years ago by SOHAIL

Hi everyone, I have BAM files of some ancient genomes (for example: 17X), authors mentioned in the paper that they genotyped all genomes individually using …
My output of variant calling is different by samtools and platypus
written 8.7 years ago by Anop Singh Ranawat •

samtools command samtools mpileup -u -f /usr/lib/cgi-bin/GenomaticMYELOID/ReferenceGenome/hg19/hg19.fa /home/anopperl6/Exomesoft/bam/sample_sort.bam | bcftools call -vmO z -o - > outputstudy.vcf.gz Platypus command python Platypus.py callVariants --bamFiles=/home/anopperl6/Exomesoft/bam/sample_sort.bam --output=/home/anopperl6/Exomesoft/bam/sample_sort.vcf --refFile=/usr/lib/cgi-bin/GenomaticMYELOID/ReferenceGenome/hg19/hg19.fa …
bcftools call segmentation fault
written 10.2 years ago by AP

Hi, I am currently trying to run bcftools call on a large bcf file and I keep receiving the "Segmentation fault" message. I use the …
Problem for pindel2vcf to create VCF file from pindel created _D files
written 11.7 years ago by hongbing.li •

Hi I am new pindel user, and use `../pindel -f hg19.fasta -i p1_config.txt -c ALL -o P2_Index_95_P319` to create `P2_Index_95_P319` with `_SI`, `_D`, and so …