This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Probability of finding a 16S read in 1000 bp fragments

I have a statistical problem I need some help with. I was wondering how could I calculate what the probability of finding a 291 bp long V4 region of the 16S rRNA gene in any 1000 bp fragment of a whole genome dna extract would be. If the genome is 3,5 Mbp long.

In other words how many of the 1000 bp fragments would be expected to contain the said region?

Considering two cases: on average 2.7 copies of the gene in a genome or 5.4 copies.

sequence gene

You should take distance from origin of replication as a factor as well, depending on replication rate.

0 answers

No answers yet.

Log in to answer this question.