This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Comparing genomes against randomized genomes

Hi there!

I would like to have some advice in comparing a composition between a real genome against a randomized genome.

The question is about the randomization (as a background or Null model). When I randomize(schuffling the genome) a genome for comparing kmer composition It is important to keep the base frequency or a dinucleotide sequence?

I did read some papers but they used a expected value instead!

But I would like to look for a count method just similar to the kmer count method.

Any tip or paper using a similar approach would be appreciated!

Paulo

sequence genomes
• 2,382 views • link
written by schlogl schlogl

Why not compare actual genomes? DNA sequence in genomes is not random. What kind of significance do you expect to gain from unnatural comparisons?

• 0 views • link
written by 5heikki 5heikki

For transcription factor sequences or CpG islands, say, you cannot treat nucleotide frequencies as independent or shuffle bases. You might investigate hidden Markov model (HMM) approaches for generating simulated sequence, based upon categories of background regions.

• 0 views • link
written by Alex Reynolds Alex Reynolds

Got it. Thank you @Alex Reynolds

• 0 views • link
written by schlogl schlogl

1 answer

This may be of interest:

https://github.com/guma44/ushuffle

I used it to shuffle metagenome assemblies at a tetranucleotide level, and it does a great job of preserving the overall composition while completely changing the sequence.

• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Thank you @Mensur Dlakic 8)

• 0 views • link
written by schlogl schlogl

Log in to answer this question.

More posts like this