Questions tagged [dna]

Acronym for deoxyribonucleic acid, the unit of inheritance for all eukaryotic and prokaryotic organisms

Filter by
Sorted by
Tagged with
1 vote
1 answer
110 views

Problem with classification model of genomic data: every machine learing model predicts wrongly almost always the same subset of dataset

First of all, I'd like to apologize for any spelling or grammar mistakes. I'm having a problem using R for a classification problem. My dataset contains ~300.000 genomic data, and the features are ...
0 votes
1 answer
74 views

Help me in understanding the PDB file

Can anybody please explain the below line to me. "We took the structure and coordinates of nogalamycin from the X-ray structure determined in PDB code = 1D17" What do I need to download from ...
1 vote
1 answer
39 views

Output of allelic association doesn't write the rsID

I'm using plink (1.9b5) to do allelic association. My problem is that my output does not write the rsID (SNP), it just writes a dot. Output: ...
4 votes
1 answer
65 views

Counting hexamers in fasta sequence and identify its structure (and interruptions)

I have a lot of fasta files, each one with thousand of reads containing the hexameric motif "CCCTCT". The hexameric motif is highly continuous in most cases but interruptions may occur. I ...
0 votes
0 answers
40 views

What is the best tools to find all the somatic mutations?

I am a beginner in this field. I would like to know which is the best way to get all the mutations from a certain sample. At this point, I am considering using GATK (here) and Maftools (here). I don't ...
5 votes
1 answer
1k views

How to install DnaSP on a Mac

I need to install DnaSp, a nice GUI originally written for Windows for the population genetics of nucleotide sequence data. However, its not working on my mac. I'm using wine to install it using these ...
2 votes
1 answer
21 views

Where to get BED files for regions that contain binding site motifs for specific transcription factors?

Quick question, is there a place where one can download BED files for the binding motifs of certain transcription factors. In my case, I'm looking for a BED file for the genomic regions (enhancers) ...
6 votes
1 answer
83 views

Why does the SARS-Cov2 genome has letter t [duplicate]

ATTAAAGGTT TATACCTTCC CAGGTAACAA ACCAACCAAC TTTCGAT... is part of the 5'UTR of genome of an RNA virus SARS-Cov-2. RNA contains letters ...
3 votes
1 answer
26 views

Tajima-Nei Distance estimate with BioPerl

I have been trying to estimate Tajima-Nei distance for my data (link below). The output is s phylogenetics distance (pairwise) matrix which accommodates a detailed transition transversion and the ...
0 votes
1 answer
21 views

Independent Subset of Rectangles Problem for Nonoverlapping local alignments

I'd like to implement a graph-based algorithm for DNA comparison, more especially by solving a problem that can be formulated as a Maximum Weighted Independent Set problem. I have found the article in ...
1 vote
2 answers
274 views

Generating PDB file for custom DNA sequence

I asked this question in biology.stack exchange and was referred to this community. I am going to request migration for this question. https://biology.stackexchange.com/q/105024/67415 I want a PDB ...
0 votes
0 answers
17 views

Generating intercalation site in DNA for custom sequence

Thank you for your help. Can anybody please tell me how I can generate an intercalation site in DNA at the base pairs I am interested on. I appreciate if you can please guide me to any tutorials or ...
2 votes
1 answer
39 views

Simulation of DNA sequences through substitution rates

I'm looking for a little bit of guidance. My question is regarding the simulation of DNA sequences with a fix substitution rate. The majority of the programs for simulating sequences use Continuous ...
1 vote
1 answer
27 views

building complex drug-dna for AMBER software

I will appreciate if you can please clarify some of my doubts about drug-DNA complex. I want to study the drug-DNA simulation using AMBER. I did go through all the tutorials video on youtube but ...
0 votes
0 answers
26 views

Drug-DNA complex simulation using AMBER

I want to carry out the simulation of drug-dna complex by placing the drug molecule at particular site in my DNA sequence (intercalating site). I have PDB files of both drug and DNA sequence. I will ...
0 votes
0 answers
31 views

getting PDB file of interested sequence

