Questions tagged [sequence-analysis]

Umbrella term for understanding any given nucleic acid sequence code, either as a locus, loci or genome with itself or as a comparison to comparable nucleic acid sequence code either intraspecifically or interspecifically

Filter by
Sorted by
Tagged with
2 votes
1 answer
448 views

Why do ten rows (Figure_1) correspond to 2 bits (Figure_2) in a sequence logo?

Following this question, I'm confused with the computation of sequence logo Following data comes from the book "Machine Learning - A Probabilistic Perspective (Figure_1)" here is the corresponding ...
czlsws's user avatar
  • 169
3 votes
2 answers
3k views

Samtools sort: most efficient memory and thread settings for many samples on a cluster

We're preparing to analyze thousands of .bam files beginning with re-alignment, sorting, etc. Complicated question: Has anyone investigated the optimal thread and ...
Mark Ebbert's user avatar
  • 1,354
13 votes
1 answer
7k views

Better aligner than bowtie2?

Bowtie2 is probably the most widely used aligner because of it's speed. Burrow-wheeler (BW) algorithms (including bwa) tend to ...
user345394's user avatar
5 votes
2 answers
846 views

How can I locate duplicated regions in a sequence?

I am facing an issue when trying to align short reads against a region in human chr5. The two Sensory Motor Neuron genes, (SMN1 and SMN2) are almost 100% identical and this causes the aligner to fail ...
terdon's user avatar
  • 9,846
5 votes
7 answers
3k views

calculating nucleotide frequency per column

I have some sequences shown below ...
user3138373's user avatar
3 votes
2 answers
125 views

missed cleavages combinations generation

I am looking for an R based solution to the following problem: Trypsin does not cleave perfectly in proteomics experiments. So starting from a protein sequence FASTA I want to generate all possible ...
Klemens Fröhlich's user avatar
2 votes
1 answer
73 views

How do I measure the proportion of different bacteria in a sample from a high-throughput sequencer?

[this question is based on a question that was asked on Reddit] I have sequenced some [mostly cell-free] DNA from a sputum sample on a nanopore sequencing machine, and would like to know what is in it ...
gringer's user avatar
  • 13.9k
2 votes
1 answer
390 views

Samtools sort: most efficient approach to sort a single sample after aligning split .fastq files

Related to my other question (Samtools sort: most efficient memory and thread settings for many samples on a cluster), we need to optimize samtools sort as we ...
Mark Ebbert's user avatar
  • 1,354
2 votes
2 answers
85 views

Multiple Alignment cost application

Is there any biological application/case where someone would be interested in the estimated total cost of the alignment between a set of sequences (genes or amino acids) without the aligned sequence ...
gcorso's user avatar
  • 159
1 vote
1 answer
784 views

What is "aligned sequences" and "consensus sequence" in the context of sequence logo? How to compute these?

In bioinformatics, a sequence logo is a graphical representation of the sequence conservation of nucleotides (in a strand of DNA/RNA) or amino acids (in protein sequences). so, What is "aligned ...
czlsws's user avatar
  • 169
1 vote
1 answer
59 views

How to allow "half windows" in biopython?

I am using Biopython to examine protein sequences including the Kyte and Doolittle hydrophobicity. However, for short sequences, this returns an empty value; the value is not allowed if the window ...
James's user avatar
  • 409
1 vote
1 answer
287 views

reverse translation from amino acid string to DNA strings

what is the opposite of .translate() function calls ? I mean let's say I am given an amino acid string CYCLIC, how do I obtain all the possible combinations of DNA ...
kevin's user avatar
  • 141
0 votes
2 answers
417 views

How to make Venn diagram for the Macs peak calling output of two data sets?

I have two outputs of macs2 peak calling for two of my data sets. I wanted to plot the Venn diagram to see how many peaks are shared. I mainly work on Unix/Linux. Do you know any way that I can have ...
Mariam's user avatar
  • 115