Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

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.

1
vote
1answer
35 views

How to convert enrichment/depletion to frequency for comparing deep sequencing to sequence profile?

I have two datasets, from different sources, that I need to compare. The first set is deep sequencing results of a directed evolution experiment, where I have the naive library and selected library ...
3
votes
0answers
53 views

Are there any examples of weighted edit distance genome alignments?

It has been said that weighted edit distance is a preferred way to compare/align genomes in practice, e.g., when identifying genetically similar patients. I wonder if there are some concrete examples ...
3
votes
0answers
108 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
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
0answers
119 views

Difference between Classic Satellite DNA and Microsatellites/Minisatellites

I am very new to this community and have a basic question about repetitive DNA sequences. How is different (classic) satellite DNA from microsatellites/minisatellites? All textbooks give concrete ...
2
votes
0answers
59 views

How does one apply Bayesian inference to quantify a read the deeper you sequence?

For NGS sequencing technology, the "deeper" you sequence given fragments, the more certain you are of what is being sequenced. This sounds like a simple application of Bayes's Rule. What is the ...
2
votes
0answers
72 views

Genotype-Phenotype databases?

Beyond the Stanford HIV database, what other databases out there provide a dataset linking virus/bacterial genotype to quantitative phenotype? I'm looking for high quality datasets to test machine ...
2
votes
0answers
118 views

Efficient algorithm to calculate various population divergence statistics

Intro and description of the data I am simulating the evolution of very long DNA sequences. The model works well, is performant and will output data in the following kind of fasta format ...
1
vote
0answers
33 views

How to know which residues belong to the variable or constant regions of a Fab (fragment antibody)?

I have a fragment of antibody (Fab) structure, and I know its amino acid sequences of its heavy chain and light chain. As each chain of Fab can be divided into variable and constant regions, how can I ...
1
vote
0answers
124 views

How to calculate RPM?

I´m just starting my postdoctoral fellow and my new project involve sRNA Deep Seq. I´m having problems with the data normalization. I´ve been reading papers about it and the variety of normalization ...
1
vote
0answers
27 views

What are the exact reason causing missing calls in a VCF file from exome sequencing?

My data is a VCF file from exome sequencing variant call. I'm not very familiar with the sequencing process and variant calling process. I noticed that there are some missing genotypes, which is ...
1
vote
0answers
133 views

Why do we use klenow fragment in pyrosequencing?

Why do we use klenow fragment in pyrosequencing instead of polymerase
1
vote
0answers
59 views

Batch several sequences for absent restriction sites

I have a collection of about 120 7kB sequences I would like to check for ether a list of specific restriction sites, or what restriction sites might be absent in all of them. Is there a app or ...
1
vote
0answers
474 views

What are constrained genetic elements?

I am somewhat of a newbie in evolutionary biology currently taking my first steps in bioinformatics. I was reading a paper when I came across the term "constrained genetic elements", referring to ...
1
vote
0answers
227 views

Demultiplexing symmetric PacBio reads

I have PacBio CCS.h5 and the corresponding fasta and fastq files and I would like to demultiplex them. Does anyone know of how this can be done in the absence of bas.h5 files. Thanks for your help! ...
1
vote
0answers
537 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 ...
0
votes
0answers
11 views

Better algorithm than Burrow-Wheeler, for aligning short reads with high number of mismatches

Burrow-wheeler (BW) algorithm is probably the most widely used aligner out there. Largely it is popular because of it's speed and tendency to not take much memory. However, it has limitations when it ...
0
votes
0answers
11 views

Quantifying loss of mtDNA using genotyping arrays

I recently read that mitochondrial DNA (mtDNA) levels can be measured by relative real-time PCR. Would it also be possible to estimate the degree of loss of mtDNA from genotyping array data?
0
votes
0answers
70 views

C to T Transition Mechanism

I am seeing C to T transitions in my sequencing data, but also some strange signatures that follow along with that. I am trying to understand the biology of what is occurring. In the data I see more C ...