Calculating dinucleotide frequencies for each sequence in a multi fasta file
Bonjour,
I am looking for an easy way to calculate the dinucleotide frequencies of every sequence in a multi fasta file, and to have these saved with the sequence names associated with each frequency. Is there an easy command, tool, or way to do this using bash, or R, etc.?
I can think of some complicated methods of doing this, but if something is already programmed to do this then that would be amazing.
Merci pour votre temps,
Emily
2 answers
The kent utils from the UCSC genome browser have a rapid tool for this
faCount -dinuc genome.fa
Precompiled binaries are here.
There is a solution in this thread that should work: A long run time problem
Log in to answer this question.
More posts like this
-
Calculating empirical allele frequency
written by drowl1 •
-
InterProScan java issue
written by emilyc •
-
Gather sorting by parameter in the metadata of input file
written by grange •
-
Calculate allele frequency for structural variant breakpoints
written by nanana
-
Generate Random Dna Sequence Data With Equal Base Frequencies python
written by elisheva
-
Ensembl Perl API - How to get genotype frequencies
written by entheologist33 •
-
Script To Calculate Dinucleotide Frequency For Many Sequences
written by biolab
-
Estimating Probability Of Differing Allele Frequencies From Pooled Samples
written by Ketil
-
Split Data From A Csv Cell Into Smaller Cells
written by RossCampbell
-
Best Way To Calculate Variation/Snp Density In The Ensembl Database?
written by Jeremy •
Script To Calculate Dinucleotide Frequency For Many Sequences
Biostrings package in R (see page 72).
My understanding of this was that it's calculating the dinucleotide frequencies for all of the sequences over all, and not each. Am I incorrect in thinking this? I did not test this one but the OP said "I should note that I do not need to calculate each sequence's dinucleotide frequency. I need to get the di-nt frequency from all sequences. Thank you for your attention." in a comment.
Thank you.
You can iterate through your sequences as input and get those for each sequence. Or are you asking something else?
I mean yeah I can do this in a round about way, however, I am posting here to see if someone has/knows of a nicely done tool that fits my needs already with some organised outputs such as a dataframe with all frequencies for each seq with the original ID, etc. I can do the work, but if I don't need to then that is helpful.
Jellyfish is a popular tool for counting kmers in a multi-fasta file.