Bioinformatics connects life sciences and quantitative sciences, typically involving the application of software and algorithms to solve computationally intensive questions, such as those in genomics, sequence analysis, and systems biology. Questions tagged Bioinformatics should have a direct ...

learn more… | top users | synonyms (1)

3
votes
3answers
52 views

Needleman Algorithm for Optimal Alignment of two Amino Acid Sequences

I want to compute the optimal alignment of two amino acid sequences as per the following definition from a patent: "The percentage of identity between two peptidic or nucleotidic sequences is a ...
3
votes
2answers
158 views

Is there a database of cell images?

We're working on an algorithm for processing images of cells, similar to but much more basic than Cell Profiler, and we are looking for a large database of cell images to test our software. Can anyone ...
3
votes
2answers
137 views

Construct Phylogeny form Edgelist

I am working on large scale biodiversity simulations and have collected a lot of data on simulated phylogenies. I have a list of child-parent edgelist species in the form: ...
3
votes
1answer
4k views

What does phasing mean?

What does phasing mean in genetics/informatics? I've heard that a phased file is a file that has genes separated by chromosome, but can someone give a concrete definition of what phasing actually ...
3
votes
2answers
84 views

Microbiome Data

How is it possible to access microbiome data like that found here? I'd like to perform analyses on similar type of data, but cannot find OTU data like that at the Human Biome Project. Thanks for any ...
3
votes
1answer
67 views

Understanding common identifier codes (IDs) in biology databases

I know this is very basic but I would like to understand it and it's difficult to know where to begin with a new database. I have the output: ...
3
votes
1answer
223 views

How to perform a DNA structural alignment in pymol

How can I "fit" two DNA structures having different nucleotide sequences in pymol? I would like to use the structure of a DNA binding protein in pdb (1h9t), which is bound to DNA in the pdb file, ...
3
votes
2answers
97 views

Sequence compression method / format

As someone working with sequencing data and sequence files I was wondering whether is there a compression method / format for sequence / read files? Of course one can use general compressors like rar ...
3
votes
1answer
816 views

meaning of the “reads” keyword in terms of RNA-seq or next generation sequencing

I'm an undergraduate student at computer science and currently, I'm interested in bioinformatics. Today, I've started to read a paper about clustering and classification of non-coding RNAs can be ...
3
votes
1answer
737 views

What do the Clustal Alignment Symbols Mean?

Occasionally I will run protein alignments on peptide families and I can never remember what the symbols mean to show degrees of identity. What do they mean?
3
votes
1answer
53 views

Number of autoregulation and FFL motifs in a network

Suppose we have a network with N* nodes (N* is the number of internal nodes). Every directed link in this network exists with probability p. What would be the number of: motifs that are auto-...
3
votes
1answer
204 views

Any tool to align whole genome sequence data to another genome and give exon regions a higher mark?

I want to align WGS data from several mammals to one reference such as the human genome sequence. Since in most cases exons should be conserved and spliced in the same way and introns should vary, I ...
3
votes
1answer
1k views

Difference between “transcriptional activity” and “RNA expression”

I am working with a network object where the nodes are different events in molecular biology. Specifically; Catalytic activity of X, RNA expression of X, GTP-bound activity of X, Kinase activity of X,...
3
votes
2answers
72 views

Computing Percent Identity between DNA / Amino Acid Sequence

What are some tools where I can input a pair of DNA sequences (or alternatively a pair of Amino Acid Sequences) and compute a percent similarity identity metric between them? Is BLAST the right ...
3
votes
2answers
68 views

What is the most appropriate way to normalize gene expression data?

This question comes because reading a paper about normalization of gene-expression data, is not clear if the method for normalize the data is just for RNA-Seq data or could be applied also for ...
3
votes
1answer
56 views

constitution of read and gene region (IGV)

I work with fastq files containing NGS reads for some human DNA regions. The reference genome is hg19. I had two fastq files (pair-ended). I generated alignment BAM files. I used "bwa" and samtools to ...
3
votes
1answer
106 views

Why perform imputation?

