This is a test version of Biostars. For the public version, visit https://www.biostars.org.
subset of paired end reads from bam file or fastq file of different depth

Hi,

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?

Thanks

depth shortread samtools read
• 1,725 views • link
written by rthapa • rthapa

1 answer

get the mean depth of coverage using , for example, mosdepth and then subsample the bam using the option -s <float> of samtools view

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Thank you. Do we need to specify anything for the paired end reads? The bam file is after mapping the paired end reads.

• 0 views • link
written by rthapa • rthapa

Log in to answer this question.

More posts like this