Question: BCFtools variant calling format
gravatar for amy16
2.9 years ago by
amy1640 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 • 969 views
ADD COMMENTlink modified 2.9 years ago by RamRS30k • written 2.9 years ago by amy1640
gravatar for RamRS
2.9 years ago by
Baylor College of Medicine, Houston, TX
RamRS30k 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 2.9 years ago by RamRS30k

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 2.9 years ago by amy1640

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 2.9 years ago • written 2.9 years ago by RamRS30k
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: 854 users visited in the last hour