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.
255
questions
2
votes
2
answers
51
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
52
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
52
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
45
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
22
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
14
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
24
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
1
vote
1
answer
39
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
25
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 ...
- 21
0
votes
0
answers
29
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
38
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
40
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
73
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
33
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 ...
- 391
1
vote
1
answer
53
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
5
votes
0
answers
41
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 ...
- 443
2
votes
1
answer
71
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 ...
- 45
0
votes
0
answers
20
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
3
votes
1
answer
94
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,234
4
votes
1
answer
76
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 ...
- 385
2
votes
1
answer
69
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,234
2
votes
2
answers
291
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,234
0
votes
1
answer
23
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
24
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-...
- 131
-1
votes
2
answers
61
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
88
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
116
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
20
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
1
vote
0
answers
62
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 ...
- 11
0
votes
1
answer
42
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
29
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
160
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
148
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
14
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
156
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
45
views
1
vote
0
answers
25
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
44
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
2
votes
1
answer
29
views
mpiblast alternatives?
I am looking for a sequence alignment software capable of pairwise alignment of a large number of protein sequences (10^7).
I tried to go for mpiblast: https://wiki.canterbury.ac.nz/display/UCHPC/...
- 249
0
votes
0
answers
38
views
Does the alignment need to be redone when using SAW method to test long branch attraction?
I am looking into different methods to test the long-branch attraction problem. In this chapter they talk about the SAW method - excluding one and the other taxa from the dataset and see if they ...
- 5,287
4
votes
0
answers
61
views
How to get phylogenetic tree from multiple genes?
I constructed a phylogenetic tree using a gene (example - secA). I had to gather the same gene sequence for all the required species from public database-NCBI and then constructed the tree after ...
- 73
0
votes
0
answers
49
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 ...
- 320
-2
votes
1
answer
43
views
Where to find rat genome for STAR alignment?
I'm trying to find the rat genome to carry out STAR alignment for several FASTQ files. Where would I find these files?
-1
votes
1
answer
48
views
Alignment statistics
Given a multi-fasta multiple sequence alignment file, is there a quick way/tool to calculate average pairwise sequence identity?
Even better, given a whole genome alignment file (...
- 758
1
vote
0
answers
97
views
How to modify dot plot in MUMmer 3 for bacteria comparative genomics?
I was trying to make a dotplot to visualize genome-genome sequence alignment by using the MUMmer 3 software (it is typically used to compare whole genome sequences of bacteria). I runned the same ...
2
votes
2
answers
102
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:
...
4
votes
1
answer
73
views
How do you convert Raw Alignment Score to Bit Score?
I'm coding a pipeline where I make a lot of pairwise alignments, and I end up with raw alignment scores. But, I really need to look at my results in terms of bit scores.
I know that the formula is:
$𝑆...
0
votes
1
answer
83
views
Smith-Waterman traceback with large gaps
I am currently trying to understand how the traceback algorithm is supposed to work for the smith-waterman algorithm as my current understanding breaks down in case of a large alignment gap.
Assume ...
- 111