Downsample BAM file to specific amount of reads
I’m attempting to downsample a bam file I have to 10 Million reads using seqtk.
For example:
seqtk sample -s100 myfile.bam 10000000 > myfile_10m.bam
However, when I run this, there are 0 lines in my bam file. This doesn’t seem to be correct.
I’m curious if there is another way to downsample a BAM file to 10 Million reads?
• 30 views
Read more here: Source link