I have a bam file for the human whole genome. I want to split the bam file into multiple bam files, each having a fixed number of bases. How can I do this efficiently? I know I can generate a bed file with fixed regions, then run samtools view. Any other better way?