Entering edit mode
dsilencio
•
0
@dsilencio-20971
Last seen 5.4 years ago
Hello,
I have several hundred .cram files that I would like to analyze with cn.mops and would like to avoid converting to a .bam intermediate. Does anyone have experience using .cram as input or using an alternative tool to generate read count matricies from .cram files?
Thanks!
if cn.mops uses Rsamtools, then unfortunately CRAM is not yet supported but this is very much a feature for this release cycle.
Thank you all. I will convert to .bam for the time being. Out of curitosity, is it possible to merge the GRanges objects containing count matrices generated from getReadCountsFromBAM? I would like to batch my files, but might be nice to run the cn.mops analysis on a combined dataset.
Yes, you can stack them column-wise (samples) or row-wise (genomic segments). This is exactly what GRanges objects are meant for: easy handling of genomic information data.