Genetics datasets contain measurements for millions of single nucleotide polymorphisms (SNPs). Some (usually small) percentage of these values are of low confidence, and are labeled as missing values. ...
3
votes
1answer
45 views

Which hydrophobicity scales are best for detecting transmembrane regions, and why?

There are many hydrophobicity scales for protein analysis. Broadly, I gather the differences between them are from the experimental method to acquire the data and the normalisation (or lack thereof) ...
3
votes
4answers
660 views

Identifying the origin of replication of an unannotated *E. coli* plasmid

I have attempted a few searches for a list of origins of replication for plasmids in E. coli, but I was only able to find a list of origins, but not their individual sequences. The available plasmid ...
3
votes
1answer
2k views

Hamming distance between two DNA strings

By definition from Wikipedia, the Hamming distance between two strings of equal length is the number of positions at which the corresponding symbols are different. In other words, it is the number of ...
3
votes
1answer
104 views

PSI-BLAST website algorithm parameters

http://blast.ncbi.nlm.nih.gov/Blast.cgi In this website, when I want to apply the psi-blast algorithm on a sequence, under the section of algorithm parameters , what does PSI-BLAST threshold mean? ...
3
votes
2answers
534 views

How to use Autodock just for (re)scoring

I have some protein-ligand complexed that I have been docking with some other software and just want to use Autodock to evaluate those complexes. So, basically I just want to use it as a scoring ...
3
votes
1answer
166 views

Should the sum of RPKMs be constant over experiments

I have cufflinks data from 18 different RNA-seq experiments. I've noticed that if the sum all the RPKM values for a particular experiment, are vastly different between experiments. How could this be? ...
3
votes
1answer
2k views

RNA-seq analysis - q-values in cuffdiff

I'm using cuffdiff 2.1.1 to look for differential gene expression between two conditions. Each condition has 2 biological replicates. The results I get look promising from a log fold change ...
3
votes
1answer
206 views

How to determine the direction of regulation of a gene by comparing gene expressions?

I am just learning about the gene expressions and regulation. Several researches focus on finding the genes of altered gene expressions on a microarray to claim that they have a correlation to a ...
3
votes
1answer
180 views

Can compounds be predicted from Polyketide synthases gene clusters?

I've seen a few talks on the subject. From a gene cluster, there is some decent ability to determine which domains contribute to certain functional groups in the final metabolites. It sounds like this ...
3
votes
1answer
679 views

Displaying nucleotide at a single position from RNA-seq reads in a BAM file

How do I display a single nucleotide position from reads in a BAM file? I have been looking at variation using samtools mpileup, but I want to actually just display the nucleotide at the position I am ...
3
votes
1answer
62 views

How to predict a mRNA secondary structure with a large sequence?

When I use some web servers to predict a mRNA secondary structure, I find they always required in a small size sequence. If I use a long sequence and cut it into small parts, do these small parts ...
3
votes
1answer
60 views

blastn: What substitution matrix is used?

I'm currently working aligning sequences, and I need to compute similarity between pairs of DNA 'words' of a particular length. For amino acids I am able to use the substitution matrices in Biopython ...
3
votes
2answers
149 views

Deciding a reasonable threshold for copy number variation in a CNV (SNP array) TCGA dataset

Is there a methodology to select a reasonable threshold for copy number variation in a CNV (SNP array) TCGA dataset, to define when there is a significative alteration? Can I download CNV data for ...
3
votes
1answer
36 views

Gene Sequencing and Plasmapper

Is there anything similar to this in Java (especially the circular map sequencing along with hover effect)? For information I would like to convey that I am using Plasmapper and BioJava for achieving ...
3
votes
1answer
900 views

Why is Cysteine and Tyrosine used to calculate a sequence isoelectric point?

Why are the amino acids - cysteine and tyrosine used in isoelectric point calculations for a protein sequence, yet neither of them are positively charged molecules? and are not used in net charge ...
3
votes
1answer
117 views

SEREX serological analysis of cDNA expression library

