bamcollate2 - split bam by readgroup

[keiranmraine]

Would it be possible to add the ability to split a BAM by the readgroups durring the collate phase similar to the way that bamtofastq does?

outputperreadgroup=<[0]>                     : split output per read group (for collate=1 only)
outputdir=<>                                 : directory for output if outputperreadgroup=1 (default: current directory)

Thanks,

[keiranmraine]

I mean to a set of BAM files rather than fastq.

[gt1]

Yes, but it may take some time.