Error "Failed to read from standard input: unknown file type" with Gzipped BAM Files
Nobody has claimed this yet.
Assessment
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Newbie friendliness
- 25/100
- Issue type
- Bug
- Clarity
- Needs clarification
- Activity status
- Stale
- Domain
- bioinformatics
Research direction
The report only identifies bcftools mpileup and call, the -b bam_files.txt option, and example .bam.gz inputs; it names no source file or test. Start by checking how the -b input list is parsed and how BAM paths are opened, then reproduce the reported error with the stated command pattern. Done means establishing whether the documented input format supports these files and recording the verified behavior or a narrowly scoped fix.
Written by the indexing model from the issue text.
Description
I’m using bcftools (v. 1.20) for variant calling with mpileup and bcftools call, but I’m encountering the following error message:
"Failed to read from standard input: unknown file type".
I’m using the -b bam_files.txt option, which includes a list of my BAM files (e.g., a.bam.gz, b.bam.gz, etc.). I suspect the issue is related to my BAM files being gzipped. Unfortunately, due to their large size, I don’t have enough space to unzip them.
Is there a way to use gzipped BAM files directly in the variant calling analysis? Any guidance would be appreciated.
Thanks!
- Dominant language
- C
- Stars
- 891
- Forks
- 277
- Avg merge
- 2d 23h
- Merged PRs (30d)
- 2
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
More from samtools/bcftools
-
needs-feedback
Difficulty 3/5 1-2 days Newbie friendliness 55/100
-
enhancement
Difficulty 5/5 Over a week Newbie friendliness 45/100
-
htslib-dependent P2: Desirable
Difficulty 3/5 1-2 days Newbie friendliness 48/100
-
Difficulty 5/5 Over a week Newbie friendliness 30/100
-
needs-feedback requires-test-case
Difficulty 3/5 1-2 days Newbie friendliness 55/100
All issues in samtools/bcftools
Similar issues
-
[adam] AdamNet network read doesn't cap to MAX_ADAM_PACKET_LEN, overflows client receive buffers Open
Difficulty 2/5 1-3 hours Newbie friendliness 88/100
FujiNetWIFI/fujinet-firmware#1649 · 2 comments ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 68/100
HarbourMasters/Shipwright#7229 ·
-
Difficulty 2/5 1-3 hours Newbie friendliness 72/100
riscv-software-src/riscv-isa-sim#2435 · 1 comment ·
-
bug Self Built Image SNAPSHOT Supported Device target/ramips
Difficulty 2/5 1-3 hours Newbie friendliness 88/100
-
Difficulty 2/5 1-3 hours Newbie friendliness 76/100