I want a PDB file of DNA sequence I am interested in (10 base pairs) for docking. I will appreciate if any of you can please clarify some of my doubts. Do you think, I can generate PDB file by using ...
1 vote
2 answers
125 views

Compressing the human genome to few megabytes

Multiple sources (see for instance this or this) discuss how genetic data will have scalability problems, given the huge file size of the human genome. The most straightforward encoding (see here) of ...
-1 votes
1 answer
27 views

What is meant by transcriptional changes executed by the cell over a time period?

I read the following line in the research paper - The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells: During differentiation, for example, each ...
1 vote
1 answer
64 views

Translate all reads in a .fastq into protein sequence from deep mutational scanning experiment

I'm working with paired-end NGS reads from an Illumina platform. The sequences I have are all of the same gene, but have one or more substitution mutations each. Here is a rough workflow for ...
1 vote
1 answer
21 views

Encoding different DNAs in pySBOL3?

I want to encode two different DNAs using pySBOL3. First dna is an engineered region usually composed by a transcriptional unit (TU) or a set of TUs. The second is vector which is a basically a ...
0 votes
0 answers
48 views

How likely is it to find primers sequences in already-trimmed reads?

So, I am analysing 18S amplicon data (fastq files) to be able to eventually investigate the taxa composition. After removing primers (using cutadapt) from both R1 ...
0 votes
1 answer
46 views

How to know if the DNA sequence has been assembled and why is it important to know how it was assembled?

I have downloaded my FASTA format files, that have the DNA sequences of the coding region of the genes and the DNA sequence of the complete genome, from NCBI. How can I recognize if these sequences ...
4 votes
1 answer
121 views

Are there computational tools to extract features of DNA sequences?

I am looking for tools to extract features from short DNA sequences. For example, entropy, complexity, GC-content, etc. I have found the generateFeatures.py script from the PyFeat repo, but is there ...
2 votes
2 answers
93 views

What is the file .mums created by mummer?

I am doing genome sequence alignment using MUMmer, in particular I want to do a dotplot with mummerplot. So the passages that I did are: 1.create a file .mums with the following command line: ...
7 votes
2 answers
346 views

Theoretical limit of human genome compression

How small can a compressed file containing the human genome be? I'm aware that this question cannot actually be answered, since it is asking for the Kolmogorov complexity of the human genome, which is ...
1 vote
2 answers
121 views

How to identify the DNA sequence of a gene in the complete genome sequence FASTA format file?

I need to identify the sequence of a gene in the complete genome sequence . I thought it was simple, instead it is not a straightforward task ! My method was the following: I downloaded the FASTA ...
0 votes
1 answer
53 views

How to compute a p-value for the similarity measurement of two DNA sequences/contigs?

I'm attempting to come up with a statistical model to determine the p-value for the similarity value $ 0 \le s \le 1 $ of two DNA contigs. The idea is to get a numerical measure to estimate if two ...
1 vote
0 answers
16 views

How to interpret complete genome sequence in NCBI

I have a very basic and quick question: does the complete sequence of nucleotides of an organism (for example this one that I am considering in my homework of Pyrococcus Fuoriosus) reported in NCBI ...
2 votes
1 answer
82 views

Shannon's Entropy and DNA

I have been working with Shannon's entropy (SE) and DNA lately, and referring to the formula and concept of SE, I just wondered whether one should in the case of DNA use each nucleotide as having a ...
1 vote
2 answers
198 views

Fastq: how can I check if they are from DNA or RNAseq data?

I have (gave me) Illumina fastq files, which I want to use for variant calling, and I do not know if they are DNAseq or RNAseq data. How can I check this? I do not have any report or who to ask. Many ...
2 votes
4 answers
264 views

What are the applications of DNA or RNA pattern matching? [closed]

I'm assuming we don't do pattern matching in DNA or RNA for the fun of it. So I'd like to know what are the applications of pattern matching or where does it fit in in larger applications? I'm a ...
2 votes
1 answer
116 views

How can I use my Myheritage DNA results file for further analysis?

I had my DNA tested by Myheritage and they sent me a csv file with RSID, Chromosome, position and result (which base) with about 700,000 rows. I understand most analyses of DNA use VCF files but is ...
0 votes
0 answers
180 views

