The Stack Overflow podcast is back! Listen to an interview with our new CEO.

Questions tagged [sequence-analysis]

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 software.

Filter by
Sorted by
Tagged with
35
votes
3answers
72k views

What is the difference between local and global sequence alignments?

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 not only between software but varies from software version to version. ...
12
votes
1answer
329 views

Which is the reference 16S rRNA?

Recently, I've stumbled upon a fact, which hasn't bothered me for many years. The fact is that all universal 16S primers are written as "[FR][0-9]+" (in regex notation), that is they have a position ...
11
votes
1answer
744 views

Sequence evolution simulation tool

I'm looking for a tool to simulate sequence evolution given a specific mutation model and birth-death model. I'm aware of tools and packages like INDELible, Seq-Gen and PhyloSim, but they simulate ...
10
votes
1answer
758 views

Standard practice for generating rarefaction curves from Next Generation Sequencing data

We have a few million 18S reads from a particular environment. The reads have been clustered into Operational Taxonomic Unit (OTU), and the OTUs annotated against a reference database. To generate a ...
9
votes
4answers
2k views

If we sequenced the genome of every species, would all phylogenies agree?

The Tree of Life is still up for debate. Most of this debate seems to be due to a lack of genomic information, but that deficiency is decreasing rapidly with advances in technology and sequencing ...
9
votes
1answer
209 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 ...
8
votes
2answers
1k views

RNA-Seq library construction challenges: the biases of RNA fragmentation vs cDNA fragmentation

I recently watched a presentation on RNA-seq that covered some of the choices one can make along the way, and I didn't fully understand one of the choices in particular. Near the beginning of the ...
7
votes
4answers
7k views

Significance of upper-case, lower-case and Ns in UCSC DNA files

I have downloaded human chromosome's data from UCSC FTP. Some part is in small alphabets and some is in large alphabets. Does it show the coding and non-coding region? Here is an example from the ...
7
votes
1answer
642 views

How can I generate a random DNA sequence?

I've found this paper which involves the construction of 19-bp random DNA sequences, but I don't know enough biology to understand how this method works. Could someone explain it to someone who is ...
6
votes
2answers
82 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
85 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
94 views

Comparative evolutionary study: is amino acid or nucleotide comparison more useful?

I am a high school student and am currently learning about evolutionary relationship study in biology. My teacher said that a comparative study of amino acid sequences is more useful than a ...
5
votes
2answers
1k views

Why is the quality range of fastq format so broad?

Referring to fastq format, it is clear that in fastq format, there are 94 quality value for a sequenced Nucleic Acid of a DNA sequence read and they are: ...
5
votes
1answer
136 views

Is it possible to establish a newer animal model completely based on Bioinformatics studies?

I wanted to know, if it is possible to prove a new organism as an animal model for any human disease by only Bioinformatics methods like NGS and Structural bioinformatics. For Example, let's say ...
5
votes
2answers
153 views

List of proteins by number of amino acids / chain length

Is there any protein database online where I could obtain a list of proteins ordered by the length of their chains / number of amino acids, starting from the shortest, as well as to see their amino ...
5
votes
1answer
922 views

Where can I find the pMON7124 plasmid sequence?

I need the map of the pMON7124 plasmid. As this is not easy to find, I welcome any hints.
5
votes
1answer
4k 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 ...
5
votes
1answer
904 views

What is the biological significance of k-mer counting?

There are many tools developed to compute the counts of k-mers present in a gene sequence. Jellyfish, Bloom Filter Counter, DSK Kmer Counter, KAnalyze, KMC 2 etc are some efficient software developed ...
5
votes
3answers
549 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 ...
5
votes
1answer
643 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 ...
5
votes
1answer
344 views

How do biologists infer correct ORF of a DNA sequence?

Each DNA (RNA) sequence has 6 possible Open Reading Frames(ORF). My question is: What are the theoretical bases of in vitro or in silico tries to find correct reading frame of a sequence? Is it just ...
5
votes
3answers
1k 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 as:...
4
votes
1answer
971 views

How to do multiple sequence alignment?

