Message-ID: <20101220142704.34515wb2r69i2h7s@webmail.unibas.ch>
Date: 2010-12-20T13:27:04Z
From: Daniel.Berner at unibas.ch
Subject: [Bioc-devel] accessing gapped alignment
dear developers
I have performed a gapped alignment of my illumina reads. I now want
to analyze the output by accessing the CHROMOSOME and POSITION
information contained in the BAM file. I therefore uploaded the BAM
using readAligned(), but found that
all those reads aligning to the reference genome with a GAP are
missing. However, when using readGappedAlignments(), I loose
information crucial to my analysis (i.e., chromosome, position). Is
there any way to upload into R GAPPED
alignments while preserving read ID, chromosome, position etc?
thanks,
Daniel Berner