Why is the total read number still more than the paired in sequencing after removing the duplicate in samtools flagstat output?
Bwa-mem may produce chimeric alignment, where different parts of a read are mapped to distinct loci. flagstat counts them as two reads.