I want to subsample the reads of different depth from a bam file. For eg. 10X, 20X, 30X and so on. I tried with ShortRead R package but I can input only certain number of reads. I couldn't know how depth I am subsampling. Does any one have suggestion on how to subsample paired end reads of desired depth from fastq file or bam file?