Question: BCFtools variant calling format
gravatar for amy16
4 months ago by
amy1620 wrote:

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 • 261 views
ADD COMMENTlink modified 4 months ago by Ram13k • written 4 months ago by amy1620
gravatar for Ram
4 months ago by
New York
Ram13k wrote:

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.

ADD COMMENTlink written 4 months ago by Ram13k

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

ADD REPLYlink written 4 months ago by amy1620

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.

ADD REPLYlink modified 4 months ago • written 4 months ago by Ram13k
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: 1329 users visited in the last hour