Questions tagged [homework]
Use for questions that are associated with canned problems, where either the solution is typically known in advance, or data are simulated to create a specific situation.
44 questions
2
votes
1
answer
68
views
Randomized Motif Search - Bioinformatics Algorithms Week 2 - Cogniterra
I’m struggling to understand this problem in Week 2 of Cogniterra’s Bioinformatics Algorithms course. It is the Randomized Motif Search problem, which prompts the user with the following pseudocode:
<...
3
votes
2
answers
82
views
Multi-pattern search in aligned sequences
I am currently working on a bioinformatics problem where I need to lookup and count the location and count of occurences of 4000-ish 5 character long patterns in each sequence of a fasta file of 700GB....
5
votes
2
answers
95
views
Help Understanding the Greedy Motif Search Algorithm from Textbook
I am trying to go through this text book on Bioinformatics: https://www.bioinformaticsalgorithms.org. It has coding challenges and problems which I hadn't had too much trouble with until I got to the ...
2
votes
0
answers
51
views
Using Python for Monte Carlo simulation for predicting protein nucleation aggregation
I want to write a code for a Monte Carlo simulation to determine the maximum stable complexes in a protein nucleation reaction. Also, need to graph the results for the same.
Here is the problem ...
4
votes
2
answers
114
views
probability of finding a 5 amino acids in a row within a proteome
How to calculate the probability of finding two proteins that share a 5 amino acid long motif from a proteome of around 1067 proteins that have an average length of 65 residues.
The probability of a ...
3
votes
1
answer
317
views
1.What size PCR product will be generated using the primers (bold and underlined) in the sequence below?
I would like to know the steps of this question and how to determine the forward and the reverse primers and also the length of the sequence, without using NCBI. If you could explain and give an ...
-1
votes
3
answers
2k
views
Find open reading frames in a DNA sequence
I have a fast file containing DNA sequences. I've already cleaned the file and has only DNA sequences with no gaps or errors or white space characters. I want to find open reading frames in the fasta ...
1
vote
1
answer
62
views
SNPs location out side the gene
I downloaded exome data from gnomad to get the allele frequency of some SNPs that I'm analyzing. The SNPs are from external resource. Some SNPs (in my data) are not present in the exome data (in ...
0
votes
0
answers
134
views
Help sought with GWAS and vcf files, lack phenotype labels
This question has also been asked on Biostars
Hi, I am very new to this area, and I am taking a class about bioinformatics. For an independent project assignment, I need to do a GWAS. I am using the ...
1
vote
1
answer
117
views
finding the frequency of a motif-matching peptide
This question has also been asked on Biostars
I want to find the frequency of each motif-matching peptide.
The sequences are:
...
1
vote
1
answer
103
views
DESeq2 and EdgeR
I am new to using both DESeq2 and EdgeR in Bioconductor used for transforming my RNA expression data.
However, I am struggling to understand their specific purpose, differences between them and ...
1
vote
0
answers
111
views
When does a 2D histogram / heatmap take a circular shape?
The above is a 2D histogram of y and z internal coordinates describing the positions of C-alpha atoms involved in hydrogen bonds between antiparallel Beta-strands.
There is no color coding. It is ...
0
votes
0
answers
101
views
Relationship between molecular weight and linear peptide
Sorry for this very beginner homework question:
Which of the following linear peptides is consistent with Spectrum = {0 71 99 101 103 128 129 199 200 204 227 230 231 298 303 328 330 332 333}? (Select ...
1
vote
0
answers
614
views
DNA genome string reconstruction from k-mer
I have the following quiz question, but the Pattern1 for both (ACC|ATA) and (CGA|ACT) are unique (just grep for ...
2
votes
1
answer
100
views
how to create "sample file" for the qAlign() function after trimming the reads in R
I'm an absolute beginner trying to solve this question "Align the trimmed and untrimmed reads using QuasR and plot alignment statistics, did the trimming improve alignments?"
I did trim the ...
0
votes
0
answers
44
views
Problem of writing scilab code on a project - Morris-Lecar model
This is a computational neuroscience project. We are asked to plot the Morris-Lescar model (question 3 below).
The model should be plotted with scilab or R.
however, I don't know where to start? For ...
0
votes
1
answer
68
views
How can I draw a diagram of hydrogen bonds only on the basis of C-alpha backbone?
We are working on an algorithm to assign proteins only on the basis of alpha-carbons.
My teacher told me to draw a diagram to represent hydrogen bonds only on the basis of alpha-carbons.
I showed him ...
7
votes
2
answers
620
views
Is there any function in Biopython to convert a DNA sequence from ambiguous to unambiguous?
I have a project to write in Python that requires to write a function that given degenerate DNA sequence (for example: KKGTACACCAG) sequence and a molecular weight interval, returns a list of all ...
1
vote
1
answer
143
views
Generating sequences from HMM
I encountered this question in a textbook (source: https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.721.2540&rep=rep1&type=pdf) and was trying to workout the solution for it. Was ...
1
vote
1
answer
156
views
How to reconstruct a target string from a large set of overlapping reads
I got this question as a homework. Does anyone know how the it can be solved?
Consider a short string $r$ as a read from a string $S$ if $r$ is converted from a substring of $S$ by at most $k$ swaps ...
0
votes
1
answer
75
views
Which sequencing technologies are considered short read
I couldn't find any information on this online. Which of the following sequencing technologies are considered short read technologies? Illumina, Oxford Nanopore, Ion Torrent, Pacific Biosciences
2
votes
2
answers
574
views
Construct the Overlap Graph of a Collection of k-mers
I am doing a course on graph algorithms in genome sequencing and the current assignment involves building an overlap graph from a bunch of k-mers. Formally the problem is as follows:
Input: A ...
1
vote
1
answer
493
views
Create GFF3 feature exporter
Problem Statement
The GFF3 format is a commonly-used one in bioinformatics for
representing sequence annotation. You can find the specification
here:
http://www.sequenceontology.org/gff3.shtml
The ...
3
votes
2
answers
422
views
Sequence allignment with suffix array?
I'm writing the code from scratch and not using libraries for the actual indexing/search since it's a project for school. Any blackbox advice is fairly useless.
I'm planning on aligning some short ...
1
vote
1
answer
393
views
Find species from FASTA files
In a school project we got the task of finding the species from two FASTA files. We got some hints of what to look for, like it is a unicellular eukaryote and that we can use blast. We don't know if ...
0
votes
1
answer
58
views
How to determine a correct promoter sequence? (Bisulfite sequencing)
I'd like to ask you for help with an university assignment. First let me post the entire question:
Design primers for bisulfite sequencing of the promoter of the ACE2
gene in microbat Myotis ...
1
vote
1
answer
102
views
BioPython - Retrieve sequence records from pubmed database
I've been assigned a task to fetch sequence records from "pubmed" database:
fetch 30 record of type "fasta" from (pubmed database, with term hemoglobin AND alpha) as a sequence ...
1
vote
0
answers
55
views
Comparison of phylogeny of SARS2 whole genome vs spike genome
I’ve built two phylogeny at minimum evolution in MEGA-X and extracted these as newick files.
My next step is to compare the two phylogeny to look at if the mutations within SARS2 are more prevalent ...
2
votes
4
answers
382
views
I have really skewed RNA-seq data, what's the best way to normalise it? Preferably in R!
My data frame compares the RNA-seq reads from many genes in different tissues.
The reads look as follows.
I tried using log to make it better but still looks pretty skewed.
Are there any better ...
1
vote
2
answers
80
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 ...
2
votes
1
answer
5k
views
how to calculate Pairwise Alignment Scores for blosum62
I have a bioinformatics exam coming up. I can understand the difference between the global algorithm and the local algorithm, but I have a problem with gap opening penalty and gap extension penalty. ...
2
votes
1
answer
239
views
Coverage required
I was came across a problem during an exercise in a book and I don't really know how to solve it. I feel like something's missing.
"coverage, c = $NL/G$ (N=number of reads, L=read length, G=genome ...
0
votes
1
answer
50
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 ...
1
vote
3
answers
2k
views
Calculate the occurrence of motif in sequences (including overlaps) using Python
I need to calculate the occurrences of a motif (including overlaps) in sequences (motif is passed in the first line of standard input and FASTA sequences in subsequent lines). The sequence name starts ...
3
votes
0
answers
193
views
How to assign the best gap penalty and gap extension penalty using BLOSUM65 [closed]
For an assignment I must do a pairwise optimal local alignment using BLOSUM65 and five protein sequences. The algorithm I want to use is the Smith-Waterman.
Context protein sequencing using Blastp: ...
2
votes
1
answer
632
views
Generating basic dna sequences in R
Looking for some guidance.
In an introductory bioinformatics class I am tasked with doing the following in R:
Create the sequence ACTGACTG
Generate a random string of nucleotides with equal ...
4
votes
1
answer
865
views
How to map short sequences to long reads, recovering all multiply-mapped high-quality matches
The dilemma:
I have around one million short sequences (21 bp to several 100s of basepairs) for which I need to identify all occurrences of in 20-30x coverage noisy long reads (both pacbio and ONT). ...
4
votes
2
answers
160
views
Writing a perl script to holding information for two genes
Basically I have a perl script in which I have an array (where each element of the array references a hash) and need to be able to print the array with a dumper function.
Thus I need to be able to ...
0
votes
1
answer
87
views
Find number of possibilities [closed]
We have a column which has following format : (1, 1, 1, 0, 1).
Here is the definition of what means two columns are compatible.
Using the notation Oi to denote the collection of rows possessing a ...
4
votes
1
answer
1k
views
Inspection of gene expression in scRNA-seq data
I am running the data preprocessing pipeline for scRNA-seq data presented here.
3.8.6.1 Gene expression
In addition to removing cells with poor quality, it is usually a good idea to exclude genes ...
4
votes
2
answers
821
views
Finding a single open reading frame with ribosomal binding site, using Biopython
I'm given a Fasta file, containing a large DNA(over 115,000 long) sequence, and I am tasked with finding a single large open reading frame contained within the DNA sequence using Biopython.
I'm ...
1
vote
1
answer
61
views
Finding out proteins that have a similar domain to a specific toxin
I am trying to find out that which DNA sequences(in the below) have a similar domain to a specific toxin which is Clostridium botulinum
You can find the fragment of the DNA sequence that belong to ...
2
votes
1
answer
454
views
How can I calculate the silent mutation for each position in codon?
I am trying to find the number of silent mutations in for each position in codon.
In the homework, I am asked:
For the codons and for each position, count the number of silent mutation. There are ...
5
votes
4
answers
2k
views
Range overlap python error with genomic regions
I have two files
s3.txt :
1 10 20
1 5 20
2 20 30
2 25 30
1 10 50
2 20 60
1 14 17
s4.txt:
1 10 20
2 20 30
I am trying to match col0 ...