Hello,
I obtained a BAM file from BWA and would like to keep only mapped reads. Samtools view -q 1 removes reads with mapQ 0 and I would like to conserve them.
Does anyone have a solution ?
Many thanks
I obtained a BAM file from BWA and would like to keep only mapped reads. Samtools view -q 1 removes reads with mapQ 0 and I would like to conserve them.
Does anyone have a solution ?
Many thanks
Comment