Skip to content

bamutil diff - Get record for every read being compared #66

@DavidStreid

Description

@DavidStreid

Hi, I am wondering if I am using the tool correctly. Will there be a record for every QNAME in the bams being compared when comparing MAPQ scores?
I am noticing my diff output includes only a fraction of the QNAME entries of one of my bams being diff'd.

$ bam diff --in1 ${BAM1} --in2 ${BAM2} --mapQual >> diffs.txt
$ cat diffs.txt | grep -v "<" | grep -v ">" | sort | uniq | wc -l.  # get only the read names, no diffs
7199
$ samtools view ${BAM1} | cut -f1 | sort | uniq | wc -l.            # qname/read names in my input bam
46367516

I am trying to determine whether one map/align method mapped reads with greater confidence than another. Thank you for any direction, correction, or help. Thanks!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions