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
2answers
43 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 ...
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
1answer
70 views

How to calculate Centimorgan lengths from run lengths in one-to-one dna profile comparison

For easier management of a database of dna profiles, I'm creating a tool to perform one-to-one matching on two profiles with snp genotype information. Basically, something that works similar to ...
1
vote
1answer
51 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
1answer
27 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
2answers
115 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
vote
2answers
1k views

Why do NEB adapters have non-complementary sequence?

I am making some infographics of different library prep types and noticed something weird about NEB's Ultra II adapters. Molecule 1 is the desired product of the ...
1
vote
1answer
19 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 ...
1
vote
2answers
51 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 ...
1
vote
1answer
17 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 ...
1
vote
1answer
38 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: ...
1
vote
1answer
61 views

How can I convert codon coordinates to genomic position?

I am looking for a given mutation in IGV, which accepts coordinates in the form of chr<X>:<Y> Where X is the chromosome number and Y is the base ...
1
vote
1answer
48 views

Why do we scale the rate matrix of a substitution model to make the average rate of substitution equal 1?

The matrix of transition probabilities for a substitution model Q over time t is found as follows: $e^{Qt}$. Since Q and t only show up as a product, we cannot distinguish time and rate of change. So ...
1
vote
1answer
1k views

What i5 index should I use on the Illumina sample sheet for an unindexed p5 primer?

I have an upcoming run on a HiSeq X and most of the libraries in the pool have both i5 and i7 indices. However, some of the libraries were made with the IS4 p5 oligo and it is unindexed. The IS4 ...
1
vote
1answer
30 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: ...
1
vote
0answers
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 ...
1
vote
0answers
37 views

"How the clustal omega can be reverse engineered, to trace ancestral inversion mutations via the guide tree?"

I apologise for very basic question but I am very new to biology and have very specific knowledge of this field, I am currently working on Bioinformatics in my machine learning project which is why I ...
1
vote
0answers
59 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
vote
0answers
20 views

mining DNA barcodes from Genebank/BOLD per location

Is there a method to download all the sequences for a particular geographic region (Panama) for 'DNA barcode' sequences (ITS2, rbcL or trnL). Hopefully, the specimen collection location will include ...
1
vote
0answers
40 views

Is there a software that can predict protein-DNA binding?

Is there any software that, given a DNA sequence as input, can predict what protein/ type of protein can bind to the DNA. The sequence of interest is a potential transcriptional regulator, I want to ...
1
vote
1answer
91 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 ...
1
vote
0answers
41 views

Unmapping Alternative reads in BAM file

I want to do some HLA typing, most of the tools require the bam file is aligned to primary genome without alternative read handling. From the International Genome Sample Resource project, I can get ...
1
vote
0answers
25 views

What are the meanings of these transcript ids?

I have three types of transcripts for Rosa chinensis "Old Blush" homozygous genome v2.0 from GDR genome browser. ...
1
vote
0answers
19 views

How can I find genes of fundamental metabolic processes for an organism?

I want to find genes of fundamental metabolic processes in the Phytophthora cinnamomi genome. I'm using the genome sequence deposited in the NCBI and I'm searching ORF by ORF, using the ORF finder (...
1
vote
0answers
24 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 ...
1
vote
0answers
354 views

Aligning nucleotide sequences in APE software

I am very new to APE software. I am trying to align complementary oligo-DNA strands using APE on macOS. I have both the forward and reverse sequence of an oligo-DNA insert in two separate files. But I ...
0
votes
2answers
102 views

How can I find the relevant pathway map from gene-gene or protein-protein interaction list? [closed]

I have some difficulties to understand/interpret the pathway map and how a gene-gene interaction list or DNA sequencing can map into pathways. In addition what's the difference between MARK/ERK ...
0
votes
1answer
38 views

Encoding a 15bp DNA sequence into a shorter than 15 digits number

I have a series of 15bp DNA sequences from single-cell barcodes and I would like to encode them into a shorter than 15 digits number so that each unique 15bp DNA sequence will create a unique numeric ...
0
votes
1answer
29 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 ...
0
votes
2answers
59 views

How to modify DNA evolution model to fit actual data?

I'm working on understanding the evolution of a gene in a phylogenetic tree. I know the rates of evolution of each organism (from the tree). I am taking a random DNA sequence with my gene, evolving it ...
0
votes
1answer
42 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 ...
0
votes
1answer
87 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: ...
0
votes
1answer
123 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 ...
0
votes
1answer
45 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?
0
votes
1answer
46 views

How can I extract all known mutations of my BAM (or SNP/INDL files)?

I am using a Genome Explorer tool to see all the mutations on my own DNA. My particular interest is on listing the variants and get their names/ids. Here are a few screenshots: You can see that ...
0
votes
1answer
126 views

Get regions' information from DNA sequence data (bsgenome.hsapiens.ucsc.hg19)

I have a problem in R. I have the following dataSet (the first three rows shown) (the 5th number is the methylation level in its region): ...
0
votes
1answer
44 views

Function to find ATP binding sites

I'm currently doing a python project where I need to write code that extracts the ATP binding site from a FASTA file using a regular expression. I've tried the following code but it does not work. It ...
0
votes
1answer
269 views

single stranded cDNA vs double-stranded cDNA

Single stranded cDNA is used by race pcr and rt-pcr while next-gen sequencing (Illumina) requires double stranded cDNA. What is the reason behind this ?
0
votes
1answer
78 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....
0
votes
0answers
13 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 ...
0
votes
1answer
51 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 ...
0
votes
0answers
24 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
0answers
22 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 ...
0
votes
0answers
44 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
1answer
29 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 ...
0
votes
0answers
112 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
0answers
42 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
0answers
23 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 ...
0
votes
0answers
25 views

Greedy Motif Search [ Reasoning ]

I was able to execute the algorithm correctly. But I'm still trying to get a hold of the result I'm achieving. This is the algorithm I'm trying to understand - http://rosalind.info/problems/ba2d/ If ...
0
votes
2answers
50 views

Transcribing DNA exons and append transcripts into a list

I am trying to find a way to read and transcribe a list of DNA sequences (list of lists) only when the for loop finds a start codon (triplet of the list items) and until it finds a stop codon, over ...