Does anyone know of any existing tools/scripts which can score the quality of MSAs? I have a couple of quite long alignments which I've done with Lastz, MAFFT, Clustalo and Kalign, but its basically impossible to check these by eye.
Crude options are valid too! Just after a metric or some kind of score/comparator I can use to choose the theoretically 'best' alignment.
•
written
by
Joe
1 answer
I think the program MUMSA might be useful for this sort of thing. It has been designed to compare different alignments of the same sets of sequences and identifies regions of consistency between different alignments as well as identifying problematic regions (see ref).
Also in the PRANK paper, they use different methods for comparing alignments.
Log in to answer this question.
-
Tool: New Illumina error mode, new BBTools release (39.09) to deal with it
written by Brian Bushnell
-
How to use Nextflow to call scripts from different environments?
written by O.rka
-
Understanding bowtie2 for short reference sequences
written by Sina •
-
Extracting CDS from fragmented assemblies
written by ebaldwin •
-
Tutorial: Yes .. BBMap can do that! - Part III clumpify (mark (and dedupe) duplicates without align…
written by GenoMax
-
Tutorial: Yes .. BBMap can do that! - Part I : bbmap (aligner), bbduk (scan/trim), repair (fix PE r…
written by GenoMax
-
Bowtie2 local setting & score threshold
written by MTW •
-
Is proper feature selection possible before all imperatively hidden objects required for conceptual…
written by tfhahn •
-
News: Workshop: Advanced Python for biologists
written by Physalia-courses
-
Tutorial: Mapping Reads With Bwa And Bowtie
written by MSU course 2011 •
Perhaps some pairwise distance metric? Levenshtein/Hamming/... averaged,... just thinking out loud
Yeah that's the best idea I've had so far, and I've written a tool in the past to do a similar thing. The only other thing I can think of is the MSA with the lowest average entropy over all columns perhaps. Though I think gaps would need 'special treatment' as a column which was mostly gaps would have low entropy, despite being a fairly uninformative site.
See if this paper helps. Another one here.
I found SuiteMSA, but I don't think it was happy with the alignments I'd given it. It took ages to open them (via Xterm admittedly) and just ground to a halt. I might try it again locally... first paper looks good though thanks!