-
Notifications
You must be signed in to change notification settings - Fork 607
Samtools mpileup [E::hts_open_format] Failed to open file .bam #1228
Copy link
Copy link
Closed
Description
Are you using the latest version of samtools and HTSlib? If not, please specify.
(run samtools --version)
Yes
Please describe your environment.
- OS (run
uname -sron Linux/Mac OS orwmic os get Caption, Versionon Windows)
Linux 5.3.0-46-generic - machine architecture (run
uname -mon Linux/Mac OS orwmic os get OSArchitectureon Windows)
x86_64 - compiler (run
gcc --versionorclang --version)
7.5.0
Please specify the steps taken to generate the issue, the command you are running and the relevant output.
I'm mapping sequences and then doing SNP calling, using the following steps:
minimap2 -ax map-ont ref.fas mydata.fq | samtools view -b | samtools sort - > outfile.sam
samtools mpileup -q 10 -aa -d 8000 -uf refs.fas outfile.sam | bcftools call -m > outfile.vcf
And get the error:
[E::hts_open_format] Failed to open file outfile.bam
[mpileup] failed to open outfile.bam : No such file or directory
Failed to read from standard input: unknown file type
I checked, my bam file is in the right folder and I can read it with samtools view.
I recently updated to samtools1.10, so I thought this might be the problem, but checked with an older version (samtools1.7) on files that used to work and it isn't working anymore. What can have happened?
And I know the u option is deprecated but bcftools mpileup gave me a similar error.
Thanks a lot!
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels