I met problem when I work on a 7M bacterium genome. I mapped 140x of paired-end reads (151x2) from illumina to a Pacbio reference genome by BWA mem. Then used samtool to get a mpileup file. Finally using varscan to call indel and SNP. Below is commands I used: