I selected properly paired (and mapq>=4) reads from BWA alignment by running:
Flagstat shows that this did remove most of the reads that have mates on different chr. However, several remain:
Why is this?
Thanks.
samtools view -h -f 2 -q 4 -b
5 + 0 with mate mapped to a different chr
4 + 0 with mate mapped to a different chr (mapQ>=5)
4 + 0 with mate mapped to a different chr (mapQ>=5)
Thanks.