Hi there,
I'm using samtools 1.2, when I try to run this command:

samtools mpileup -B -f hg19.fasta file1.bam file2.bam file3.bam >
cohort.mpileup

I have the following error:

[fai_load] fail to open FASTA file.

After that, I tried to use the same command on reference chr1.fa (from ucsc
hg19) and it works perfectly. I downloaded ucsc.hg19.fasta already indexed
and I tried to build index again but faidx command doesn't work on
hg19.fasta reference and I have the following error:

[fai_build] fail to open the FASTA file path_to/ucsc.hg19.fasta

Do you have any idea? Why this command works perfectly on chr1.fa but not
on the entire reference.fasta?

Thank you for your help,

Matteo
------------------------------------------------------------------------------
_______________________________________________
Samtools-help mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/samtools-help

Reply via email to