This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bedtools shuffle with adjustment of GC ratio

Hi All,

Is there any script which can be used for bedtools shuffle and can be applied to ‘GC' ratio control/adjustment ?

I want to random re-sampling same genomic regions with same length distribution and similar 'GC' contents.

Perl script will be preferred.

Thanks.

bedtools shuffle sampling gc contents

Hi Dr. Guo,

Wondering if you have found a feasible solution? Or how to define the GC% range? I met the same question too.

Thank you for your time.

1 answer

I would encourage you to look in to the below links for coming up with a solution. There is a bash script which you can convert to perl script to be honest.

Picking random genomic positions

Perl script: generate pseudo-CDSs from intergenic region

Log in to answer this question.