Hi guys, My samtools version is: 0.1.19-44428cd
I am trying to use the '-s' option.. My bam file has: samtools view Example.bam | wc -l 2000000 samtools view -s 0.5 Example.bam | wc -l 1001528 samtools view -s 1.5 Example.bam | wc -l 998410 samtools view -s 2.5 Example.bam | wc -l 998744 I cant understand the results.. -s 0.5 should produce exactly 1000000 reads and 1.5 should be 3000000 etc.. Am i missing somthing? Moreover, samtools view -s 0.5 Example.bam > file1 samtools view -s 0.5 Example.bam > file2 file1 file2 are identical.. Isnt -s should be random subsampling??? Actually i saw this issue in some threads as solved, but it doesnt solve for me.. https://github.com/samtools/samtools/issues/160 (for example) Thanks, Pap -- ----------------- Dror Hibsh 0507-669599 ------------------
------------------------------------------------------------------------------
_______________________________________________ Samtools-help mailing list Samtools-help@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/samtools-help
