bcftools command not found

[sadiexiaoyu]

Hi, All,

I am trying to generate a consensus sequence using samtools. And I used the command as follow:

samtools mpileup -uf ref.fa aln.bam | bcftools view -cg - | vcfutils.pl vcf2fq >cns.fq

But I get the error:

-bash: bcftools: command not found
Use of uninitialized value $l in numeric lt (<) at /Users/samtools-0.1.19/For consensus sequence/vcfutils.pl line 545.
Use of uninitialized value $l in numeric lt (<) at /Users/samtools-0.1.19/For consensus sequence/vcfutils.pl line 545.


I saw bcftools sub-file in samtools file, but there is no bcftools file in the sub-file. Could anyone tell me what is wrong?

Best,

Sadiexiaoyu

[GenoMax]

It appears that the directory containing samtools programs is not in default $PATH.

Have you tried to provide full file system path to the "bcftools" program in the above command?

[sadiexiaoyu]

It appears that the directory containing samtools programs is not in default $PATH.

Have you tried to provide full file system path to the "bcftools" program in the above command?

Hi, GenoMax,

Yes, I tried. I build PATH first to go to samtools directory. And I also tried to build PATH to go to bcftools sub file. But neither worked out...

[GenoMax]

Hi, GenoMax,

Yes, I tried. I build PATH first to go to samtools directory. And I also tried to build PATH to go to bcftools sub file. But neither worked out...

What version of samtools are you using? What do you mean by sub-file?

Can you post the output of

$ echo $PATH

and

$ which samtools
$ which bcftools

[dan]

Seems vcfutils.pl isn't part of samtools / bcftools 19... (it was present in release 18). Has it been wrapped into some other script or options for bcftools itself?

Cheers,

[GenoMax]

vcfutils.pl is in a directory called "bcftools" inside the main samtools directory.

[dan]

Thanks GenoMax, for some reason it's missing in this 'systems' install (bcftools is an exe, not a dir), but I see it when I re-download here:
http://sourceforge.net/projects/samtools/files/

I notice that it's the same file in release 18 and 19.