Changing gap open and extension penalties in muscle (or any other command line tool that works on mac, really)
Entering edit mode
3 months ago

I'm trying muscle after learning that clustal omega doesn't allow for changing gap opening and extension scores. I am looking at the documentation, and I don't see a way of changing gap opening and extension scores (except for the termini, etc.). When I run the compiled executable on the command line, it prints out some options, which includes:

 Generate resampled ensemble from existing ensemble by sampling columns with replacement:
muscle -resample ensemble.efa -output resampled.efa

-maxgapfract F
   Maximum fraction of gaps in a column (F=0..1, default 0.5).

-minconf CC
   Minimum column confidence (CC=0..1, default 0.5).

I thus tried

muscle -align seqs.fasta -threads 6 -output test.fasta -maxgapfract 0.1

but I get an error message that says

Invalid command line
Unknown option maxgapfract

I'm surprised that this, in general, isn't a trivial thing to change. Am I using the wrong software? What command line tool I can compile on mac can allow me to change the gap open/extend scores?

Edit: I should also add that I'm trying to align proteins, but I'd like to be able to do this for dna sequences as well.

multiple-sequence-alignment • 260 views

Login before adding your answer.

Traffic: 3696 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6