Analysis of DNA, RNA or protein (amino acid) sequences. Here sequences are compared so analyze similarities and conserved regions. Sequence analysis is usually done by specialized softeware.

learn more… | top users | synonyms

6
votes
2answers
62 views

Predicting and identifying microbes and enzymes DNA sequence with metabolic prediction

Presently I am working on metagenomics of coal biomethenation by bacterial consortium. I have got the sequence result (Illumina). The sequence is huge and I can't predict anything from the sequence. ...
6
votes
1answer
68 views

(Rough) Model for DNA evolution in E. coli genes

I need a model of in-gene DNA drift. I'm not interested in bacterial phylogenies alone. Here is what I understand: Sequences corresponding to genes have both exons and introns, but in bacteria the ...
5
votes
2answers
2k views

What is the difference between local and global sequence alignments?

Obviously there are a bunch of different alignment tools out there, and I don't want to get bogged down in the maths behind them as this differs from software version to version. I am aware that ...
5
votes
3answers
244 views

Tool for nucleotide alignment with all nucleotide codes (e.g. R, Y, W, S, etc.)?

I have a vector sequence and would like to find the following nucleotide sequence in it. AASYWSRA This query sequence uses several degenerate symbols, defined ...
4
votes
1answer
53 views

Blosum matrix with probabilities instead of the positive and negative scores

I am trying to find a version of the BLOSUM matrix that has the frequencies instead of the scaled log-odds. i.e. instead of the common version that tells us that the score LEU/ASP is -4, I would like ...
4
votes
2answers
74 views

Reference sequence for defining single nucleotide polymorphisms

Single nucleotide polymorphism (SNP) or genetic variation in general, by definition are relative to a reference sequence. When we talk about databases of SNP as in ("dbSNP—Database for Single ...
4
votes
0answers
83 views

Comment on the introduction to a bioinformatics paper

I've written a paper about DNA sequence analysis. This paper attempts to use Bayesian modelling for a set of DNA sequences. It will probably end up either in a statistics journal, or, more likely, in ...
3
votes
1answer
346 views

Ka/Ks (dN/dS) analysis module for Python?

From this wiki article: In genetics, the Ka/Ks ratio (or ω, dN/dS), is the ratio of the number of Nonsynonymous substitutions per non-synonymous site (Ka) to the number of synonymous ...
3
votes
3answers
155 views

Is there a PSI-BLAST for nucleotide sequences?

I understand that one can translate a nucleotide sequence and run PSI-BLAST on the protein (proteins if you take the 6 reading frames), but I'm looking for distant homology for bacterial small RNAs ...
3
votes
1answer
64 views

What is the state-of-the-art algorithm for multiple sequence alignment?

Which algorithm or algorithms are considered the standard or state-of-the-art for multiple sequence alignment? How big is the need for better algorithms? How many sequences need to be alignment in a ...
3
votes
0answers
63 views

Data analysis of transcriptome sequencing data

I want to learn more about the data analysis and statistics on transcriptome sequencing data. I would like to read some important papers of the field and books and maybe some MOOCS, if they are ...
2
votes
6answers
161 views

Sequence analysis software suites

I have used DNAStar laser gene and VectorNTI in the past for cloning, primer design, sequence alignments, etc but no longer have access to these. Are there any comparable replacement suites or what ...
2
votes
1answer
36 views

Standard datasets for testing new multiple sequence alignment algorithms?

Are there any open and freely distributed standard datasets for testing new algorithms for multiple sequence alignment?
2
votes
1answer
112 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 ...
2
votes
1answer
158 views

How to check if a fastq file has single or paired end reads

I am trying to check if a fastq file has single or paired end reads. How can I achieve this with an error-proof method? I checked wikipedia and MAQ but I want to know if is there a reliable document ...
2
votes
1answer
40 views

Determination of gene direction

I have one sequence that composed of sequences of three genes. I should determine the direction of each gene, but two of these genes are disrupted, and I can't determine the initial codon. Does anyone ...
2
votes
0answers
101 views

Good poly-A filtering rules or tools

I am aligning a large number of ESTs. It seems poly-A tails show in many different ways. In addition to occurring at the very end, they can be flanked by the cloning sequence one one end, or have ...
1
vote
2answers
37 views

Recommended sequence clustering algorithm for transcriptome data

I'm working on a project where I'm going to analyze a large amount of transcriptome data. After assembling our RNA-Seq reads into contigs using Trinity, it looks like I'm going to have about 10GB of ...
1
vote
1answer
130 views

Codon alignment via Python? [closed]

I have pairs of coding DNA sequences which I wish to perform pairwise codon alignments via Python, I have "half completed" the process. So far.. I retrive pairs of orthologous DNA sequences from ...
1
vote
2answers
80 views

what does sequences overlap indicate

what does it mean if different protein sequences overlap upon multiple alignment, are they the same, do they have a common ancestor, are they recombinants of each other and how does one go on to ...
1
vote
1answer
47 views

What is the close and related genome used for in Gene models?

all I am a little bit confused about using the related genome or reference genome. When we have a reference genome, we can do alignment. Also we can do the assembly. Can you give some more reason ...
1
vote
1answer
129 views

Can PolyA tails be within Expressed Sequence Tags?

Can polyA tails occur within (rather than at the end) of a sequenced tag? Consider, for example the following two sequences from NCBI: DY008075 ...
1
vote
2answers
36 views

Hiding identical sequences in NCBI web interface

Sometimes there are identical sequences in Genbank/Genpept, most often coming from the same species but different strains. Ie if I look up "bacteria chitin synthase", the 624-aa protein from Dickeya ...
1
vote
0answers
20 views

Is the default kmer 25 value in Trinity still ok for assembly of 150bp reads from Illumina HiSeq 2500?

I will soon be receiving Illumina HiSeq 2500 data (150bp PE stranded reads). In the past we have been using Trinity as our assembler of choice, but it uses a default kmer value of 25. I believe this ...