There is a bam file with the pair-end reads alignment.
I want to exclude the reads only when neither of pair reads are mapped onto the genome.
What command I should type?
I have tried samtools view -f 2, which just removed the pair if either one is unmapped.( The same as samtools view -F 12 )
Besides, samtools view -F 4 or samtools view -F 8 will remove a read while the other is mapped.
Thanks!
Best regards!
I want to exclude the reads only when neither of pair reads are mapped onto the genome.
What command I should type?
I have tried samtools view -f 2, which just removed the pair if either one is unmapped.( The same as samtools view -F 12 )
Besides, samtools view -F 4 or samtools view -F 8 will remove a read while the other is mapped.
Thanks!
Best regards!
Comment