Hello everyone:
I have a problem with samtools' mpileup.
Im trying to map a set of reads (only 1 library for now) against a reference genome, to find SNPs. The command line I used is:
$ samtools mpileup -uD -f ref_genome.fasta library.sorted.bam | bcftools view -bvcg library.raw.bcf
the output I get is this:
open: No such file or directory
[mpileup] 1 samples in 1 input files
Segmentation fault
i dont know what's wrong. I used this same command line with other data and it worked just fine, so I could use some help.
thank you!
I have a problem with samtools' mpileup.
Im trying to map a set of reads (only 1 library for now) against a reference genome, to find SNPs. The command line I used is:
$ samtools mpileup -uD -f ref_genome.fasta library.sorted.bam | bcftools view -bvcg library.raw.bcf
the output I get is this:
open: No such file or directory
[mpileup] 1 samples in 1 input files
Segmentation fault
i dont know what's wrong. I used this same command line with other data and it worked just fine, so I could use some help.
thank you!
Comment