Question: BCFtools variant calling format
0
gravatar for amy16
13 months ago by
amy1630
amy1630 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 • 564 views
ADD COMMENTlink modified 13 months ago by RamRS19k • written 13 months ago by amy1630
1
gravatar for RamRS
13 months ago by
RamRS19k
Houston, TX
RamRS19k 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 13 months ago by RamRS19k

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 13 months ago by amy1630

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 13 months ago • written 13 months ago by RamRS19k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1524 users visited in the last hour