[BioC] export a GappedAlignment object to a BAM file
gorillayue at gmail.com
Tue May 8 15:24:05 CEST 2012
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.
More information about the Bioconductor