3
votes
Read Clustal file in Python
pyMSAviz may help with this problem.
pyMSAviz is a Python tool with a convenient CLI/API implementation for visualizing MSA files.
For example, you can output MSA figure in PDF format with the ...
3
votes
Accepted
File format of substitution matrix in clustalw
This example seems to have information on this, looks like it should have this format:
...
2
votes
Multiple Alignment cost application
I'm not a fan personally, but it is context specific. The problem is that it is heading towards optimisation, particularly algorithmic optimisation. If the alignment is the final calculation - that is ...

M__♦
- 10.2k
2
votes
find conservative regions relative to a specific accession
You can opt for a similarity matrix as output (here is an example). The fact that you would like to judge conservation "relative to a given sequence", you can use this matrix to rank ...
2
votes
ClustalOmega alignment
Have you had a look at your alignment? I personally use Jalview but any alignment viewer should work. Colour > Percentage Identity (or Nucleotide) should help with visualization.
Does the alignment ...
2
votes
Select representetive sequence within MSA (protein)
I have two suggestions that I used in the past (both using Python/BioPython though). They are arguably the same, but obtained through different paths.
Calculate the (artificial) consensus sequence of ...
2
votes
multi-sequence alignment of samples with multiple contigs each
in the end, I just loaded the .bam files into artemis and, by inspecting the depth (heat-map), I could check which samples had the genes I was looking for:
2
votes
multi-sequence alignment of samples with multiple contigs each
If you have a reference genome (or are willing to designate one of your de novo assemblies a reference), you may find QUAST helpful.
QUAST will perform an alignment of all genomes against the ...
2
votes
multi-sequence alignment of samples with multiple contigs each
MUMmer4 is a versatile alignment tool for DNA and protein sequences. It supports one reference genome and up to 32 query genomes. MUMmer4 will align every contig in each genome to the reference genome....
2
votes
Accepted
Read Clustal file in Python
Ok, figured out a way, not sure its the best one,
nedd to install fpdf2 (pip install fpdf2)
...
2
votes
Accepted
How to identify sequence origin in DNA shuffle reads?
This is the sort of thing that LAST has been specifically designed to work well for:
LAST is designed for moderately large data (e.g. genomes, DNA reads,
proteomes). It's especially good at finding ...
2
votes
Read Clustal file in Python
adding, a second answer because the OP request about
I need to import the file and then highlight some specific word
was making me uneasy. I kept using the Pyfpdf2...
2
votes
Accepted
How do I plot clusters based on DNA sequence alignment?
I think one of the reasons you struggle is that clustering DNA sequences is not a clearly defined task. In general intention of clustering is to reconstruct, or approximate the relatidness of the DNA ...
2
votes
Remove gaps from alignment?
If you want an ungapped alignment, you could do an ungapped alignment (the older read aligners were ungapped, e.g. bowtie). However, I suggest that you don't want an ungapped alignment, you just want ...
2
votes
Remove gaps from alignment?
You may also find useful the BMGE tool (stands for Block Mapping and Gathering with Entropy). BMGE selects regions in a multiple sequence alignment that are suitable for phylogenetic inference by ...
1
vote
How to highlight the specific peptide sequences after performing multiple sequence alignment for the fasta file?
Checkout an MSA visualisation and analysis package that you can load the alignment into. It may not be the newest tool, but I've used Jalview (https://www.jalview.org/) in the past
1
vote
Remove gaps from alignment?
Gblocks is a tool for removing gaps/poorly conserved positions from alignments. Available here. Here is a webserver that runs it for you.
It's a bit old, possibly there are others at this point.
1
vote
drawing multiple sequence alignment as tree in R
Firstly, do check tree is not an empty object, i.e. it uploads the tree file. Secondly, check the path to the alignment is correct ...

M__♦
- 10.2k
1
vote
multi-sequence alignment of samples with multiple contigs each
It won't work on more than 3 genomes at a time, but SynMap3D from the Comparative Genomics webserver (https://genomevolution.org/coge/) should do the trick. I'm mentioning it here because the ...
1
vote
Generate consensus protein sequence from relatively gappy alignment?
Given the alignment you've presented, there are obvious groups of sequence similarity, but I'm seeing a lot of different peptide sequences in this (e.g. ...
1
vote
Accepted
Sequence alignment with MUSCLE
HMM Super5 is a hidden Markov model (HMM) where gap penalties and the substitution matrix are the parameters at the centre of the model. That is the answer, thats why you can't change it - it does it ...

M__♦
- 10.2k
1
vote
How to make MSA for contact prediction?
You might have some luck with Parsnp. From the description:
Parsnp was designed to align the core genome of hundreds to thousands of bacterial genomes within a few minutes to few hours. Input can be ...
1
vote
How to make MSA for contact prediction?
That's a LOT of sequences! I'm doubtful you'd be able to get TripletRes or CCMPred to run at all with an MSA depth of 1.5 million. I suggest pruning your dataset to remove redundant sequences (...
1
vote
Remove gaps from alignment?
This new tool, CIAlign, talks about some automated removal of gaps in multiple sequence alignments
https://twitter.com/thekatybrown1/status/1306215132670889984
https://github.com/KatyBrown/CIAlign
1
vote
Multiple Sequence Alignment
It is likely the sequence data in this instance, the more sequences the more chance an issue with the alignment where variations in the gap penalty are observed.
MAFFT outperforms Clustal Omega to ...

M__♦
- 10.2k
1
vote
Multiple Sequence Alignment
The length of the alignment will depend on the model. For example, high gap costs will typically result in shorter alignments, whereas low gap-costs will do the opposite.
Consider the sequences below:
...
1
vote
Accepted
Generating sequences from HMM
[Comments from other post migrated as answer]
The parenthetical statement in the first bullet point says without corrections. In which case, if an unseen nucleotide has a probability of zero, the ...
1
vote
BEAUti not accepting taxa tip dates
I found the solution, one of my Seq IDs had a , in one of the meta-data fields, this appears to have thrown off the software and all subsequent Seq IDs could not be ...
1
vote
Median string problem & multiple sequence alignment
Why is this approach not used? Are there any algorithms that use it?
As shown in this paper, the MEDIAN STRING problem is NP-Complete on alphabets of carnality greater than 3. While the problem isn't ...
Only top scored, non community-wiki answers of a minimum length are eligible
Related Tags
multiple-sequence-alignment × 29sequence-alignment × 11
phylogenetics × 5
fasta × 4
phylogeny × 3
sam × 3
r × 2
python × 2
proteins × 2
sequence-analysis × 2
shell × 2
awk × 2
ggtree × 2
rna-seq × 1
bioconductor × 1
gene × 1
protein-structure × 1
genomics × 1
variant-calling × 1
clustering × 1
visualization × 1
linux × 1
reference-genome × 1
wgs × 1
variation × 1