I would like to find reads in one bam files which are not in a second bam file.
Is there an easy way to do this?
My way for now is:
- sort both bam files by read name
- extract read headers of both files (using
cut -f1) for easier handling
diffon both files
- manually search for differences
But this seems to me very tedious. Is there a better way?