Convert DNA raw data between a different formats [duplicate]

I have .vcf file with partly decoded genome provided by some service. How can I convert this file into other popular formats which may be used in other services? I'...
0 votes
1 answer
30 views

Dataset showing evolution of gene

As a way to exemplify how the hierarchical clustering of particular DNA sequences can be used to derive phylogenetic trees, I'm looking for a (even very small) dataset of the same genes (or other ...
9 votes
3 answers
4k views

What is a quick way to find the reverse complement in bash

I have a DNA sequence for which I would like to quickly find the reverse complement. Is there a quick way of doing this on the bash command line using only GNU tools?
-3 votes
1 answer
47 views

Does anyone know how I can convert DNA code into FASTA for this TTGAAACACTGGATGAATGAAAAGCCCTGCTTTGCAACCCCTCAGC [closed]

TTGAAACACTGGATGAATGAAAAGCCCTGCTTTGCAACCCCTCAGC this is the DNA code Sequence. I have tried converting each into the amino acid and ended up with this LKHWMNEKPCFATPQX but I was told that this isn't ...
0 votes
0 answers
52 views

Pairwise alignment using CIGAR String

I want to perform pairwise alignment and get the resulting CIGAR string. The channel can emit insertions, deletions and substitution errors. For example, if the input was AAGCT Then the output could ...
0 votes
0 answers
24 views

DNA sequence error annotation

Suppose I generate a DNA sequence of the following pattern. AAGTC And after being passed through a channel with insertion, deletion and substitution errors obtain the following sequence AAAGGC Where ...
5 votes
3 answers
2k views

How to translate amino acid sequences to Nucleotide sequences

I want to convert a list of fasta ( protein sequences) in a .text file into corresponding nucleotide sequences. A Google search gives me result of DNA to protein conversion but not vice versa. Also, I ...
1 vote
0 answers
27 views

How can I get the estimate of lmFit function form limma package

This is the first time that I used lmfit function from the limma package I'm a little bit confused on how to interpret the result. is there a way to get the ...
0 votes
1 answer
150 views

How to concatenate two fasta sequences by py

I have three backbone and i want to concatenate 70 sequence into these backbone. such like: fasta file 1: ...
1 vote
2 answers
49 views

Multiple Alignment cost application

Is there any biological application/case where someone would be interest in the estimated total cost of the alignment between a set of sequences (genes or amino acids) without the aligned sequence ...
20 votes
8 answers
9k views

What is the fastest way to get the reverse complement of a DNA sequence in python?

I am writing a python script that requires a reverse complement function to be called on DNA strings of length 1 through around length 30. Line profiling programs indicate that my functions spend a ...
1 vote
1 answer
51 views

Get raw genome sequence

I'm currently working on testing some matching algorithms for strings. I would like to do some tests on raw genomic data as I expect different results from random strings given the lower entropy. How ...
0 votes
1 answer
153 views

Reading fasta sequence regions with faidx

Im coding in c++ and reading in different reference genomes to examine regions across the chromosomes a few hundred basepairs at a time. To do it in c++ i use the library htslib and the command ...
2 votes
1 answer
57 views

Mining documents for nucleotide and peptide sequences

The text mining literature has an emphasis on identifying and normalizing gene names, mutations, pathways, concepts, and so on. I haven't been able to find much, however, on methods for extracting ...
0 votes
1 answer
79 views

Help Adjusting Bedtools Output in R

I have been asked to write a script in R to edit, using data.frame, this .fa file that is the result from extracting gene sequences from a repeat masker coordinate file on the genome I am working with....
1 vote
1 answer
30 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 ...
0 votes
1 answer
50 views

BioJava - Protein to DNA

Using Biojava 5, I am trying to convert a sequence of Amino Acids to a Nucleotide sequence. Any idea how to do this?
6 votes
1 answer
374 views

Biohackers Netflix - DNA to binary and video

I'm not sponsored or anything, just interested in their challenge to decipher their DNA code. They encoded their first episode of "Biohackers" video/binary file to DNA code and said if we ...