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

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