I am trying to subsample from a bam file using the samtools view -s command. This is working when sampling 50% or lower (-s 42.50, 42 being the seed), but anything higher fails (returns an empty file).
He are the exact commands i use
samtools view -s 0.25 -b chr6_all.bam > chr6_25p.sam
samtools view -s 0.50 -b chr6_all.bam > chr6_50p.sam
samtools view -s 0.75 -b chr6_all.bam > chr6_75p.sam
I also made sure that 49% is working, but 51% is not. Any ideas, suggestions, or is this an intented mechanic ? There doesn't seem to by any documentation about the subsampling parameter in samtools docfile.