How to read .bam file in Rsamtools R package?

I have an accepted_hits.bam file which is from tophat output, and I need to use Rsamtools "R package" to do some analysis for the .bam file. So any one know how to read this file in R.

• 21,875 views • link

updated 4.9 years ago by Ram

• written 12.3 years ago by M K

1 answer

Did you try the package vignette?

http://www.bioconductor.org/packages/2.13/bioc/vignettes/Rsamtools/inst/doc/Rsamtools-Overview.pdf

• 0 views • link

updated 6.6 years ago by Ram

• written 12.3 years ago by Danielk

Hi Danielk, I used the Rsamtools package and I follow the the example in the pdf and it works. But I need to use my .bam file not the one provided in the pdf file. Also, I need to read all information in this file (not a portion).

• 0 views • link

updated 6.6 years ago by Ram

• written 12.3 years ago by M K

More posts like this

converting bam file to bed
written 7.6 years ago by azzopardiannalise •

Hi, I am trying to convert a bam file to a bed file using R software however, I am currently unsuccessful. The below is what …
BAM file processing in chunks in R
written 9.3 years ago by Tobias

When I want to process a FASTQ file in chunks in R, I can use the "FastqStreamer" function from the Bioconductor package "ShortRead". Is there …
Using featurecounts for EdgeR input?
written 10.4 years ago by nikelle.petrillo

Hello, I have paired end reads and I am performing de novo RNA seq. I used tophat to align my reads to a reference transcript …
Concordant pairs bam file larger than accepted_hits bam file?
written 11.7 years ago by bmpbowen •

I referenced [this post][1] on how to take tophat output from paired end RNAseq data (`accepted_hits.bam`) and use samtools to make a new bam file …
GTF annotation file for Human
written 12.3 years ago by M K

I want to run Tophat and I need to use the `-G` option to provide the Human annotation file .GTF file in my command. I …
How To Calculate Read Count For Antisense Study?
written 12.3 years ago by M K

Hi All, I am going to identify antisense for RNAseq using the R package NASTI-seq, but to run this package we should calculate read counts. …
Viewing .Bam And .Bai Files In Igv
written 12.5 years ago by HNK •

Hey I have run tophat in galaxy server on sample, and i got deletion, insertion, junction and accepted_hits.bam files. Now i want to view the …
What Kind Of Inputs Does Cuffdiff Take?
written 13.5 years ago by narges •

Hi, I have run tophat over the RNA-seq data and I got the output files : accepted_hits.bam, insertions.bed, junctions.bed, deletions.bed, unmapped.bam. Now to find …
What Does The Tophat File Named 'Accepted_Hits.Bam' Include?
written 13.8 years ago by narges •

Hi all, I wanted to ask about <code>accepted_hits.bam</code> file from TopHat. Is it correct that it does not contain all the valid alignments? If …
Empty Cuffdiff Output Files
written 13.9 years ago by narges •

Hi I have ran the cuffdiff over my BAM files derived from Tophat but all the output files are empty and I did not receive …