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
1
vote
2answers
150 views

Are codons in RNA layered? Are we misinterpreting RNA codons?

I am analyzing nucleotide base-pair patterns in RNA and DNA, and had a thought about RNA and DNA (Let me first state though, I am not a biologist; I am an algorithmatician, employing a sort of ...
0
votes
0answers
9 views

What format does Neutralization test data have? [closed]

Is there a template for virus neutralization assay data, as in what the information in the matrix is and how it should be presented?
0
votes
3answers
82 views

Finding motifs: fasta file with 10,000 sequences

I am new to Python. I am trying to parse a fasta file containing 10,000 sequences to look for motifs (microsatellites in particular). I tried using Seq Utils to parse my sequences for a particular ...
0
votes
0answers
10 views

Methods to analyze MeDIP-seq data

I have MeDIP-seq data. I have 6 groups with 15 replicates and a total of 90 human samples. I am trying to choose the best method to analyze them. I would highly appreciate if you can give me your ...
1
vote
1answer
16 views

How can I run the fdnadist command from the EMBOSS package?

My operating system is Mac OS X. I would like to compute DNA distances between different sequences. For this, the dnadist package from the Phylip software would be good but hard to automate. EMBOSS ...
1
vote
1answer
40 views

How to understand methylation level?

For Bisulfite Sequencing data, we can calculate methylation level for a specific locus, $$ \mbox{methylation level} = \frac{\mbox{# methelyated reads}}{\mbox{# total reads}} $$ Here I think, for this ...
1
vote
0answers
30 views

How to implement SEG (Wootton and Federhen ,1993) algorithm?

I want to implement SEG in my MATLAB environment and I'm relying on Wootton and Federhen (1993) Reading the article I cannot succeed to understand what kind of process I have to implement in order to ...
0
votes
2answers
26 views

Separating peaks of chip-seq with specific length

My data contain several chip-seq results. I have the peaks called by MACS2.I wanted to only look at those peaks that their size is e.g 500bp to 1000bp. How can I separate those peaks efficiently? I ...
2
votes
1answer
18 views

How to calculate the number of peaks that are upstream/downstream of some other peaks

I have 3 histone marks,I have used Macs2 for peak calling and diffBind to analyze the peaks. I was wondering if you know any way to calculate the peak numbers of one specific histone mark that occur ...
1
vote
1answer
23 views

DiffBind, diffferentially binding site

I have data for 3 histone marks (2 for silencing and 1 for activation) each mark has three replicates. when I run the diffBind package I have three contrast: ...
4
votes
1answer
504 views

Calculating average coverage for .bam files (sequence data)

(Full discolosure that this is my first time working with sequence data, and with the bash scripting.) I need to calculate the average coverage for any .bam file. After some searching I wrote the ...
1
vote
1answer
56 views

Do you know a program/script to count how many sequences have mismatches on each position of my primers?

I'm working with primer analysis and design, and so far I've performing it in Excel but: a) I don't like Excel, b) I do it manually, c) I don't like Excel, d) it gets tedious when working with ...
0
votes
2answers
36 views

differential analysis of chip-seq data

I have several sets of chip-seq data. I called the peaks using Macs2. I am pretty new to the field and I will appreciate any help. I wanted to annotate the peaks and see which peaks are shared between ...
0
votes
0answers
17 views

Generating taxonomic hierarchy by species/genus name

I am attempting to create a reference library of DNA plant barcodes in the ITS2 barcode region for plants from a specific region (Panama). I downloaded all the sequences for plants that resulted in a ...
0
votes
0answers
16 views

How to calculate overlapping between different chip-seq peaks

I have 3 sets of data for peaks called by Macs2 as below: 1-DNA methylation 2-Chip-seq (repression marks) 3-Chip-seq (activation marks) I need to show how many peaks in DNA methylation have the same ...
1
vote
0answers
15 views

Residue coevoution from position wise residue preferences rather than multiple sequence alignments

Hello experts in residue coevolution analysis of proteins, I am looking for a way to calculate residue coevolution (coupling scores) between all pairs of residues of a protein. The issue is that ...
1
vote
0answers
47 views

What steps should I follow for DNA analysis, for a classification problem?

I have a bed file which contains DNA sequences information as follow: ** ...
-1
votes
1answer
56 views

Compare nucleotide mutation with amino acid

To begin with I have two nucleotide sequence to compare I want to see how they might have evolved in terms of sequence similarity and differences. One of the way i can think of is align the sequences. ...
0
votes
1answer
52 views

How to change sequence format in my alignment file?

I have fasta file with alingned several sequences (from MUSCLE), when I open it (e.g. notepad++) they look like: And I want dot format of identities like and then save it in txt file. I failed ...
1
vote
0answers
26 views

Available Protein sequence alignment dataset and HMM model

It may better to move the question here. I am new to biology and I find my algorithm may be used in the Protein sequence alignment, since it is a henced HMM model. I find that people use HMM to ...
0
votes
0answers
19 views

Determining a Gene's Unknown Function

I'm a beginner in this field (so pardon my rudimentary knowledge/explanation) trying to understand ways on how I could potentially identify genes of unknown function. I am trying to get as close as I ...
1
vote
1answer
49 views

Evaluating amateur bioinformatics results

I've might have created a new algorithm for finding patterns in DNA. The technique has not been used before, and it managed to find a 33-mer with corresponding complement that exists both in E. Coli ...
0
votes
1answer
27 views

Arlequin files not able to converge beyond 2000 steps for some .arp files

When I use .arp files in arlequin to run mismatch distribution it says that can’t converge beyond 2000 steps and also no graphs are produced.can someone suggests what can I do to overcome this problem
1
vote
2answers
60 views

How to combine two Genome-wide Association Study (GWAS)? [closed]

I did a GWAS analysis in the past for antibiotic resistance of E. Coli and the results were interesting (matching the literature). I did a new GWAS analysis for some new samples, but the results are ...
0
votes
1answer
50 views

Is there a way to get the code from “github.mit.edu”?

everyone. I am a postdoc who just begins to analyze single-cell RNA-seq data. Nowadays, I found a really interesting paper (https://doi.org/10.1016/j.celrep.2018.10.047, PMID: 30404002). So, I ...
0
votes
2answers
44 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 ...
1
vote
1answer
25 views

After completing the evaluation, how to predict PTM lysine or any other modification with unlabeled data for making a web server?

I have created a model with SVM classifier and evaluated it with 5-fold cross-validation. The dataset was sequence data containing lysine modified sites. Now I want to build a web server where anyone ...
1
vote
2answers
62 views

Is it possible to go from a digital sequence to an actual nucleotide?

I'm trying to learn about what the state of technology is at presently. It seems like we clearly can go from nucleotide to digitally stored sequence, but can we transcribe something from the digital ...
1
vote
2answers
63 views

how to use list of gene id to get cds sequence(cds fasta file have many annotation, only gene id: is same to query id)

i have a question when i want to extract cds sequence using gene id. but cds file is not just start with >gene is, it has many other annotation. the only same is star with gene: cds fasta: ...
1
vote
1answer
54 views

Significance of k-mer length in COVID-19 sequence analysis?

I'm getting started in biology and bioinformatics with sequencing the SARS-Cov-2 or Coronavirus genome. I'm interested in this code which identifies k-mers in the genome: ...
1
vote
1answer
66 views

plot correlation between several bam files

I have several mab files from medip-seq. they supposed to be replicates and I have to draw their correlation. First I used multiBamSummary from deeptools to merge all the bam files based on bi/bed. ...
0
votes
2answers
360 views

merging two/or more bed file into one bed file

I am trying to merge two bed files (more in future) to one. my bed files are something like : . I need to merge them in a way to have the shared chromosome location. Is there a way to do that ?
8
votes
2answers
150 views

Understanding some of the computational bottlenecks of Covid-19 research

I am a researcher in high-performance computing (with very little bioinformatics background), and I am trying to understand what are the current biggest computational bottlenecks of software used for ...
0
votes
0answers
24 views

Coding object structure properties into a sequence

Hopefully, I found the right place to ask. Please, note that I'm not a specialist in the current field. Is there an algorithm to code information about object structural properties into a sequence (...
3
votes
1answer
43 views

Overlap in Paired-end Reads for Sequencing?

Regarding Sequencing: Do/can the paired-end reads have overlaps?
1
vote
1answer
47 views

Location of Reads in Sequencing

I have questions regarding sequencing: Are paired-end reads from different stands? What about single reads, are they coming from both strands?
1
vote
1answer
37 views

Warning in fastqc

I am checking the quality control of my sequences using the Fastqc tool. For some steps, like "Per base GC content", etc., I received a warning. So, I was wondering whether I Should also take care of ...
1
vote
1answer
48 views

SRA run data access

Trying to download run data of Run accession SRR2155174. Here's the link: https://trace.ncbi.nlm.nih.gov/Traces/sra/?run=SRR2155174 Please go to the link and go to the "Data Access" tab. Hereunder ...
0
votes
1answer
39 views

NCBI SRA database sample: control vs test

I was trying to download some data from NCBI SRA (SRA059451). There are 27 samples available for SRA059451. But i am unable to understand which samples are 'control' and 'test' samples. Please help me ...
2
votes
2answers
50 views

Tools or softwares for clustering evaluation

I am using a few clustering algorithms as well as my own tool for protein sequences. Also, I have benchmarked clustering results (i.e. ground truth) to compare with the results. Is there any tool/...
1
vote
2answers
45 views

Download sequences of isoforms

I want to collect all isoforms of all genes (Fasta/Fastq file of nucleotide and protein sequences) Arabidopsis Col-0. I am wondering if there is a straightforward way to download the file from any ...
0
votes
0answers
26 views

What kind of error-correction code is used in the error-correction code DNA sequencing?

In 2017, Chen et al reported a DNA sequencing technology called error-correction code (ECC) sequencing [1]. They borrowed ideas from communication and coding theory and encoded information redundency ...
2
votes
1answer
107 views

Any user friendly way to find rare mutations in whole genome raw?

Is there any user friendly way to find rare mutations in the individual human whole genome sequencing raw data? (from Dante, 30x coverage). To be more specific, I want to find mutations from this ...
0
votes
0answers
42 views

working with cutadapt

I'm working with ion torrent data where I apply the program cutadapt. I'm analyzing ITS seq data, as well as Matk. When I'm using cutadapt, I search the information for this genes, and built new ...
2
votes
1answer
184 views

How to convert VCF format file into Fasta(.fa) with the Reference data

I need a Software or Python Program for Converting VCF(.vcf) to FASTA(.fa) format with the help of reference (.fa)
0
votes
1answer
75 views

How to reproduce this specific sequence logo from 10 000 human mRNAs?

this is A sequence logo showing the most conserved bases around the initiation codon from 10000 human mRNAs. I would like to reproduce this sequence logo. How can I get the sequences? which database ...
0
votes
1answer
141 views

Is there a convention or standard of the highest bits in a sequence logo? [duplicate]

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). wiki says The total ...
1
vote
1answer
589 views

How to compute the Shannon entropy for a strand of DNA?

I'm confused by the computation of sequence logo. Wikipedia gives a process about this without a concrete example. Let's just consider DNA, so there are 4 bases (nucleic acids). The following data ...
1
vote
1answer
135 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 ...
1
vote
1answer
243 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 ...