I'm just trying to use samtools mpileup and it is aborting with out showing any error for some bacterial genomes, with the file given there are two variants and the program dies.
The error is just:
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-file depth to 8000
Aborted (core dumped)
The reference file (indexed with bwa) is: https://drive.google.com/file/d/0BzPjZ1hM-XWPMmxucnpDUjdpUFE/view?usp=sharing
The sorted bam file is: https://drive.google.com/file/d/0BzPjZ1hM-XWPUWw4NkdhM0E3Vm8/view?usp=sharing
The call is:
samtools mpileup -uf SaureusTCH1516_nt_genome.fasta 446.srt.bam
To obtained the sorted bam file:
First I trimmed pair end reads by quality and if adapters were found.
Second the trimmed reads are splitted in two different files, one for each end of the sequencing data, resulting in files without quality information.
Third I aligned the reads files to the reference with bwa aln and bwa sampe
Fourth I created the sam file with samtools view | sort
And finally i indexed the sorted bam file resulting in the file given in the link.
The version of samtools i'm using is 1.1 using htslib 1.1
Any hints on what i could be doing wrong, or what may be happening?
Thanks in advance