Folks:
I have a mouse sequencing project with multiple .bam files generated.
samtools version 0.1.19-44428cd
bcftools version 1.1-137-g1f79bbd (using htslib 1.1-124-g9b1cb94)
I tried to do the following:
samtools mpileup -g -f mm10.fa mybam.bam > mybam.bcf
The program runs for several hours then generated a 5Gb file mybam.bcf
However,when I tried to use
bcftools view mybam.bcf
it has error:
[E::hts_hopen] fail to open file 'mybam.bcf'
[E::hts_open] fail to open file 'mybam.bcf'
Failed to open mybam.bcf: Success
same error message if I use 'bcftools call -v mybam.bcf' command
anyone has the same issue and is there a solution? the file is just there but seems the program fail to see it
I have a mouse sequencing project with multiple .bam files generated.
samtools version 0.1.19-44428cd
bcftools version 1.1-137-g1f79bbd (using htslib 1.1-124-g9b1cb94)
I tried to do the following:
samtools mpileup -g -f mm10.fa mybam.bam > mybam.bcf
The program runs for several hours then generated a 5Gb file mybam.bcf
However,when I tried to use
bcftools view mybam.bcf
it has error:
[E::hts_hopen] fail to open file 'mybam.bcf'
[E::hts_open] fail to open file 'mybam.bcf'
Failed to open mybam.bcf: Success
same error message if I use 'bcftools call -v mybam.bcf' command
anyone has the same issue and is there a solution? the file is just there but seems the program fail to see it
Comment