Jan. 25, 2013, 7:03 p.m. by Rosalind Team
Topics: Bioinformatics Tools
Let's Be Practical
If you are an accomplished coder, then you can write a separate program for every new task you encounter. In practice, these programs only need to be written once and posted to the web, where those of us who are not great coders can use them quickly and efficiently. In the Armory, we will familiarize ourselves with a sampling of some of the more popular bioinformatics tools taken from "out of the box" software.
To be equitable, we will focus mainly on free, internet-based software and on programs that are compatible with multiple operating systems. The "Problem" section will contain links to this software, with short descriptions about how to use it.
This initial problem is aimed at familiarizing you with Rosalind's task-solving pipeline. To solve it, you merely have to take a given DNA sequence and find its nucleotide counts; this problem is equivalent to “Counting DNA Nucleotides” in the Stronghold.
Of the many tools for DNA sequence analysis, one of the most popular is the Sequence Manipulation Suite. Commonly known as SMS 2, it comprises a collection of programs for generating, formatting, and analyzing short strands of DNA and polypeptides.
One of the simplest SMS 2 programs, called
DNA stats, counts the number of occurrences of each
nucleotide in a given strand of DNA. An online interface for
DNA stats can be found
Given: A DNA string
Return: Four integers (separated by spaces) representing the respective number of times that the symbols 'A', 'C', 'G', and 'T' occur in
20 12 17 21
Our default choice for existing functions and modules to analyze biological data is BioPython, a set of freely available tools for computational biology that are written in Python. We will give you tips on how to solve certain problems (like this one) using BioPython functions and methods.
Detailed installation instructions for BioPython are available in PDF and HTML formats.
BioPython offers a specific data structure called
Seqfor representing sequences.
Seqrepresents an extension of the "str" (string) object type that is built into Python by supporting additional biologically relevant methods like
In this problem, you can easily use the built-in Python method
.count()for strings. Here's how you could count the occurrences of 'A' found in a
Seqobject.>>> from Bio.Seq import Seq >>> my_seq = Seq("AGTACACTGGT") >>> my_seq.count("A") 3