BCFtools variant calling format
1
0
Entering edit mode
7.0 years ago
amy16 ▴ 40

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 • 1.8k views
ADD COMMENT
1
Entering edit mode
7.0 years ago
Ram 44k

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 COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 1104 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6