[BioC] export a GappedAlignment object to a BAM file

[Yue Li]

Dear List,

I wonder if there's a way to export a GappedAlignment object to a bam file.

rtracklayer 'export' function exports gff, bed, bed15, bedGraph but not bam. It also does not recognize GappedAlignment class.

One could reformat the GappedAlignment object into a SAM-like table, write.table it to a local disk, and then use 'asBAM' to convert that file to BAM. But that seems a little cumbersome.

Thanks much,
Yue