I have a DNA sequence that makes protein 1. but now I have asked to: compare the amino acid sequence of protein 1 with nine homologous proteins and make a multi sequence alignment (MSA) of the ...
4
votes
3answers
885 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 (...
4
votes
3answers
4k views

What is the difference between sequence alignment and sequence assembly?

I read the wikipedia page about sequence alignment and sequence assembly but I have not been able to find any difference between the two. What is the difference between sequence alignment and sequence ...
4
votes
2answers
198 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 ...
3
votes
1answer
1k 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
1answer
237 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
2k 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
2answers
56 views

An efficient and free method of assembling a large number of Sanger sequencing reads?

I have a rather large number of Sanger Sequencing reads.  Is there a automated method in either Ugene or a similiar freeware that allows you to quickly assemble forward and reverse reads. Ideally ...
3
votes
2answers
979 views

difference between motifs, domains, patterns, signatures and profiles

I can't get clear the difference between those terms, I see them a lot while browsing on Prosite, Pfam, Expasy e.t.c. However, I can find documentation about them, but It still not clear what the ...
3
votes
2answers
667 views

Basic Genotyping

As a disclaimer, I'm from a computer science / statistics background, working on some bioinformatics problems. I've got some genotyping data (VCF from Exome Sequencing), and I'm struggling to get my ...
3
votes
1answer
359 views

Is sequencing error a function of the nucleotide being read?

Checking out on Google Scholar, I can see that for Illumina (just to consider one example) the sequencing error rate is of the order of 0.001-0.01 per nucleotide. Talking about sequencing error, let'...
3
votes
0answers
111 views

How can I create R8 (homopolymer repeat) filter without using illumina pipeline?

illumina instruments have built-in -or online- analysis software for variant analysis (CASAVA). This software can filter out the false positive variants near the homopolymer repeats (AAAAAAAA) and ...
3
votes
0answers
102 views

Data analysis of transcriptome sequencing data [closed]

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 ...
3
votes
0answers
161 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 ...
2
votes
2answers
478 views

Change Genbank entry date with Biopython

I can create a new Genbank record in Biopython with the following code: ...
2
votes
2answers
44 views

Where I can download the genome sequence of various species of Fusarium?

I'm looking for genome sequence of various species of Fusarium (a kind of fungus!) on net, but I found no suitable information. Could you please help me out for this purpose? Thanks in advance
2
votes
2answers
255 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 ...
2
votes
6answers
477 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
44 views

Lysozyme amino acid sequence: N-terminal extension

I looked up the amino acid sequence of lysozyme here: http://www.biology-pages.info/L/Lysozyme.html Then I crossed referenced that with the lysozyme sequence on UniProt: https://www.uniprot.org/...
2
votes
2answers
75 views

Methods for microbial identification in soil

I'm trying to perform a survey of the life in a soil sample. I want to know what species of bacteria, fungi, and other organisms are in it. I've heard that ITS sequencing and 16S rRNA sequencing are ...
2
votes
1answer
129 views

Which DNA sequence will have higher melting temperature: CCCCCC… or GCGCGC…?

I tried this service to calculate temperatures. The nearest neighbour method consistently gives larger results for sequences like GCGCGC over sequences like ...
2
votes
2answers
1k views

Creating a Phylogenetic tree using DNA sequences [closed]

I was wondering given DNA samples, how can we create a phylogenetic tree? I mean I don't even want to have actual DNA sequence. Can I generate small samples, like reads of a length 100 or 1000 of ...
2
votes
1answer
99 views

Through which mechanisms did the Drosophila Histone Cluster evolve?

The Drosophila Histone Cluster is a gene cluster that contains more than a hundred copies of a sequence that encodes several histone proteins. These copies are very similar. My question is what ...
2
votes
1answer
227 views

Whole genome sequence analysis software

Please help to choose Bioconductor R packages and other software for the whole genome sequence data analysis and, in particlular, the goals of false discovery mutation rate, mutations exclusion, ...
2
votes
1answer
636 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
59 views

How can sequence analysis establish whether a genomic sequence is a functional gene or a pseudogene?

BLASTX of a query gene shows 5 hits in my subject, how can I find out which of these are expressed and which are pseudogenes (using DNA sequence or by sequence analysis)??
2
votes
2answers
71 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 ...
2
votes
1answer
6k views

What are Codominant vs Dominant Genetic Markers?

When talking about types of genetic markers, the adjective "dominant" and "codominant" are often used. I don't fully understand their definitions and found contradicting definitions. Foll and ...