Posts

Showing posts from October, 2018

NCBI BLAST Graphical Summary using HTML/CSS

Image
This is a simple tutorial for designing a webpage for displaying the graphical summary of biological sequence alignment similar to the output of the NCBI BLAST+ program using HTML and CSS languages. In this tutorial, I have given few example hits (in bars) that match all range of bit scores (mentioned by color key).NCBI BLAST Graphical SummaryNCBI BLAST is a most popular bioinformatics framework for finding local similarity between two or more biological sequences. It addresses a fundamental problem in bioinformatics research. The heuristic algorithm it uses is much faster than other approaches, such as calculating an optimal alignment. This emphasis on speed is vital to making the algorithm practical on the huge genome databases currently available, although subsequent algorithms can be even faster.In NCBI BLAST, the graphic is an overview of the database sequences aligned to the query sequence. These are represented horizontal bars colored coded by score and showing the extent of th…

ATCG Content of Multiple Gene Sequences using C++

Image
This is a simple tutorial for computing ATCG contents of multiple gene sequences using the classes, objects, and members in C++. In this tutorial, I have used Dev C++ v5.11 software for compiling the C++ program.Program Design A. Define a class to represent a gene sequence data. Include the following members: Data members:Gene nameGene IDLengthA, T, G, C contentMember functions:To read data for a geneTo compute A, T, G, C contentTo display all the details of a gene B. Write a main program to test the program by reading n gene sequences data. Source Code/* Computing ATCG Content of Multiple Gene Sequences */ #include <iostream> #include <iomanip> #include <conio.h> #include <string.h> #include <stdlib.h> class gene { char r, gene_name[20][20], gene_seq[20][200], gene_id[20][20], id[20]; int gene_length, a, t, c, g, o, n, i, j, m, substr_eq, substr_rem; std::string seq, s; public: void get_sequence(); void find_length(); void find_at…