This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Biopython - Extracting Nucleotide Frequency On Aligned Sequences

I was wondering how I might be able to use the BioPython package to extract nucleotide frequencies over an aligned set of sequences. Does anybody know how I might be able to do this? I am assuming that I have aligned sequences, including the appropriate spacers at places where base pair deletions have occurred.

Is there a single function in the BioPython package that allows me to do that? I've tried browsing through the API, but I haven't been able to find anything. (It might be I'm not looking at it correctly, or that the terminology used is different, perhaps?)

nucleotide-frequency alignment biopython
• 5,403 views • link
updated by Ram Ram • written by ericmajinglong ericmajinglong

2 answers

did you tried this Bio.Align.AlignInfo.Summary Calculate summary info about the alignment.

This class should be used to caclculate information summarizing the results of an alignment. This may either be straight consensus info or more complicated things.

• 38 views • link
written by Medhat Medhat

I was trying to do the same thing recently, and like you couldn't find anything. In the end I used pycogent which has quite a rich API for alignments. Take a look at this page for some toy examples of doing this.

• 33 views • link
written by cts cts

Log in to answer this question.

More posts like this