Question: Get aligment quality / identity from CRAM file (for PacBio and ONT data)
gravatar for mschmid
19 months ago by
mschmid110 wrote:

I have PacBio and Oxford Nanopore reads mapped to a reference sequence. Since I have very long reads I had to use CRAM format.

Do you know a nice tool which calculates the aligment quality / mapping identity for every read from the CRAM alignment file?

Idealy it gives me also the the mapping length of the respective mapping.

ADD COMMENTlink modified 19 months ago • written 19 months ago by mschmid110

How did you align the reads? Most tools able to parse BAM can also parse CRAM, so maybe Picard CollectAlignmentSummaryMetrics may do what you want?

ADD REPLYlink written 19 months ago by h.mon27k

Partially. But I would like to have it per read. Not an average, to be able to plot a distribution of the alignment quality.

ADD REPLYlink written 19 months ago by mschmid110
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1721 users visited in the last hour