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.
263
questions
0
votes
0
answers
34
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 ...
- 143
1
vote
1
answer
40
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,...
2
votes
0
answers
25
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 ...
- 121
3
votes
2
answers
67
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
204
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 ...
- 137
1
vote
0
answers
35
views
Most reads no_feature with htseq-count
I have mapped reads to the yeast (R64-1-1) genome with STAR like this:
...
- 281
2
votes
1
answer
28
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
40
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 ...
- 307
4
votes
3
answers
102
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-...
- 243
2
votes
1
answer
28
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
72
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
154
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 ...
- 41
3
votes
2
answers
58
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
1
vote
1
answer
68
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?
- 41
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
1
vote
0
answers
27
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 ...
- 21
1
vote
0
answers
28
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
201
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 ...
1
vote
1
answer
31
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 ...
- 31
0
votes
0
answers
55
views
Extract SNPs from multiple sequence alignment
I’m wondering if folks have recommendations for tools/scripts to extract SNP sites from a multiple sequence alignment of consensus bacterial genomes? Specifically, I am interested in a multiple ...
- 11
1
vote
2
answers
58
views
Phylogenetic tree building from proGenomes database for shotgun metagenomics
For some weeks I'm fighting with an issue about phylogenetic tree building to use in a phyloseq object in order to calculate beta-diversity metrics that takes into account phylogenetic distance.
I’m ...
- 41
0
votes
1
answer
53
views
Why must a maximal non-branching path be a contig?
The following is from Bioinformatics Algorithms:
Fortunately, we can derive contigs from the de Bruijn graph. A path in a graph is called non-branching if in(v) = out(v) = 1 for each intermediate ...
- 127
2
votes
1
answer
159
views
What is the different between the I-TASSER, phyre2, SWISS-model in the 3D tertiary structure?
What is the difference between the I-TASSER, phyre2, and SWISS-model in the 3D tertiary structure?
How do they get the results? When I did a prediction, I got a similar result for the highest template ...
2
votes
0
answers
82
views
How to convert Foldseek .m8 alignment files to FASTA format?
I want to use sequence alignments found by Foldseek in PROSS. I know there are other ways to get alignments (also tried with results of HHPred), but I had in mind that Foldseek may find remote ...
- 401
1
vote
1
answer
126
views
Bowtie2 mapping with multiple indexes
Background:
We performed NGS using cells collected from mice in a xenotransplantation study.
As such, the FASTQ files contain reads of DNA from both mice and human cells.
I expect ~30% of reads are ...
- 97
6
votes
1
answer
78
views
Where to find asymmetric nucleotide substitution matrix with IUPAC encodings?
I am using the pairwiseAlignment function from the Biostrings package to calculate the distance between a consensus sequence and ...
- 483
2
votes
1
answer
118
views
Sequence alignment with MUSCLE
For protein multiple sequence alignment (MSA), I am using MUSCLE (v5.1). The default setting without any specifications gives a 'bad' output even though my proteins are very similar (UniProt sequence ...
- 75
0
votes
0
answers
26
views
Help me to calculate Heaps Alpha value from the roary pangenome pipeline result?
I need to know whether my pan-genome is open or closed. For that, I need to calculate the ...
- 109
4
votes
1
answer
178
views
BWA mem | samtools view: Intermittent parsing error
Update
The issue was that bwa was running out of memory and failing, but that error wasn't floating to the top (see @Steve's answer, below). I was getting an error ...
- 1,254
4
votes
1
answer
133
views
Programmatic way of accessing the FASTA sequence similarity tool (or similar) in Python
I am looking for a tool that performs a sequence similarity algorithm for proteins. More specifically, I am looking for something that would be usable in Python (or anything else usable in the command ...
- 405
2
votes
1
answer
88
views
Samtools sort: most efficient approach to sort a single sample after aligning split .fastq files
Related to my other question (Samtools sort: most efficient memory and thread settings for many samples on a cluster), we need to optimize samtools sort as we ...
- 1,254
2
votes
2
answers
569
views
Samtools sort: most efficient memory and thread settings for many samples on a cluster
We're preparing to analyze thousands of .bam files beginning with re-alignment, sorting, etc.
Complicated question: Has anyone investigated the optimal thread and ...
- 1,254
0
votes
1
answer
25
views
Are there different PAM matrices for different protein families? - How is the inferred common ancestor known to start with?
I am confused about the sequences one starts with to generate the PAM matrices.
The sequences to start with are supposed to diverge by no more than 1% (in sequence identity) wrt a common "...
3
votes
0
answers
45
views
PSSM Matrix in PSI-BLAST
After running PSI-BLAST one obtains a profile matrix with this header:
A R N D C Q E G H I L K M F P S T W Y V
How does PSI-...
- 257
-1
votes
2
answers
68
views
ClustalOmega alignment
I am performing multiple sequence alignment using Clustal Omega for 600 RNA covid-sequences, i.e. cDNA on Genbank, with input characters around 30574 characters for each sequence.
I am running it on ...
1
vote
2
answers
105
views
What is the most appropriate way to find the most recent common ancestor between two distantly related species
I want to specifically find the common ancestor between a lobster and a humans. I suspect it was an aquatic worm of some description. But I want to know about the nervous system of this common ...
1
vote
2
answers
128
views
Help with MinION sequencing data species identification
Hi I'm new to bioinformatics and have just completed my first run on the MinION (long read sequencing Oxford Nanopore Technologies). I was hoping someone could direct me towards R packages, workflow, ...
- 13
0
votes
0
answers
37
views
Why does Nucmer Genome alignment produce faulty Dot plot?
I am using Mummer v.3.23 to make a dot plot between two genomes of the same species. One genome was assembled using PacBio and one from Dovetail. None of them is ...
- 75
2
votes
1
answer
95
views
MarkDuplicatesSpark failing with cryptic error message. MarkDuplicates succeeds
I have been trying to follow the GATK Best Practice Workflow for 'Data pre-processing for variant discovery' (https://gatk.broadinstitute.org/hc/en-us/articles/360035535912).
This has all been run on ...
- 21
0
votes
1
answer
51
views
Extracting Sequence Unique to a Certain Genome
I have two genomes with high similarity. But I found out that one genome has a longer sequence of about 200 kbp. I try to align the sequence with Mauve. How can I extract this unique 200 kbp?
- 121
1
vote
1
answer
32
views
replacing periods with hyphen
I am a software engineer and not a bioinformatician. I am looking at some code where they are looking at aligned proteins that are in the A2M alignment format (https://compbio.soe.ucsc.edu/a2m-desc....
- 113
1
vote
4
answers
300
views
Remove gaps from alignment?
I have an MSA (protein sequence) and I have used various programmes (Clustal, Aliview, MEGA11 etc) to align. However, in all programmes I get many gaps which is not ideal as I am trying to construct a ...
- 85
0
votes
2
answers
197
views
How to compare sequences of genes obtained after whole-genome sequencing to reference genome?
My idea is to align whole-genome sequencing data (as fastq files after, 30× coverage, gDNA) to the reference mouse genome (NCBI), extract the immunoglobulin locus and compare it to the reference. I ...
- 101
0
votes
0
answers
17
views
Global vs Local alignment scoring matrix
When creating scoring matrices for global and local alignments is there are difference regarding the "highest score of the matrix". For example, for the two matrices below, the correct ...
- 103
0
votes
0
answers
217
views
pysam "Exec format" error
I am a beginner and trying to read a bam file in Python.
The lines below throw the error
...
- 11
2
votes
0
answers
52
views
1
vote
0
answers
26
views
Colouring amino acid substitutions by scoring matrix
I have a list of site-specific amino acid substitutions from a mutagenesis experiment which I am visualising by colouring the structure of the protein in PyMOL. I would like a way to colour the ...
- 347
0
votes
0
answers
18
views
Type and number of gaps in sequence alignments
I used Kalign and Muscle to align given sequences, the results are presented below. My question is, based on the two alignments which tool seems better? I would personally believe that the Muscle tool ...
- 103
2
votes
1
answer
33
views
blastn returning an inferior alignment
In my work I've found that in circumstances where there is a mismatch near the edges of the query sequence, blastn prefers to return a shorter contiguous alignment, rather than allowing for a mismatch ...
- 249
1
vote
1
answer
81
views
Given two ADN sequence, can there be multiple optimal alignement? And how to optimally find them?
Given two sequences, I know we can determine optimal alignment of sequence with the Needleman–Wunsch algorithm. But is it possible that two alignments have multiple optimal alignment (I guess it would ...
- 11