1
vote
2answers
35 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 ...
4
votes
1answer
48 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 ...
1
vote
1answer
46 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 ...
4
votes
0answers
77 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 ...
2
votes
1answer
143 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 ...
4
votes
2answers
72 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 ...
1
vote
1answer
127 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 ...
3
votes
1answer
330 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 ...
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 ...
3
votes
3answers
147 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 ...
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 ...
2
votes
1answer
108 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
0answers
98 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 ...
5
votes
3answers
231 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 ...