[Bioc-sig-seq] Rsamtools scanBam output to Bam

Sean Davis seandavi at gmail.com
Mon Apr 12 20:13:45 CEST 2010


On Mon, Apr 12, 2010 at 12:35 PM, burak kutlu <burak_kutlu at yahoo.com> wrote:
> Hi
> Is there a way to create a SAM file from the output of scanBam from Rsamtools?
> The reason we would like to do is use IGV to visualize a subset of regions in a large Bam file using load2newIGV function from SRAdb.
> For this, we will need to put the list structure back to SAM format.
> Please let me know if this is enough to give you an idea what we are trying to achieve.
> Thanks in advance

Hi, Burak.

You shouldn't need to subset the BAM file into regions.  IGV will load
HUGE BAM (we have loaded 200x coverage human genome) files without a
problem since it reads only the part of the file necessary and does
not store large files in memory.  And you will definitely want to use
BAM and not SAM for IGV.

Sean



More information about the Bioc-sig-sequencing mailing list