Hi Ben, This looks like it might be a bug. When I ask for the filespec for the 'merge' command: my @filespec = $new_bam->filespec; print join "\n", @filespec, "\n"; I get: obm *ibm (note the leading '*'). Could you please submit this as a bug? http://www.bioperl.org/wiki/Bugs Thanks, Dave