Hi all,
We have been using Rsamtools to import BAM files also. On previous sequencing
runs the BAM files were around 800 megs and I could import the entire BAM with
.readAligned_bam(). Our new sequencing data creates larger BAM files and I am
unable to import them like before. So, as Steve pointed out to me offline, I
probably will need to import the reads from each chromosome individually.
Unfortunately it does not seem to be recognizing my chromosome names, error msg
below. Has anyone had a similar problem or recognize something that I might be
doing incorrectly? Thanks in advance.
> which <- RangesList(chr1=IRanges(start=1,end=247249719))
> params <- ScanBamParam(which=which)
>
> chr1reads <- scanBam("../aln.sorted.nodupes.bam",param=params)
Error in function (bam, tmpl, space, start, end) : failed to scan BAM
file: �
last record: 0
In addition: Warning message:
In function (bam, tmpl, space, start, end) : 'space' not in BAM header
file: ../aln.sorted.nodupes.bam
space: chr1
-------------------------------------------
Karl Dykema
Computational Biologist
Van Andel Research Institute
-----Original Message-----
From: [email protected]
[mailto:[email protected]] On Behalf Of Steve Lianoglou
Sent: Tuesday, January 26, 2010 9:59 AM
To: Aurelie Laugraud
Cc: [email protected]
Subject: Re: [Bioc-sig-seq] Rsamtools
Hi,
On Tue, Jan 26, 2010 at 9:41 AM, Aurelie Laugraud <[email protected]> wrote:
> Dear all,
> I have read about Rsamtools being in dev repositories and advices about
> waiting for it to be stable to use it ...
> my question is therefore how long before it comes as a "regular" R package ?
> Or is that that awful to make the dev version working ?
> I am using mac os X 10.6.2
For what it's worth, I'm using the version of Rsamtools from SVN in
order to read BAM files, and it works "as advertised". Note that this
means you (should be) using R 2.11 (you can get it from here:
http://r.research.att.com/).
-steve
--
Steve Lianoglou
Graduate Student: Computational Systems Biology
| Memorial Sloan-Kettering Cancer Center
| Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact
_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
The information transmitted is intended only for the person or entity to which
it is addressed and may contain confidential and/or privileged material. Any
review, retransmission, dissemination or other use of, or taking of any action
in reliance upon, this information by persons or entities other than the
intended recipient is prohibited. If you received this in error, please contact
the sender and delete the material from any computer.
_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
