This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Random Genomic Regions Of A Given Length Distribution

Hi,

I'd like to know if someone has a script for generating a bed file with random genomic regions over the whole human genome. The regions should also follow a given length distribution as input (i.e., mean length = 3000 and some standard deviation,...). In general the length distribution should also be an empirical vectors with all the desired given lengths.

Thanks!

bed genome python
• 5,022 views • link
updated by Irsan Irsan • written by dfernan dfernan

Not sure what you want these regions for, but have you considered that some areas of the genome cannot be sequenced or have large sequence biases, such as occurrence of repeat regions, GC content, etc. So consider whether the random sequences are biologically relevent to you. You may want to consider filtering by repeatmasker or mappability tracks from UCSC.

• 0 views • link
written by Ian Ian

makes sense, thanks for the suggestion

• 0 views • link
written by dfernan dfernan

how do you control for GC content as well? Do you calculate the GC content for all the random regions of a length and see if it matches both the length and the GC content?

• 0 views • link
written by epigene epigene

I have the same question. How do you control GC content ?

• 1 views • link
written by Ming Tommy Tang Ming Tommy Tang

1 answer

Try random-genome-fragments from RSA tools or random.intervals from seqbias.

Seqbias allows you to give a vector of desired lengths. You can create that vector by sampling from a (e.g. normal) distribution like

lengths<-rnorm(mean=100000,n=10,sd=1000)

• 0 views • link
written by Irsan Irsan

thanks a lot. sounds like a great tool.

• 0 views • link
written by dfernan dfernan

BTW, you can also give random-genome-fragments a file with desired lengths, you can make that file by sampling from a distribution in R (like I previously said) and then writing the result to "desired_lengths.txt"

• 0 views • link
written by Irsan Irsan

Log in to answer this question.

More posts like this