When using samtools 1.19, I would index my .bam, and use idxstatst to make a file of read counts per sequence in my reference, and the unmapped reads would be at the bottom of the list.
The line would look like
* 0 0 1027820
To show 1027820 reads did not align
When using samtools 1.0 to index, now, the idxstats file always claims there is exactly 1 read that mapped to nothing, which is clearly wrong.
Can anyone else replicate this? This can't be the desired behavior, can it?
The line would look like
* 0 0 1027820
To show 1027820 reads did not align
When using samtools 1.0 to index, now, the idxstats file always claims there is exactly 1 read that mapped to nothing, which is clearly wrong.
Can anyone else replicate this? This can't be the desired behavior, can it?
Comment