cn.mops on CRAM files
1
0
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!

cn.mops • 1.1k views
ADD COMMENT
0
Entering edit mode
@gunter-klambauer-5426
Last seen 3.8 years ago
Austria

Havent worked with CRAM files, but it should be possible to directly get read counts from them. Will check for solutions.

Maybe anyone from the community could help with this?

Regards, Guenter

ADD COMMENT
1
Entering edit mode

if cn.mops uses Rsamtools, then unfortunately CRAM is not yet supported but this is very much a feature for this release cycle.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

Traffic: 768 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6