This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BCFtools variant calling format

I run the command line below (as a job script on HPC):

samtools mpileup -uf Ca_v1.0_kabuli_ref.fasta B1_C.bam
bcftools view -bvcg B1_C.raw.bcf

I find the output file is generated in the slurm-out (the HPC specific job run output file). I am not sure if I am going wrong with the script somewhere.

rna-seq variant call
• 2,300 views • link
updated by Ram Ram • written by amy16 • amy16

1 answer

That's because you're not redirecting the output of the bcftools command, so it writes to STDOUT, which HPC systems dump to the file in their -o or -oo parameter. Add a >out.vcf to your bcftools command and everything will be fine.

• 0 views • link
written by Ram Ram

bcftools view -bvcg B1_C.raw.bcf >B1_C.vcf I tried it.. it is still not generating a separate output file

• 0 views • link
written by amy16 • amy16

The command is incomplete - -c needs an INT to work off of, and I'm not sure if -g and -v have defaults. Also, what is -b? What version of bcftools are you using?

It would also help if you looked at the STDERR file from the job.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this
  • Variant Calling (Samtools Mpileup): Error
    written by venuraherath • venuraherath
    Hi, I am trying to update a variant calling script came from the predecessors. I am getting the error "bash: mpileup : No such file …
  • Failed to open var.raw.bcf: unknown file type
    written by Ric Ric
    Hi, I used the following commands to call SNPs samtools merge ${1}/all.sorted.dedup.bam ${1}/*.sorted.dedup.bam samtools index ${1}/all.sorted.dedup.bam samtools faidx $2 samtools mpileup -uf $2 ${1}/all.sorted.dedup.bam | …
  • Failed to open var.raw.bcf: unknown file type
    written by Ric Ric
    Hi, I used the following commands to call SNPs samtools merge ${1}/all.sorted.dedup.bam ${1}/*.sorted.dedup.bam samtools index ${1}/all.sorted.dedup.bam samtools faidx $2 samtools mpileup -uf $2 ${1}/all.sorted.dedup.bam | …
  • samtools mpileup to bcftools mpileup
    written by Ric Ric
    I got samtools mpileup option `u` is functional, but deprecated. Please switch to using bcftools mpileup in future. How could I convert the following commands …
  • Identification of homozygous SNP
    written by olga_viktorovskaya • olga_viktorovskaya
    Dear Biostars Community, I am trying to identify SNP/indels in my WGS samples and faced an issue (I am just a newbie in bioinformatics). I …
  • SNP call using bcftools
    written by blur blur
    Hi, I am trying to do SNP call using bcftools, I tried running this in samtools: samtools mpileup -u -f CBS138.fasta 11_304.sorted.bam | bcftools view …
  • Calling SNPs with 1X depth (samtools/bcftools)
    written by Coryza • Coryza
    Hi, I am trying to call SNPs with at least 1X depth (My estimation is that most of the SNPs I want to see have …
  • Reference Genome Path In Samtools
    written by win win
    <p>Hi all, One of the classic commands in <a href='http://samtools.sourceforge.net/'>samtools</a> is as follows:</p> <p>samtools mpileup -uf ref.fa aln1.bam aln2.bam | bcftools view -bvcg - &gt; …
  • Mpileup -D And -S
    written by GPR GPR
    <p>I am calling out SNPs from a high coverage (&gt; 20fold) 75bp paired-end, &gt; 400 million reads data set. I have used a standard BWA&gt;MARK/REMOVE …
  • Bcftools Mpileup Removes Marked Duplicates?
    written by Angel Angel
    <p>Hey,</p> <p>My question is does samtools mpileup automatically removed duplicates from test.marked.bam output from Picard/MarkDuplicates? Because I am not following any other step to remove …