Downsample BAM file to specific amount of reads

Downsample BAM file to specific amount of reads

0

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?


bam

• 30 views

updated 1 hour ago by

▴

260

written 2 hours ago by

▴

30

Read more here: Source link