This is a test version of Biostars. For the public version, visit https://www.biostars.org.
mpileup Error cannot open file type

Hi! I have been working with next-gen data for a few months now (no training) and I have been trying to develop vcf files from a BAM and FASTA file. On one set of data I found that this command worked well and produced the file and data I am looking for: samtools mpileup -uf Reference_slow_base_Slow_Cmac-32w.fa Cmac_L_S_A_baseline_slow_mapping.bam | bcftools call -cv > var.raw2.vcf

However when I swap in a different FASTA and BAM file (that are aligned to each other), I get errors like unknown file type or failed to load. I have tried indexing the FASTA and BAM files and running them with the index and it still comes back with errors. Any ideas?

next-gen sequencing alignment
• 1,810 views • link
updated by Biostar Biostar • written by lewissechrome • lewissechrome

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • convert genome assembly to vcf file (Variant Call Format)
    written by cabraham03 • cabraham03
    is any way to convert an assembly of a genome in contigs to vcf file (variant call format) ?? I have some genomes in fastq …
  • Samtools consensus for WES files
    written by sidjain • sidjain
    Hi all, I am new to bioinformatics analysis. I have a sorted BAM file which is **WES**. In order to get chromosome 1 from the …
  • Failed to open var.raw.bcf: unknown file type
    written by Ric Ric
    Hi, I used the following commands to call SNPs samtools merge ${1}/all.sorted.dedup.bam ${1}/*.sorted.dedup.bam samtools index ${1}/all.sorted.dedup.bam samtools faidx $2 samtools mpileup -uf $2 ${1}/all.sorted.dedup.bam | …
  • Failed to open var.raw.bcf: unknown file type
    written by Ric Ric
    Hi, I used the following commands to call SNPs samtools merge ${1}/all.sorted.dedup.bam ${1}/*.sorted.dedup.bam samtools index ${1}/all.sorted.dedup.bam samtools faidx $2 samtools mpileup -uf $2 ${1}/all.sorted.dedup.bam | …
  • plink .bim file missing variant ID
    written by genomics Newbie • genomics Newbie
    Hello, First I generated a pileup from my BAM sequences and converted the pileup to vcf. samtools mpileup -uf ucsc.hg19.fasta sample_rna-seq.bam | bcftools call -mv …
  • Error in SNP calling using samtools and bcftools
    written by wu.zhiqiang.1020 • wu.zhiqiang.1020
    Dear all, I am using the newest version of samtools and bcftools (Version: 1.2 (using htslib 1.2.1)) to call the snps from sorted Bam files, …
  • Parallelized mpileup script
    written by Lee Katz Lee Katz
    Hi, I am wondering if anyone has made a script to make `samtools mpileup` run on one cpu per contig. I saw this idea on …
  • Error in vcf, samtools
    written by wu.zhiqiang.1020 • wu.zhiqiang.1020
    ``` ../software/samtools-1.1/samtools mpileup \ -uf ./S.int.fas \ ./S_new_F4_pe_mem.sort.bam | \ ../software/bcftools-1.1/bcftools view -vcg - > S_new_F4_pe_mem.sort.bam.raw.bcf [mpileup] 1 samples in 1 input files <mpileup> Set …
  • Error in SNP calling, samtools
    written by Adrian Pelin Adrian Pelin
    Hello, I have been following what most tutorials are suggesting (https://wikis.utexas.edu/display/bioiteam/Variant+calling+using+SAMtools), and I get: ``` @NCS ~/heterozygosity/TrimReads/bwa_Mapping_9iso_REF_GDR18/vcf/samtools_diploid $ samtools mpileup -uf gdr18_k75-85_NHC_conc.fasta ../../GDR-18.sort.grp.md.bam | bcftools …
  • Varscan Error On Mpileup Input
    written by GPR GPR
    <p>Hello, I am trying to filter low quality data off my mpileup files using Varscan but keep getting stuck with the error message: Invalid format …