How to randomly extract pair reads from bam file？

Question

[Tao]

How to randomly extract pair reads from bam file？

Answer 1

[Pierre Lindenbaum]

see option -s of samtools view

  --subsample FLOAT      Keep only FLOAT fraction of templates/read pairs
  --subsample-seed INT   Influence WHICH reads are kept in subsampling [0]
  -s INT.FRAC                Same as --subsample 0.FRAC --subsample-seed INT