Questions tagged [sequence-alignment]
Use this tag to refer to the alignment of nucleotide or peptide sequences, producing an arrangement that maximizes similarity between sequences.
410
questions
0
votes
1
answer
32
views
Mapping List of Probes/Primers/Short Oligos to a Reference Fasta/q
Could you help me bioinformatics SE people.
So I have been duckduckgo-ing up a storm and even consulted the AI overlords and haven't come up with a solid way to do this, but it is so fundamental to ...
- 51
0
votes
0
answers
15
views
Error while opening xmfa generated by progressiveMauve
I used commandline progressiveMauve to align two genome .gbk files. The alignment was successful. But when I try to visualize the final .xmfa file with Mauve, I get the following error.
...
- 105
1
vote
1
answer
50
views
Pal2nal translation of large multi-fasta files produces a codon translated file where some sequences half length of the average
I did sequence alignment of a large peptide multi-fasta (n= 4991 sequences). The peptide alignment has sequences with the same length and pal2nal went through just fine... except some of the codon ...
- 217
1
vote
0
answers
63
views
Pacbio HIFI pbmm2 alignment metrics
I am new to pacbio sequencing data. I just did some alignment of pacbio HiFI data using pbmm2. I have the bam now, I would like to collect some metrics on the alignment. I used to work on the illumina ...
- 119
1
vote
1
answer
100
views
How to add bootstrap values to the phylogenetic tree generated by OrthoFinder?
When we run the OrthoFinder analysis tool on a group of genomes to get the orthologues shared by them one of the output files include a folder named 'Species_Tree' that contains a text file named '...
- 53
1
vote
0
answers
116
views
Mapping statistics from the bam files
I would like to find the mapping statistics from the sorted bam files.
Samtools flagstat gives the output only for a single file. What is the easiest way to find the mapping statistics for all ...
- 11
0
votes
0
answers
13
views
Automated method to figure out whether a recognition site is conserved or not using the results of simulated restriction digestion
Restrict from EMBOSS Suit version 6.3.1 with the following
parameters: snucleotide1, sitelen = 4, rformat = table, enzymes =
enzymes.txt.
From the 4379 enzymes present in REBASE, we selected the 650 ...
0
votes
0
answers
18
views
Seeking Guidance on Designing DNA Sequences with Specific Criteria
Dear Bioinformatics Community,
I hope you're all doing well. I am currently working on a project that requires designing a DNA strand or more, every 10 nucleotides long (preferably), with specific ...
- 1
2
votes
2
answers
46
views
RNAseq alignment: best practices for aligning to multiple isoforms?
I have Illumina RNAseq data and would like to maximize my power to find candidate genes that are differentially expressed genes between experimental conditions.
Many of my (de novo assembled and ...
- 23
0
votes
1
answer
61
views
Is my reference sequence too small?
I'm trying to map ONT long reads to a portion of a gene I'm looking at. The region is about 25bp long. When I search for the region in the document it pulls up the sequence in every read but when I ...
- 844
3
votes
0
answers
23
views
Add or Simulate Noise in Short-Read Paired End Data
For a project I'm working on, I need to figure out how to model noise that may be happening in real genomes due to alignment errors, contamination, etc. Specifically, short-read paired-end data either ...
Suraj Rajendran
3
votes
1
answer
62
views
Sanger Sequencing Knitting Error
I am doing a project where I am reproducing the analysis from the article "sangeranalyseR: Simple and Interactive Processing of Sanger Sequencing Data in R". Below is the example chunk for ...
- 31
1
vote
1
answer
74
views
Phylogeny building in R from FASTA files:
Im working on building a phylogeny from scratch with downloaded FASTA sequences from GeneBank. I think Im doing alright up until the multi sequence alignment in the ...
- 111
0
votes
1
answer
67
views
How do I build the MinusB database for Kraken2? (Taxonomy issues)
I am attempting to build my own custom database for Kraken2. I have two questions:
If I have the MCPyV genome in a file called MCPyV.fasta, how do I build a database with just this?
How do I build ...
3
votes
2
answers
64
views
How is the "canonical" version (`_entity_poly.pdbx_seq_one_letter_code`) obtained in the PDB?
I encountered the dichotomy in the context of PDBx/mmcif files, say, 6OSQ: each chain has a
_entity_poly.pdbx_seq_one_letter_code...
- 342
3
votes
3
answers
193
views
Read Clustal file in Python
This question was also asked on StackOverflow
I have a multiple sequence alignment (MSA) file derived from mafft in clustal format which I want to import into Python and save into a PDF file. I need ...
1
vote
1
answer
64
views
Maximum likelihood estimation for tree construction
Could someone direct me to a resource (textbook chapter, lecture notes, online video, etc.) that demonstrates, in a mathematically or computationally rigorous way, how to use maximum likelihood ...
1
vote
1
answer
106
views
Matches, mismatches and indels
I have a bam-file of reads and reference genome. How can I located matches, mismatches and indels in three different files with Phred score in each data of the files?
- 21
3
votes
1
answer
146
views
heatmap of specific sequence motifs in aligned fasta files
I have a collection of fasta files, each containing three aligned sequences. I am interested in understanding the distribution of a specific sequence motif PGP, RGP and KGP in all the alignments. I ...
- 31
1
vote
2
answers
110
views
How can I align two novel sequences to a reference genome and build a phylogenetic tree?
I have two new amino acid gene sequences that I want to align to a reference genome. Is it a good idea to have a combined phylogenetic tree? what will this type phylogenetic analysis called?
- 11
1
vote
1
answer
81
views
How to identify sequence origin in DNA shuffle reads?
The data
I have a large number of reads from sequences that were generated by randomly shuffling regions of two parent sequences together. See the following image:
The regions that are shuffled are ...
- 13
2
votes
1
answer
247
views
Why does the Smith-Waterman algorithm have a complexity of $O(m^2 \times n)$?
I recently read that original the Smith-Waterman algorithm to do a local pairwise alignment of sequences has a complexity of $O(m^2 \times n)$, where $m$ and $n$ would be the lengths of the sequences ...
3
votes
2
answers
150
views
What can be the bias of aligning paired-reads in a single-end mode?
I don't know if I'm in the right place but I have a technical problem to fix. I would have to align paired-end reads from Illumina sequencing to compare a normal genome with a tumor one.
When I align ...
- 43
3
votes
0
answers
386
views
Aligning PacBio HiFi reads to reference genome using pbmm2
I am trying to align a yeast strain sequenced by PacBio HiFi reads to the reference genome S288C (https://www.ncbi.nlm.nih.gov/data-hub/genome/GCF_000146045.2/) using pbmm2 but for some reason I am ...
- 844
2
votes
1
answer
250
views
How to calculate average BLOSUM62 scores?
I can understand the motive behind the BLOSUM62 matrix, this being a pairwise mutation matrix describing aggregate mutations between the 20 amino acids.
However how would you calculate the average ...
0
votes
0
answers
66
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 ...
- 141
1
vote
1
answer
91
views
How to chop fasta / bed /peak file genomic segments into smaller fixed or custom genomic intervals
I am trying to split my .bed/.fa file genomic segments into arbitrary smaller overlapping intervals:
Consider a typical line/row in my *.bed file as follows:
chrom. $\ \ \ \ \ $ start $\ \ \ \ \ \ \ \...
- 323
2
votes
1
answer
54
views
What is the best way to process yeast genomes?
I have obtained several hundred raw, unassembled yeast genomes from NCBI and I am looking for advice on how to process the genomes for downstream analysis.
I have a reference genome (S288C) to use for ...
- 844
1
vote
0
answers
20
views
The difference between sequence-independent and sequence-order independent structural alignment?
I am confused about the difference between sequence-independent and sequence-order independent structural alignment. TM-align for example claims to be sequence-independent; however, I've seen it is ...
- 11
2
votes
2
answers
1k
views
How to convert fasta alignment to nexus or phylip format?
I generated a multiple sequence alignment using mafft, but I would like to convert it to nexus or phylip format (required by a downstream software). What's the easiest way to do the conversion?
- 5,497
2
votes
1
answer
478
views
PyMol alignment script
Is there a way write a script in python to use "super" to align proteins from a 2D array? Each row of the array has 2 columns, and for each row, I want to "super" the two proteins (...
- 21
1
vote
1
answer
67
views
Lower case vs. upper case nucleotids in sequence vs. dots at the end
What is the difference between lower case and upper case nucleotides in a sequence? My other question is what are the dots at the end of the sequence? Some examples are shown below:
GGgG,GGgG,GGgG,...
3
votes
0
answers
35
views
How are BLAST record scores calculated in biopython
Can someone explain how the score is calculated for Bio.Blast.Record.Description? For example, if Record.Description.num_alignments is >1, is ...
- 131
3
votes
2
answers
135
views
How can I find the maximum percent identity between two sets of (unaligned) sequences?
As the title states, I've got two sets of unaligned amino acid sequences (~25k sequences in one, ~3k in the other). I want to find the minimum distance between a sequence in the first set and a ...
- 163
3
votes
2
answers
2k
views
STAR vs Bowtie2
What is the fundamental difference between STAR and Bowtie(2). Specifically, what is the difference in their final output (regardless of run-time differences, speed, memory usage etc.). Both seem to ...
- 323
1
vote
0
answers
210
views
Most reads no_feature with htseq-count
I have mapped reads to the yeast (R64-1-1) genome with STAR like this:
...
- 335
1
vote
1
answer
33
views
Multifasta protein alignment
I have around 4000 orthologous proteins for each of 5 species, I'd like to align them and to concatenate alignments. Which tools can I use? And how?
Marco
- 151
2
votes
1
answer
40
views
Statistical significance of pairwise protein sequence alignment
I have identified a conserved region in a protein I'm studying, by doing multiple sequence alignments against a dozen of its orthologs. I've then found that this region also aligns relatively well ...
- 183
1
vote
1
answer
128
views
Aligning scRNA-seq fastq to .bam without cell barcodes
I would like to convert the fastq files from this dataset: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE98679 to .bam files. The group did not provide the scripts they used to to their ...
- 376
1
vote
1
answer
38
views
Problem with sequence alignment
I'm trying to find if the peptides in column A are present within the longer peptides of column B.
Pairwise alignment tools take the data from both columns as one string (i.e., TIIDYTNNHYTNNHLEII... ...
6
votes
4
answers
636
views
multi-sequence alignment of samples with multiple contigs each
I want to perform a multi-sequence alignment on 12 samples that clustered based on cgMLST. Ultimate goal is to find out whether they differ by the presence of a few genes.
I have performed multi-...
- 283
2
votes
1
answer
47
views
What tool I can use to map short-reads sequences to reference genome and get specific mapped size
I have about 300 90-bp sequences which I would like to map to a reference genome to make it longer to a 300bp sequence, wherein the 90bp is at the middle.
Anyone knows what bioinformatics tool I can ...
- 21
2
votes
2
answers
158
views
BWA-mem and sambamba read group line error
this question has been asked [and answered] on Stack Overflow
This is a two-part question:
help interpreting an error;
help with coding.
I'm trying to run bwa-mem ...
3
votes
1
answer
900
views
How to remove the unpaired reads in sam/bam files?
I have sam and bam files for the chimeric reads, which come from two different parts of the genome (For example, the first half of the read from part of Chromosome 1 and the second half of the read ...
- 101
3
votes
2
answers
78
views
How to align integer (non-DNA/protein) sequences?
I am looking for an algorithm to find the "best" alignment between two sequences of integers similar to how one aligns nucleic acids or amino acids for homology comparisons. For example, the ...
- 131
2
votes
1
answer
257
views
How does split reads look like in sam files?
I used bwa mem to align the DNA with the reference genome. If there are split reads (chimeric reads, come from two different parts of the genome), will it be split into two lines rather than one line?
- 101
0
votes
0
answers
24
views
How to identify most important Amino Acid residues in a given protein?
Continuing from the title, in my context, important means that changing the set of AAs (either or one more many) changes the domain and/or fold and/or function and/or family etc.
I noticed that there ...
- 1
3
votes
1
answer
125
views
Identification of unique chains in proteins by sequence similarity
I'm trying to find unique chains or proteins within their PDB files. Many proteins have multiple chains, but very often they are identical (say, the PDB file consists of 4 homodimers, for example, so ...
- 43
1
vote
0
answers
42
views
Multiple sequence alignment with gaps
I'm working on a problem where the elements of multiple arrays need to be aligned with each other. While doing some research on algorithms that could possibly be applicable to this task, I stumbled on ...
- 11
3
votes
1
answer
865
views
Perform protein structure-based sequence alignment in Python
I am looking for a Python package that performs pairwise structural alignment of protein structures (i.e., PDB files) and returns a sequence alignment. PyMOL is able to do this through the GUI, for ...
