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.

Filter by
Sorted by
Tagged with
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 ...
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 ...
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 ...
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: ...
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 ...
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 ...
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 ...
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?
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 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 ...
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 ...
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 ...
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 ...
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 ...
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

Bowtie2 gets stuck on alignment

I am aligning a fastq file as follows: ...
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 ...

1
2 3 4 5 6