Entering edit mode
7.1 years ago
lewissechrome
•
0
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?
