Welcome to the New Seqanswers!

Welcome to the new Seqanswers! We'd love your feedback, please post any you have to this topic: New Seqanswers Feedback.
See more
See less

Subsample regions of too high coverage

  • Filter
  • Time
  • Show
Clear All
new posts

  • Subsample regions of too high coverage

    Hi all,

    Do you think there is an efficient way of downsampling only a few regions of a bam files (in my case the regions with a too high coverage).

    The idea, would be too randomly remove reads in regions where the coverage is above a given coverage.

    Indeed, in my analyses, those regions cause some steps of the pipeline to become really slow.

    Thanks for all your suggestions...

  • #2
    OK... I just ran into a great tool that seems to do exactly what I wanted. It is called VariantBam (,

    ./variant highcoveragebam.bam -m maxcoverage -o reducedmaxcoveragebam.bam -b

    I hope it will be useful for some of you.


    • #3

      Thank you for the suggestion. It works great!