Effective subsampling of BAMs
1
0
Entering edit mode
9 weeks ago
esimonova.me ▴ 10

I was trying to subsample WES BAM files with samtools view -s on coverage 75. My samples had mean coverage in target regions 195 and 130 so I used fractions 100/195 =0.51 and 100/130=0.76, when I subsampled them with samtools and recalculated coverage it was around 65-68. What I did wrong? If you have any ideas of other tools how to effectively subsample BAMs on the same number of reads or on the same coverage, please let me know.

BAM sampling • 108 views
ADD COMMENT
0
Entering edit mode
9 weeks ago
ATpoint 49k

I have some code here A: Subsampling Bam File With Samtools to use sambamba for the subsampling. You can specify a number of reads rather than a fraction, maybe that helps.

ADD COMMENT

Login before adding your answer.

Traffic: 2412 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6