What is Serological Analysis of cDNA expression library? I went through this article:http://cancerimmunity.org/serex/introduction/ but could not really make out. Can someone please explain this to me ...
3
votes
1answer
302 views

Software to calculate Fst from sequence data

I'm looking for a software to calculate Fst from 3 loci DNA data of individuals from a metapopulation. I don't have any prior on the population structure (have no idea of the number of subpopulations ...
3
votes
1answer
120 views

Regarding the NCBI FTP site

I am trying to extract information about SNP data from the FTP server of NCBI. Could someone please explain to me how the directory is organised? There are many many files and folders and I can't ...
3
votes
1answer
56 views

I2D vs STRING for protein-protein interactions

Can someone clearly explain me the differences between STRING and I2D(Interlogous interaction database) ? What does each database do best ? Which database is better for what kinds of studies etc ?
3
votes
2answers
109 views

Sequence alignment scoring

The following table is from Deonier's text Computational Gene Analysis at p. 152. This is an exercise in global sequence alignment and scoring of alternative sequences. The text proposed a solution (...
3
votes
1answer
96 views

Dnase data for GRCh38

This question below turned to be completely faulty. I don't have to do anything with DNase data for GRCh38. I asked it because of the file count difference between hg38 and hg37, which I thought to be ...
3
votes
1answer
49 views

Are there any major noticeable limitations to genome sequence compression methods that use reference templates?

Recently, I have been researching about big data analytics in biochemistry, and started wondering about how genome sequence compression could affect analysis. Of all the method listed on the ...
3
votes
1answer
77 views

Prediction of transmembrane beta barrels?

I studied that prediction of transmembrane alpha helices is more easy and accurate and also good algorithms are available for their prediction. But when we move towards prediction of transmembrane ...
3
votes
1answer
31 views

Gemini usage of --sample-filter

I am using gemini to query databases made from vcf files, which contain data from multiple samples. However, I need to query data only from one sample, (for which I currently use grep to filter the ...
3
votes
1answer
47 views

Finding confidence level of miRNA disease associations

I'm an undergraduate computer engineering student, and I have a project about bioinformatics. In this manner, I need to find prediction( or association I'm not sure the correct terminology) confidence ...
3
votes
1answer
233 views

mFOLD prediction result interpretation

I am trying to predict the secondary structure of certain predicted pre-miRNAs through mFOLD which is the generally accepted technique for structure prediction in most studies. I am finding it hard to ...
3
votes
1answer
241 views

Biological meaning/interpretation of clustering coefficience

I have been looking for some time to find an answer to the following question. I would be grateful for any help/advice or directions. I work very frequently with large number of proteins and assess ...
3
votes
2answers
96 views

Mass Spectrometry Terminology

I began reading this paper (http://www.biomedcentral.com/1471-2105/14/56) and had a few questions about mass spectrometry terminology that I couldn't find answers to elsewhere. Consider the following ...
3
votes
2answers
191 views

Why only heterogeneous SNVs for validation using genotyping arrays?

I am trying to validate the variants I found using whole genome sequencing . The standard practice, I have seen in the two publications below were to check for the number of heterozygous SNPs called ...
3
votes
1answer
102 views

What is a good source for info on disease frequency distribution among age groups?

I need information on the disease frequency distribution among age groups for an Android app I'm building... Hopefully I've come to the right place. Is there a good data source for this? Like a ...
3
votes
1answer
119 views

Superposing DNA

I have a series of protein models with DNA docked. I now want to superpose the DNA on a reference DNA molecule and extract the translational distance applied and the rotation angle used. I can ...
3
votes
0answers
29 views

filtering a BLAST search to only include a certain protein sequence

I am performing a BLAST search and I need to filter the output to a certain family of proteins. Specifically, I need to get matches within the CYP152 family of the Cytochrome P450 proteins. I can ...
3
votes
2answers
79 views

How much genomic variation one usually find within a given bacterial species?

If I find the exact starting position (say 1152471) of the coding sequence of a given gene in the genome of a bacterium, is the genome of the bacterium in general stable enough so that I can expect to ...