BAM file processing in chunks in R
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 an analogous function for BAM files (for instance, in the "Rsamtools" package)? If so, could you give some toy code?
Many thanks for your help in advance!
• 3,356 views
•
link
updated
by
ivivek_ngs
•
written
by
Tobias
•
written
by
Tobias
1 answer
Have you tried to read the entire vignette of Rsamtools? There is a function call scanBam. Take a look at this link.
Alternatively if you already have bam chunks for specific chromosomes which you can obtain by splitting the total bams for each chromosomes and then process with scanBAM. Take a look at this SO link.
• 0 views
•
link
written
by
ivivek_ngs
Log in to answer this question.
More posts like this
-
filterAndTrim : Error in add(bin) : record does not start with '@'
written by Bioinfo •
-
Convert Samtools flags to R code
written by DanK •
-
[R, ShortRead] Subscript contains out-of-bounds indices
written by a.puwanat •
-
Why is there a difference between compressed and uncompressed files when reading it with the ShortR…
written by rantanplan82 •
-
R: Convert ShortReadQ class to DNAStringSet?
written by trodriguez2011 •
-
How to convert bam to bam.bai files
written by jain.harshika •
-
Inexplicable Cufflinks output FPKM variation
written by Tobias
-
Using CSAR package with ChIP-seq data: how to fix NAs generated in mappedReads2Nhits function
written by ivoryec •
-
R: Readaligned Only Junction Reads From Bam-File
written by lydia.herzel •
-
Processing Fastq Files In R/Bioconductor Without Reading Entire Files Into Memory?
written by Ryan Thompson
