biod / sambamba

Tools for working with SAM/BAM data
http://thebird.nl/blog/D_Dragon.html
GNU General Public License v2.0
558 stars 104 forks source link

filtering uniquely mapped reads results in odd number of reads in PE data #410

Closed rikrdo89 closed 3 years ago

rikrdo89 commented 4 years ago

When using -F "[XS]==null" to get uniquely mapped reads, I get a odd number of reads, despite having paired-end data which should result in an even number of reads. is there a way to only get the pairs that are uniquely mapped? I tried adding "proper_pair" but it still resulted in an odd numbers of reads

pjotrp commented 3 years ago

It probably requires a new feature. Closing this until someone wants to write the code :)