Question: BCFtools variant calling format
gravatar for amy16
7 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 • 392 views
ADD COMMENTlink modified 7 months ago by Ram15k • written 7 months ago by amy1620
gravatar for Ram
7 months ago by
New York
Ram15k 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 7 months ago by Ram15k

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 7 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 7 months ago • written 7 months ago by Ram15k
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: 1679 users visited in the last hour