Split subreads.bam file into pieces

gistfile1.txt

# Extract all hole numbers / ZMW ids
pbindexdump test.bam.pbi --format cpp | grep basicData.holeNumber_ | sed 's/basicData.holeNumber_ = {//;s/};//' | tr ',' '\n' | uniq > zmws.list
# Count them
wc -l zmws.list
# Split into two lists
head -n XXXX zmws.list > part1.list
tail -n YYYY zmws.list > part2.list
# Create two files
zmwfilter --include part1.list test.bam part1.bam
zmwfilter --include part2.list test.bam part2.bam

# Sanity check to ensure that the merged files contain the same ZMWs
pbmerge -o merged.bam part1.bam part2.bam
pbindexdump merged.bam.pbi --format cpp | grep basicData.holeNumber_ | sed 's/basicData.holeNumber_ = {//;s/};//' | tr ',' '\n' | uniq > merged.list
diff zmws.list merged.list

# Sanity check to ensure that the merged files contain the records
samtools view test.bam > test.sam
samtools view merged.bam > merged.sam
diff test.sam merged.sam